<?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">JPEE</journal-id><journal-title-group><journal-title>Journal of Power and Energy Engineering</journal-title></journal-title-group><issn pub-type="epub">2327-588X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jpee.2018.69003</article-id><article-id pub-id-type="publisher-id">JPEE-87165</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>
 
 
  Optimization of Temperature Drop Formula for Annular Gathering Pipeline in High Water Cut Oil Region of Central Hebei Province
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Haonan</surname><given-names>Cai</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>Lining</surname><given-names>Zhao</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>Xuemei</surname><given-names>Chen</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>Xintong</surname><given-names>Chen</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>Tong</surname><given-names>Guan</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>Jinxiu</surname><given-names>Wang</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>Bolin</surname><given-names>Lv</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>Shuaibin</surname><given-names>Li</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Engineering Technology Institute of PetroChina Huabei Oil Field, Renqiu, China</addr-line></aff><pub-date pub-type="epub"><day>07</day><month>09</month><year>2018</year></pub-date><volume>06</volume><issue>09</issue><fpage>21</fpage><lpage>26</lpage><history><date date-type="received"><day>10,</day>	<month>July</month>	<year>2018</year></date><date date-type="rev-recd"><day>4,</day>	<month>September</month>	<year>2018</year>	</date><date date-type="accepted"><day>7,</day>	<month>September</month>	<year>2018</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 order to respond to the national policy of energy saving and consumption reduction, the temperature drop formula of annular pipeline should be established to optimize the water mixing temperature and water mixing flow. The traditional temperature drop formula applies the high water cut actual situation to 13.7%. On the basis of the traditional Schulhof temperature drop formula, the influence of hydraulic friction on the oil flow temperature is calculated by the calculus idea. The heat transfer coefficient of the pipeline under the complex condition is obtained by the inverse algorithm, and the water mixing flow and water temperature are simulated with MATLAB. The function surface diagram of the degree and the back station temperature. The actual error can be reduced to 4.8%, and the return station temperature can be optimized by 7 degrees. 
  
 
</p></abstract><kwd-group><kwd>Circular Water Collecting Oil</kwd><kwd> Hydraulic Friction</kwd><kwd> Inverse Algorithm</kwd><kwd> Saving Energy Cost</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>East oilfield going most of the oil gradually stepped into the stage of high water cut exploitation (moisture content of 90% or more), most of the oil collecting process adopted cricoid watered oil gathering process, and the energy consumption of oil collecting part take up little 65% - 80% of total energy consumption, for countries to provide energy at the same time, and consumes a lot of energy, as the national energy conservation and consumption reduction policy enforcement, annular watered oil gathering process saving energy and reducing consumption a difficult task, and in the annular oil gathering process in the past, water temperature and water flow are two key parameters adjusted by production supervisors with artificial experience, now need to through the establishment of theory of pipeline medium temperature drop formula, To calculate the optimal water blending temperature and water blending flow value, so as to achieve the purpose of energy saving and consumption reduction.</p><p>However, in the traditional theoretical equation of temperature drop of pipeline media, the schuhof temperature drop formula is usually adopted [<xref ref-type="bibr" rid="scirp.87165-ref1">1</xref>]. Second, in the pipeline buried underground, it is necessary to consider the heat dissipation effect between fluid and pipe wall, wall with thermal insulation layer between the heat dissipation effect, thermal insulation layer and soil between the cooling effect and so on the many kinds of complicated conditions, the total friction coefficient of pipe in the traditional engineering manual already does not apply to the complicated working condition.</p></sec><sec id="s2"><title>2. Mathematical Model</title><p>In the process of annular oil collection, the change of heat is mainly reflected in the two processes of pipeline heat dissipation and single well fluid mixing during linear transmission [<xref ref-type="bibr" rid="scirp.87165-ref2">2</xref>]. Straight pipes of heat mainly adopts the revised ShuHuoFu temperature drop formula, and in the process of single well output liquid mixture, were it not for the heat loss by the collision, the temperature of the oil pipeline is the starting point should be for a period of piping within the medium mixed with oil well produced fluid temperature, after finishing the conclusion:</p><p>T j = C W G W j T W + ∑ j = 1 n G y j C y j T y C j ( G W j + ∑ i = 1 n G y j ) (1)</p><p>For the heat capacity of the mixed liquid of J/(kg∙˚C), for mixed fluid flow kg/s, as the temperature of the mixed liquid ˚C, for the well produced fluid flow kg/s, as well the heat capacity of the produced fluid J/(kg∙˚C), for the mouth of the well produced fluid temperature ˚C.</p><sec id="s2_1"><title>2.1. Modify the Schuhof Temperature Drop Formula</title><p>In the process of transmission line, the traditional model using ShuHuoFu temperature drop formula, namely from the starting point of the fluid temperature L, starting point for the temperature, the temperature of the surrounding medium for D m pipe diameter, pipe length L, m, pipe quality of oil-water mixture flow G kg/s, pipeline total heat transfer coefficient KW/(m<sup>2</sup>・˚C).</p><p>T L = T 0 + ( T Q − T 0 ) e K π D L C G (2)</p><p>But in practice, due to the process of crude oil in pipe flow, with the inner wall of pipeline caused by friction, along the pipeline pressure is reduced, to reduce the pressure will be converted into friction heat and fluid inside the pipeline heating [<xref ref-type="bibr" rid="scirp.87165-ref3">3</xref>]. The heat loss caused by pipeline hydraulic friction is not considered in the traditional schuhof temperature drop formula.</p><p>Revised ShuHuoFu temperature drop formula, is the premise of considering the friction heat, assuming a pipeline diameter D, G throughput and the temperature of the surrounding medium for, take the oil temperature in dL infinitesimal section is T, take for hydraulic grade I, crude oil temperature drop for dT produced on the dL. Under stable conditions, the heat balance equation of dL segment is:</p><p>K D π ( T − T 0 ) d L = − G C d T + G g i d L (3)</p><p>The left part of the equation above is the heat dissipation from the pipeline to the surrounding medium, and the right part is the sum of heat release from oil flow temperature drop and heat of friction resistance [<xref ref-type="bibr" rid="scirp.87165-ref4">4</xref>]. As dL goes up, dT goes down, so minus sign is introduced [<xref ref-type="bibr" rid="scirp.87165-ref5">5</xref>]. After the above equation is deformed, the sukhoff temperature drop formula considering friction-heat generation can be obtained, that is, the modified sukhoff temperature drop formula is:</p><p>T L = ( T 0 + b ) + [ T R − ( T 0 + b ) ] e − a L (4)</p><p>a = K D π G C , b = g i C a</p></sec><sec id="s2_2"><title>2.2. Fixed the Total Heat Transfer Coefficient of the Pipeline</title><p>Most of the gathering and transportation pipelines are buried pipelines, whose heat transfer is divided into three parts: heat dissipation between the fluid and the inner wall of the pipeline; Heat conduction between the outer wall of the pipe and the insulating layer or the insulating layer; Heat dissipation between the insulation layer and the surrounding environment. The total heat transfer coefficient of crude oil pipeline can be expressed as:</p><p>K = 1 D α 1 D 1 + ∑ ​ D 2 γ I ln D ( i + 1 ) D i + D α 2 D w (5)</p><p>where K is the total heat transfer coefficient of the pipeline W/(m<sup>2</sup>・˚C); α 1 , the heat release coefficient W/(m<sup>2</sup>・˚C), α 2 , between oil flow and pipe inner wall, and is the heat release coefficient W/(m<sup>2</sup>・˚C) γ I , between pipe wall and soil. Is the heat conduction coefficient W/(m<sup>2</sup>・˚C) corresponding to each layer, and D is the calculated pipe diameter. If there is no insulation layer, the pipe outer diameter can be directly adopted. If there is an insulating layer, the inner diameter of the insulating layer m is adopted. D is the inner and outer diameter m of steel tube and asphalt insulation layer; Is the pipe outer diameter m.</p><p>However, in the actual production process, due to complex field conditions, the heat transfer coefficient K calculated by the above formula has a large error compared with the actual heat transfer coefficient [<xref ref-type="bibr" rid="scirp.87165-ref6">6</xref>]. In order to reduce this error, the inverse algorithm is usually used to calculate the heat transfer coefficient K of the pipeline.</p><p>T R as the section temperature ˚C, T Z as section temperature ˚C, L R for the section length m.</p><p>K = G C D L R π ln T R − T 0 T Z − T 0 (6)</p></sec></sec><sec id="s3"><title>3. Model Validation and Analysis</title><p>In order to verify the accuracy of the models and applied, the use of Huabei oilfield third production plant is high moisture content of the annular watered oil collecting data (month 12, water temperature, water flow rate, back to the station temperature, pipe diameter size, etc.), and the simulation of temperature drop data and compare [<xref ref-type="bibr" rid="scirp.87165-ref7">7</xref>]. First, use Pipesim to draw the well location diagram, as shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>According to the ring of the outbound parameters such as temperature, the temperature of back to the station, the inverse calculation formula to calculate the total heat transfer coefficient K, each output of the well fluid flow, temperature, pipe parameters, buried temperature constant coefficient substitution model, the water flow rate set for x G, water blending temperature T is set to y, back to the station temperature as the dependent variable is set to z. After the mesh treatment of the model with matlab, the relationship between the return station temperature and water incorporation rate and water incorporation temperature is obtained, as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p>The actual water temperature, water flow into the traditional model (using the traditional ShuHuoFu total heat transfer coefficient of the model and the direct use of engineering manual) and the revised model respectively to calculate the temperature, back to the station to get the following <xref ref-type="table" rid="table1">Table 1</xref>; The actual temperature of the return station was marked on the graph of the function surface of water incorporation rate, water incorporation temperature and temperature of the return station, as shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>.</p><p>According to the data obtained, the error between the calculated results of the traditional model and the actual return station temperature is approximately 13.7%. The error between the calculated results of the modified model and the actual return temperature is approximately 4.8%. It can be seen that the calculated</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Water incorporation rate, water incorporation temperature statistics table</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Date</th><th align="center" valign="middle" >Water flow m<sup>3</sup>/d</th><th align="center" valign="middle" >The water temperature ˚C</th><th align="center" valign="middle" >Initial return temperature ˚C</th><th align="center" valign="middle" >Traditional model calculation ˚C</th><th align="center" valign="middle" >New model calculation ˚C</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >85.44</td><td align="center" valign="middle" >66.64</td><td align="center" valign="middle" >45.12</td><td align="center" valign="middle" >51.23</td><td align="center" valign="middle" >47.23</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >85.44</td><td align="center" valign="middle" >66.49</td><td align="center" valign="middle" >45.52</td><td align="center" valign="middle" >53.26</td><td align="center" valign="middle" >46.89</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >85.92</td><td align="center" valign="middle" >66.35</td><td align="center" valign="middle" >45.54</td><td align="center" valign="middle" >50.12</td><td align="center" valign="middle" >47.68</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >85.44</td><td align="center" valign="middle" >66.03</td><td align="center" valign="middle" >45.37</td><td align="center" valign="middle" >52.21</td><td align="center" valign="middle" >47.23</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >66.06</td><td align="center" valign="middle" >45.46</td><td align="center" valign="middle" >52.33</td><td align="center" valign="middle" >48.13</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >65.94</td><td align="center" valign="middle" >45.75</td><td align="center" valign="middle" >53.11</td><td align="center" valign="middle" >48.47</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >66.17</td><td align="center" valign="middle" >45.86</td><td align="center" valign="middle" >52.32</td><td align="center" valign="middle" >47.99</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >86.64</td><td align="center" valign="middle" >66.29</td><td align="center" valign="middle" >46.04</td><td align="center" valign="middle" >52.11</td><td align="center" valign="middle" >48.24</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >66.25</td><td align="center" valign="middle" >45.96</td><td align="center" valign="middle" >53.08</td><td align="center" valign="middle" >48.38</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >66.19</td><td align="center" valign="middle" >45.89</td><td align="center" valign="middle" >51.65</td><td align="center" valign="middle" >48.24</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >86.4</td><td align="center" valign="middle" >66.15</td><td align="center" valign="middle" >45.83</td><td align="center" valign="middle" >51.22</td><td align="center" valign="middle" >48.13</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >86.64</td><td align="center" valign="middle" >66.41</td><td align="center" valign="middle" >45.89</td><td align="center" valign="middle" >51.35</td><td align="center" valign="middle" >48.29</td></tr><tr><td align="center" valign="middle" >average</td><td align="center" valign="middle" >86.16</td><td align="center" valign="middle" >66.2475</td><td align="center" valign="middle" >45.68583</td><td align="center" valign="middle" >51.99917</td><td align="center" valign="middle" >47.90833</td></tr></tbody></table></table-wrap><p>error of the temperature drop formula obtained after the modification is significantly reduced.</p></sec><sec id="s4"><title>4. Conclusions</title><p>1) After correction ShuHuoFu temperature drop formula and the total heat transfer coefficient, in view of the north China, Huabei Oilfield chu a high water cut oil region, the revised model calculated back to the station temperature compared with traditional model calculated results error is reduced from 13.7% to 4.8%, the temperature drop of the optimized model for going more in line with high water cut oil ring watered oil temperature drop model, it has application value.</p><p>2) Get a more accurate temperature drop model, generally known crude oil freezing point of 35˚C, considering the safety of 3˚C temperature coefficient, the optimal back stand temperature for 38˚C, and now back to the station temperature 45˚C. In the annular watered oil collecting the existing process, through the model can further adjust the water flow and water temperature, finally can also be made back to the station 7˚C to 38˚C temperature drop, so as to further reduce energy consumption.</p></sec><sec id="s5"><title>Cite this paper</title><p>Cai, H.N., Zhao, L.N., Chen, X.M., Chen, X.T., Guan, T., Wang, J.X., Lv, B.L. and Li, S.B. (2018) Optimization of Temperature Drop Formula for Annular Gathering Pipeline in High Water Cut Oil Region of Central Hebei Province. Journal of Power and Energy Engineering, 6, 21-26. https://doi.org/10.4236/jpee.2018.69003</p></sec></body><back><ref-list><title>References</title><ref id="scirp.87165-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Su, H. (2017) Optimization of Energy Consumption for Crude Oil Pipeline. Xi’an Petroleum University.</mixed-citation></ref><ref id="scirp.87165-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Wu, C.C. and Yan, D.F. (1989) Two Level Hierarchical Optimization Model for Steady Operation of Hot Oil Pipeline. Acta Sinica, 10, 109-117.</mixed-citation></ref><ref id="scirp.87165-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Wang, B.Q., Lin, Y.H., Dai, Y.B., et al. (2012) Current Situation and Prospect of China’s Crude Oil Pipeline. Petroleum Planning and Design, 23, 8-11.</mixed-citation></ref><ref id="scirp.87165-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Liu, Y. and Guan, X.J. (1993) Optimal Design of Oil and Gas Gathering and Transportation System. Natural Gas Industry, No. 2, 36-38.</mixed-citation></ref><ref id="scirp.87165-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Liu, Y. and Guan, X.J. (1993) Topology Optimization Design of Annular Gathering Pipeline Network. Acta Petrota Sinica, No. 3, 28-30.</mixed-citation></ref><ref id="scirp.87165-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Shi, G.Y. and Dong, J.L. (1999) Optimization Method. Higher Education Press, 53-57.</mixed-citation></ref><ref id="scirp.87165-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Jiang, J.S., He, C.X. and Pan, S.H. (2007) Optimization Calculation Method South China University of Technology Press, 63-217.</mixed-citation></ref></ref-list></back></article>