<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2015.712072</article-id><article-id pub-id-type="publisher-id">ENG-62160</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Thermodynamic Model on Paraffin Wax Deposition Prediction
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>aojun</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Wanting</surname><given-names>Sun</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Chengting</surname><given-names>Liu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Liping</surname><given-names>Guo</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>The College of Petroleum Engineering, Northeast Petroleum University, Daqing, China</addr-line></aff><pub-date pub-type="epub"><day>17</day><month>12</month><year>2015</year></pub-date><volume>07</volume><issue>12</issue><fpage>827</fpage><lpage>832</lpage><history><date date-type="received"><day>4</day>	<month>December</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>21</month>	<year>December</year>	</date><date date-type="accepted"><day>24</day>	<month>December</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In process of crude oil production and transportation, wax and other solid deposition issues have a significant impact on oilfield production. Solid precipitation not only reduces the production efficiency and increases the cost of production. Therefore, there is a need to study the rate of paraffin wax deposition and cloud point temperature in order to guide the oil field control the paraffin wax deposition. In this paper, we use the Flory theory of polymer solution to correct the liquid activity coefficients, and regular solution theory to correct for the non ideality of the solid mixture, and we consider the impact of isoparaffin. Finally, thermodynamic model is established. The actual example calculation shows that the forecast results of this model are more accurate.
 
</p></abstract><kwd-group><kwd>Paraffin Wax</kwd><kwd> Deposition</kwd><kwd> Thermodynamic Model</kwd><kwd> Activity Coefficients</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Preface</title><p>Paraffin wax deposition issues have a significant impact on oilfield production. The deposition of paraffin wax in the tubing and pipe can reduce the production efficiency. In the meantime, removing the deposition of paraffin wax will increase a lot of expenses. The process of paraffin wax deposition is a very complicated problem. On one hand, because the composition of the oil and gas system is very complex, there needs further research on the effect of various components on paraffin deposition. On the other hand, many theoretical issues are involved in the process of paraffin deposition as solubility, crystal, fluid dynamics, mass transfer dynamics, heat transfer, etc. The mechanism of paraffin deposition is not yet fully understood; there are many kinds of explanation theory as solubility theory, crystallization theory, diffusion theory, and phase equilibrium theory [<xref ref-type="bibr" rid="scirp.62160-ref1">1</xref>] . At present, phase equilibrium theory was widely accepted. According to this theory, reason of paraffin wax deposition is the thermodynamic condition of oil and gas system, such as the change of oil and gas system composition, component, temperature, and pressure [<xref ref-type="bibr" rid="scirp.62160-ref2">2</xref>] .</p><p>Model established in this paper can be used to determine cloud point of paraffin wax and the deposited amount of paraffin with the lowering of the temperature, which has important reference value for oilfield to solve paraffin deposition problems.</p></sec><sec id="s2"><title>2. The Thermodynamic Model</title><p>When the phase equilibrium of gas, liquid, solid phase is reached [<xref ref-type="bibr" rid="scirp.62160-ref3">3</xref>] .</p><disp-formula id="scirp.62160-formula786"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x6.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.62160-formula787"><label>. (2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x7.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula788"><label>. (3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x8.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula789"><label>. (4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x9.png"  xlink:type="simple"/></disp-formula><p>If considering the heat tolerance of the liquid and solid, fugacity of solid standard state and liquid standard state has the following relationship.</p><disp-formula id="scirp.62160-formula790"><label>. (5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x10.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.62160-formula791"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x11.png"  xlink:type="simple"/></disp-formula><p>where b<sub>1</sub>, b<sub>2</sub> get from the experimental fitting. In the absence of experimental data b<sub>1</sub>, b<sub>2</sub> can get</p><disp-formula id="scirp.62160-formula792"><label>. (7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x12.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula793"><label>. (8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x13.png"  xlink:type="simple"/></disp-formula><p>Take (6) into (5), (5) becomes</p><disp-formula id="scirp.62160-formula794"><label>. (9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x14.png"  xlink:type="simple"/></disp-formula><p>The non ideal of the solid mixture is corrected by using the regular solution theory, then derived the calculation formula for the activity coefficient of solid phase.</p><disp-formula id="scirp.62160-formula795"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x15.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.62160-formula796"><label>. (11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x16.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula797"><label>. (12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x17.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula798"><label>. (13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x18.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula799"><label>. (14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x19.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula800"><label>. (15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x20.png"  xlink:type="simple"/></disp-formula><p>Use the Flory theory of polymer solution to correct the liquid activity coefficients, the calculation formula as follows [<xref ref-type="bibr" rid="scirp.62160-ref4">4</xref>] .</p><disp-formula id="scirp.62160-formula801"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x21.png"  xlink:type="simple"/></disp-formula><p>where</p><disp-formula id="scirp.62160-formula802"><label>. (17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x22.png"  xlink:type="simple"/></disp-formula><p>Molecular interaction is divided into three categories: (1) PN―PN (2) PN―A (3) A―A. Among which P stands for alkane, N stands for cycloalkane, A stands for aromatic. Their molecular interaction coefficients are as follows:</p><disp-formula id="scirp.62160-formula803"><label>. (18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x23.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula804"><label>. (19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x24.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula805"><label>. (20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x25.png"  xlink:type="simple"/></disp-formula><p>Paraffin wax is mainly composed by the n-alkane, but contains a small amount of isoparaffin, cycloalkane and aromatic. For these non-normal paraffin, their melting point and other properties are very different from the n-alkane, the melting point of normal paraffin wax decreases obviously when a branch is added. For example: the melting point of C<sub>43</sub>H<sub>88</sub> is 188˚F, it is generally believed that the melting point of C<sub>52</sub>H<sub>104</sub> is relatively high. But because there is a branched chain in C<sub>52</sub>H<sub>104</sub>, its melting point will be lower than C<sub>43</sub>H<sub>88</sub> nearly 100˚F, it’s 91˚F. Therefore, it is necessary to consider the effect of non-normal paraffin [<xref ref-type="bibr" rid="scirp.62160-ref5">5</xref>] .</p><p>For normal structure paraffin wax:</p><disp-formula id="scirp.62160-formula806"><label>. (21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x26.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.62160-formula807"><label>. (22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x27.png"  xlink:type="simple"/></disp-formula><p>For non-normal structure paraffin wax:</p><disp-formula id="scirp.62160-formula808"><label>. (23)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x28.png"  xlink:type="simple"/></disp-formula><p>For isoparaffin, cycloalkane:</p><disp-formula id="scirp.62160-formula809"><label>. (24)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x29.png"  xlink:type="simple"/></disp-formula><p>For aromatic:</p><disp-formula id="scirp.62160-formula810"><label>. (25)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-8102506x30.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. The Example Calculation and Results</title><p>In order to validate the model, calculated by using the data of [<xref ref-type="bibr" rid="scirp.62160-ref6">6</xref>] (see <xref ref-type="table" rid="table1">Table 1</xref>) (SRK state equation is used in calculation). Calculation results are shown in <xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>From <xref ref-type="fig" rid="fig1">Figure 1</xref> we can see that the cloud point temperature is different under different pressures. This difference is because of the combined effect of pressure and composition. Temperature decreases with increasing pressure. Under 1.5 MPa, the measured cloud point temperature is 318.7 K, calculating cloud point temperature is 318 k, measured and calculating results are quite close. It can be seen from the <xref ref-type="fig" rid="fig2">Figure 2</xref> that the theoretical results are in good agreement with the actual results.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> The components of a kind of crude oil in the Beihai Oilfield</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Components</th><th align="center" valign="middle" >Mole fraction</th><th align="center" valign="middle" >Density</th><th align="center" valign="middle" >Molecular weight</th></tr></thead><tr><td align="center" valign="middle" >N<sub>2</sub></td><td align="center" valign="middle" >0.000</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >CO<sub>2</sub></td><td align="center" valign="middle" >0.028</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >C1</td><td align="center" valign="middle" >0.128</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >C2</td><td align="center" valign="middle" >0.240</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >C3</td><td align="center" valign="middle" >1.186</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >IC4</td><td align="center" valign="middle" >1.235</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >NC4</td><td align="center" valign="middle" >3.792</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >IC5</td><td align="center" valign="middle" >1.586</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >NC5</td><td align="center" valign="middle" >2.567</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >C6</td><td align="center" valign="middle" >5.235</td><td align="center" valign="middle" >0.6701</td><td align="center" valign="middle" >84.40</td></tr><tr><td align="center" valign="middle" >C7</td><td align="center" valign="middle" >8.152</td><td align="center" valign="middle" >0.7260</td><td align="center" valign="middle" >93.58</td></tr><tr><td align="center" valign="middle" >C8</td><td align="center" valign="middle" >9.418</td><td align="center" valign="middle" >0.7346</td><td align="center" valign="middle" >109.21</td></tr><tr><td align="center" valign="middle" >C9</td><td align="center" valign="middle" >7.024</td><td align="center" valign="middle" >0.7621</td><td align="center" valign="middle" >124.00</td></tr><tr><td align="center" valign="middle" >C10+</td><td align="center" valign="middle" >60.730</td><td align="center" valign="middle" >0.8997</td><td align="center" valign="middle" >345.63</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The diagram of cloud point temperature</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-8102506x31.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> The diagram of wax deposition content</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-8102506x32.png"/></fig></sec><sec id="s4"><title>4. Conclusions</title><p>1) The model proposed in this paper is more accurate, and the prediction results are in good agreement with the practical data.</p><p>2) The model considers the non ideality of the solid phase mixture, and the activity coefficients of solid phase are corrected by using the regular solution theory.</p><p>3) The model considers the influence of the non-normal paraffin; it makes the results more accurate.</p></sec><sec id="s5"><title>Funding</title><p>This work was supported by the National Natural Science Foundation of China (No: 51404072).</p></sec><sec id="s6"><title>Cite this paper</title><p>BaojunLiu,WantingSun,ChengtingLiu,LipingGuo, (2015) The Thermodynamic Model on Paraffin Wax Deposition Prediction. Engineering,07,827-832. doi: 10.4236/eng.2015.712072</p></sec><sec id="s7"><title>Explanation of Nomenclature</title><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x33.png" xlink:type="simple"/></inline-formula>fugacity;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x34.png" xlink:type="simple"/></inline-formula> fugacity of solid standard state;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x35.png" xlink:type="simple"/></inline-formula> fugacity of liquid standard state;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x36.png" xlink:type="simple"/></inline-formula> mole fraction;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x37.png" xlink:type="simple"/></inline-formula> activity coefficient;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x38.png" xlink:type="simple"/></inline-formula> fugacity coefficient;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x39.png" xlink:type="simple"/></inline-formula> molar enthalpy of fusion;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x40.png" xlink:type="simple"/></inline-formula> melting point temperature;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x41.png" xlink:type="simple"/></inline-formula> specific heat under constant pressure;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x42.png" xlink:type="simple"/></inline-formula> coefficient;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x43.png" xlink:type="simple"/></inline-formula> molecular weight;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x44.png" xlink:type="simple"/></inline-formula> molar volume of component;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x45.png" xlink:type="simple"/></inline-formula> solubility parameter;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x46.png" xlink:type="simple"/></inline-formula> carbon number of component i;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x47.png" xlink:type="simple"/></inline-formula> mole number of component i;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x48.png" xlink:type="simple"/></inline-formula> coefficient of molecular interaction;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x49.png" xlink:type="simple"/></inline-formula>pressure;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x50.png" xlink:type="simple"/></inline-formula> gas constant;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x51.png" xlink:type="simple"/></inline-formula> temperature.</p><p>Superscript:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x52.png" xlink:type="simple"/></inline-formula> gas phase;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x53.png" xlink:type="simple"/></inline-formula> solid phase;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x54.png" xlink:type="simple"/></inline-formula> liquid phase;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x55.png" xlink:type="simple"/></inline-formula> standard state;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x56.png" xlink:type="simple"/></inline-formula>liquate.</p><p>Subscript:</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x57.png" xlink:type="simple"/></inline-formula> pure component;</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-8102506x58.png" xlink:type="simple"/></inline-formula> constant pressure.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.62160-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Weingarten, J.S., et al. (1988) Methods for Prediction Wax Precipitation and Deposition. SPE Production Engineering, 3.</mixed-citation></ref><ref id="scirp.62160-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Brown, T.S. (1993) Measurement and Prediction of the Kinetics of Paraffin Deposition. SPE Annual Technical Conference and Exhibition, Houston, 3-6 October 1993. http://dx.doi.org/10.2118/26548-ms</mixed-citation></ref><ref id="scirp.62160-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Won, K.W. (1986) Thermodynamic for Solid Solution-Liquid-Vapor Equilibria: Wax Phase Formation from Heavy Hydrocarbon Mixtures. Fluid Phase Equilibria, 30, 265-279. http://dx.doi.org/10.1016/0378-3812(86)80061-9</mixed-citation></ref><ref id="scirp.62160-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Erickon, D.D., Niesen, V.G., et al. (1993) Thermodynamic Measurement and Prediction of Paraffin Precipitation in Crude Oil. SPE 26604, 933-948.</mixed-citation></ref><ref id="scirp.62160-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Pedersen, K.S., et al. (1995) Prediction of Cloud Point Temperatures and Amount of Wax Precipitation. SPE Production &amp; Facilities, 10, 46-49. http://dx.doi.org/10.2118/27629-PA</mixed-citation></ref><ref id="scirp.62160-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Pan, H.Q., et al. (1997) Pressure and Composition Effect on Wax Precipitation Experimental Data and Model. SPE Production &amp; Facilities, 12. http://dx.doi.org/10.2118/36740-pa</mixed-citation></ref></ref-list></back></article>