CN116738113A - 预测地热井筒中结垢位置的数学模型的构建方法 - Google Patents
预测地热井筒中结垢位置的数学模型的构建方法 Download PDFInfo
- Publication number
- CN116738113A CN116738113A CN202310405346.9A CN202310405346A CN116738113A CN 116738113 A CN116738113 A CN 116738113A CN 202310405346 A CN202310405346 A CN 202310405346A CN 116738113 A CN116738113 A CN 116738113A
- Authority
- CN
- China
- Prior art keywords
- geothermal fluid
- geothermal
- fluid
- wellhead
- flow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000013178 mathematical model Methods 0.000 title claims abstract description 16
- 238000010276 construction Methods 0.000 title abstract description 8
- 239000012530 fluid Substances 0.000 claims abstract description 230
- 239000012071 phase Substances 0.000 claims abstract description 28
- 239000007791 liquid phase Substances 0.000 claims abstract description 17
- 230000006835 compression Effects 0.000 claims abstract description 13
- 238000007906 compression Methods 0.000 claims abstract description 13
- 230000008020 evaporation Effects 0.000 claims abstract description 12
- 238000001704 evaporation Methods 0.000 claims abstract description 12
- 239000007788 liquid Substances 0.000 claims description 33
- 238000000034 method Methods 0.000 claims description 32
- 230000008569 process Effects 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000003595 mist Substances 0.000 claims description 8
- 230000007704 transition Effects 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000009529 body temperature measurement Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000005086 pumping Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 53
- VTYYLEPIZMXCLO-UHFFFAOYSA-L Calcium carbonate Chemical compound [Ca+2].[O-]C([O-])=O VTYYLEPIZMXCLO-UHFFFAOYSA-L 0.000 description 10
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 10
- 230000005514 two-phase flow Effects 0.000 description 7
- 239000000126 substance Substances 0.000 description 6
- 229910021532 Calcite Inorganic materials 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 5
- 229910000019 calcium carbonate Inorganic materials 0.000 description 5
- 229910002092 carbon dioxide Inorganic materials 0.000 description 5
- 239000001569 carbon dioxide Substances 0.000 description 5
- 239000003208 petroleum Substances 0.000 description 5
- 238000004088 simulation Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000007797 corrosion Effects 0.000 description 4
- 238000005260 corrosion Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000000630 rising effect Effects 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000003889 chemical engineering Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 240000007124 Brassica oleracea Species 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 241000237858 Gastropoda Species 0.000 description 1
- 238000009360 aquaculture Methods 0.000 description 1
- 244000144974 aquaculture Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012412 chemical coupling Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000004090 dissolution Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 239000010438 granite Substances 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- -1 second section Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000009491 slugging Methods 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Mathematical Analysis (AREA)
- Quality & Reliability (AREA)
- Pure & Applied Mathematics (AREA)
- Game Theory and Decision Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
本发明公开了预测地热井筒中结垢位置的数学模型的构建方法,步骤为:(1)测量地热井中井口处的地热流体的流量、温度、压力以及流体特性;(2)计算地热流体中的气体压缩系数;(3)计算地热流体中修正后的液相密度和气相密度;(4)确定在计算出地热流体所处的深度内所处的流型;(5)确定一次增压井深位置处的雷诺数和相对粗糙度;(6)计算地热流体与井筒之间的摩擦力;(7)计算地热流体在井筒中的闪蒸深度;(8)计算地热流体预测结垢位置的相对误差。本发明所述的一种预测地热流体在井筒中的结垢位置的数学模型的构建方法,模型简捷,不用复杂的数学推导,能够更快的预测地热流体在井筒中的闪蒸位置,适合工业化应用。
Description
技术领域
本发明涉及数学模型涉及地热流体在井筒中的结垢阻垢,具体地说,涉及一种预测地热流体在井筒中的结垢位置的数学模型的构建方法。
背景技术
地热能的储藏十分丰富,是重要的新能源之一(蔡义汉.地热直接利用[M].天津:天津大学出版社,2014.)。地热能可用于发电、供暖、水产养殖以及化学工业等领域(IngvarB,Fridleifsson.Geothermal energy for the benefit of the people(地热能造福人民).Renewable&Sustainable Energy Reviews(可再生和可持续能源评论),2001,5(3):299-312.)。我国的地热***主要分为以羊八井为代表的高温花岗岩储层、咸阳为代表的中低温砂岩储层和雄县为代表的中低温碳酸盐储层(陈俊光.羊八井地热双工质发电管理***的设计与研究[D].成都:电子科技大学,2021.)。但是在开发利用地热能的同时也遇到了各种各样的问题,其中地热流体在井筒中的结垢就是一个很严重的问题(Boch R,Leis A,Haslinger E,et al.Scale-fragment formation impairing geothermal energyproduction:interacting H2S corrosion and CaCO3 crystal growth(结垢形成损害地热能生产:H2S腐蚀和CaCO3晶体生长相互作用).Geothermal Energy(地热能),2017,5(1):1-19.)。地热流体在井筒中的结垢位置在一定程度上是基于质量守恒、动量守恒和能量守恒原理来建立模型(Arnórsson S.Deposition of calcium carbonate minerals fromgeothermal waters-theoretical considerations(地热水中碳酸钙矿物沉积的理论思考).Geothermics(地热学),1989,18(1/2):33-39.),并耦合了化学热力学的相关方程进行求解(D,Kohl T.Coupled thermal–hydraulic–chemical modelling of enhancedgeothermal systems(强化地热***的热力-水力-化学耦合建模).Geophysical JournalInternational(国际地球物理学杂志),2005,161(2):533-548.)。进而开发了许多软件和相关程序来模拟地热流体在井筒中的结垢过程(Bankoff S G.A variable densitysingle-fluid model for two-phase flow with particular reference to steam-water flow(两相流变密度单流体模型,特别参考蒸汽-水流动).Journal of HeatTransfer(传热杂志),1960,82(4):265-272.)。通过井口和井底的温度和压力可计算出地热流体在井筒中的温度梯度和压力梯度,当温度梯度和压力梯度发生明显突变时,地热流体就会在明显变化的位置处发生闪蒸,之后产生结垢。因此由地热流体在井筒中的温度梯度和压力梯度的突变位置就可以得到闪蒸点的位置,进而得到结垢位置。
井筒物理模型通常描述为在较深区域有液体流动的垂直圆形管道,随着地热流体沿井筒上升,地热流体的动量和热量都会损失,因此准确估计热损失对于模拟地热井中的流动非常重要。地热流体在地热井筒中由于压力和温度的降低而在上部区域出现闪蒸。从闪蒸位置处开始,随着混合速度和蒸汽质量的增加,产生两相流(液体-蒸汽)。而气相中水蒸气的不断增加,很容易稀释气相中的二氧化碳气体,造成液相地热流体中溶解的二氧化碳快速释放到气相中,导致地热流体中的碱度增加,而这时也就容易产生结垢现象。根据这个规律,可预测结垢位置,其实也可以通过预测闪蒸点位置来间接得到。在向上流动的过程中,根据蒸汽和液体的相对量或相速度,可能会出现几种流型模式(例如,气泡流、段塞流、搅动流或半环状流、分散环状流和环状流)。由于蒸汽(气相)和液相之间的滑动,向上垂直和倾斜管道中的相分布相当复杂。当混合物上升到井口时,蒸汽被分离并用于发电,而残留的液体(或盐水)被重新注入地下。分析地热井内的两相流需要计算生产过程参数(例如压力、温度、焓、热通量和速度剖面)和流体热物理特性。存在的问题是使用不合适的关系式确定一些两相流参数(例如,空隙率、相含率和摩擦系数)。地热流体在井筒中的具体流型如下:
(1)气泡流;气体以小气泡的形式存在,随机分布,其直径也随机变化。气泡根据它们各自的直径以不同的速度运动,液体以很均匀的速度沿管道向上流动。除密度外,气相对压力梯度几乎没有影响。管道内几乎完全充满液体,游离的气相很少。
(2)弹状流;液相仍是连续的,小气泡开始聚结并形成大小和形状几乎相同的稳定气泡,其直径几乎与管道直径相同,被液体块隔开。气泡速度大于液体速度。气泡周围有一层液体薄膜,液体速度开始不稳定,气相和液相都对压力梯度有显著影响。
(3)过渡流;在该区域发生了从连续液相到连续气相的转变。液相被夹带在气相当中,气相开始占据了优势。
(4)雾状流;气相是连续的,液相被气相夹带和携带,气相为主要的控制因素。
Hughmark和Pressburg在1961年最先将相含率应用于地热井储模拟器中(Hughmark GA,Pressburg B S.Hold up and pressure drop with gas-liquid flow ina vertical pipe(垂直管道中气液流的持液率和压降).AIChE(美国化学工程师学会).1961,7(4):677-682.),之后Duns和Ros设计了一种垂直井筒中的流体流动模式图(DunsH,Ros N C J.Vertical flow of gas and liquid mixtures in wells(井内气液混合物的垂直流动).In:Proceedings 6th World Petroleum Congress,Section II,Paper 22–PD 6,Frankfurt am Main,Germany,pp(附:第六届世界石油大会会议记录,第二节,论文22–PD 6,德国法兰克福).1963:451–465.),Pátzay等人使用Davies和Pitzer活度系数计算方法开发了用于CaCO3-H2O-CO2***中方解石结垢形成的平衡模拟算法和计算机程序(Pátzay G,Stáhl G,Kármán FH,et al.Modeling of scale formation and corrosion fromgeothermal water(地热水结垢和腐蚀建模).Electrochimica Acta(电子化学学报),1998,43(1-2):137-147.),用于确定地热井口中不同气体(例如CO2、CH4和N2)浓度下的地热流体在井筒中的闪蒸点位置和气体分压的分布情况。当模型建立之后,可以通过HOLA(Bjornsson G.A multi-feedzone geothermal wellbore simulator(多补给带地热井筒模拟器).Technical Report LBL-23546,Berkeley:Lawrence Berkeley Laboratory(技术报告LBL-23546,伯克利:劳伦斯伯克利实验室),1987:1-102.)和WELLSIM(Gunn C,Freeston D.An integrated steady-state wellbore simulation and analysispackage(集成稳态井筒模拟与分析软件包).Proceedings of the 13th New ZealandGeothermal Workshop.Auckland,NZ:New Zealand Geothermal Workshop(第十三届新西兰地热研讨会论文集。新西兰奥克兰:新西兰地热研讨会),1991:161-166.)等软件程序来计算地热井筒中闪蒸点的位置。这两种软件的区别在于HOLA不考虑地热流体中盐分和不凝气体含量对闪蒸点位置的影响,而WELLSIM则需要同时考虑这两个参数的影响。Garg等人则利用新建立的相含率经验关联式来模拟分析地热井筒中温度和压力随着井筒深度的变化规律来判断闪蒸点的位置(Garg S K,Pritchett J W,Alexander J H.A new liquidhold-up correlation for geothermal wells(一种新的地热井持液率关联式).Geothermics(地热学),2004,33(6):795-817.),进而优化了参数以便最好地代表现场数据。Akin等人利用PHREEQC软件计算出土耳其Kizildere地热田的碳酸钙在闪蒸点之上80米处左右产生结垢(Akin T,Güney A,Kargi H.Modeling of calcite scaling andestimation of gas breakout depth in a geothermal well by using PHREEQC(利用PHREEQC建立方解石结垢模型并估算地热井中的气体漏失深度).Fortieth Workshop onGeothermal Reservoir Engineering(第四十届地热储层工程研讨会).StanfordUniversity,Stanford,California,2015:1-8.)。
美国、菲律宾、意大利和土耳其等国家开发利用地热能比较早,因此积累了丰富的井筒流动和结垢模拟相关经验(Fujii Y.CaCO3 scale problems in the Nigorikawageothermal area,Hokkaido(北海道Nigorikawa地热区CaCO3结垢问题).Japangeothermal association(日本地热协会),1988,25(4):54-65.)。所建立的地热流体结垢位置预测模型通常包括地热流体中碳酸钙的溶解-沉淀规律、二氧化碳等不凝性气体的溶解规律、井筒中单相流和两相流的压降确定、结垢机理分析和闪蒸位置预测等(Rangel G,Pereira V,Ponte C,et al.Reaming calcite deposits of well PV8 whiledischarging:a successful operation at Riberia Grande geothermal field,SáoMiguel Island,Azores(排水时挖掘PV8井的方解石矿床:亚速尔群岛萨奥米格尔岛Riberia Grande地热田的成功作业).European Geothermal Congress 2019Den Haag,TheNetherlands(2019年欧洲地热大会荷兰Den Haag),11-14June 2019.)。
我国在地热流体结垢模拟方面开发利用比较晚,而且针对该问题的积累较少(卜宪标,郭志鹏,王令宝.地热流体在井筒中的流动及碳酸钙结垢过程模拟.新能源进展,2021,9(5):434-442.)。目前对于所建的地热流体结垢位置预测模型有两种方法来求解,一是根据已知的热储或地热井筒的井底参数(压力,温度,总质量流量,二氧化碳含量等),由热储或井底逐步向地热井的井口方向求解,确定各个截面的参数,称为Down-Top算法;二是已知井口参数(压力,温度,干度,总流量,二氧化碳含量),由井口逐步向井底和热储方向求解,称为Top-Down算法。
发明内容
本发明的目的在于克服已有技术的缺点,提供一种能够更快的预测地热流体在井筒中的闪蒸位置,可以确定地热流体的结垢位置的预测地热井筒中结垢位置的数学模型的构建方法。
本发明的预测地热井筒中结垢位置的数学模型的构建方法,包括以下步骤:
步骤一、测量地热井中井口处的地热流体的流量、温度、压力以及流体特性;所述的流体特性包括井口处的地热流体的气体溶解度、标准临界压力、临界温度以及黏度,其中:
气体溶解度:Rs=20.481;
标准临界压力:
临界温度:
黏度:μ=1.0×10-3Pa·s;
步骤二、通过计算得到地热流体中的气体压缩系数,具体过程如下:
步骤201,根据井口处的地热流体的压力值设定井深范围内的压力梯度的初始值,然后计算得到一次增压压力值p1,与一次增压压力值对应的距离井口一定距离的位置记为一次增压井深位置;
式中p0为井口处的压力值,Δp为假设的压力梯度的初始值;
步骤202,利用对比温度和对比压力矫正井口处的地热流体的温度和压力的测量值,计算得到井口处地热流体的无量纲温度Tr和无量纲压降pr;
无量纲温度:
无量纲压降:
式中:为井口处的地热流体温度测量值,Tpe为地热流体的临界温度,/>为井口处地热流体压降测量值,ppe为地热流体的临界压力;
步骤203,根据井口处地热流体的无量纲温度Tr和无量纲压降pr,在无量纲气体压缩系数图中查找得到井口处的地热流体中的气体压缩系数Z;
步骤三、修正井口处的地热流体中的液体体积流量和气体体积流量;然后利用根据经验参数修正地热流体的液体质量流量和气体质量流量;最后再计算地热流体中修正后的液相密度和气相密度,具体公式如下:
地热流体在井口处根据经验参数修正后的液体体积流量为:
qL=6.49×10-5q0
其中,qL为地热流体在井口处根据经验参数修正后的液体体积流量,q0为地热流体井口处的实际体积流量;
地热流体在井口处根据经验参数修正后的气体质量流量为:
其中,qg为根据经验参数修正后的井口处地热流体中气体体积流量;R是气体溶解度,z为井口处的地热流体中的气体压缩系数;
qt=qL+qg
式中:qt为根据经验参数修正后的地热流体在井口处的地热流体总的体积流量;
修正后的井口处的地热流体中的液相质量流量wL为:
wL=q0(4.05×10-8γ0+8.85×10-7γgRs)
其中:γ0为地热流体中水对油的流体比重,γg为地热流体中的气含率;
修正后的井口处的地热流体中的气体质量流量wg和地热流体中总的液体质量流量wt为:
wg=8.85×10-7q0γg(R-Rs)
wt=wL+wg
修正后的井口处的地热流体中的液相密度ρL和气相密度ρg分别为:
步骤四、确定在计算出地热流体所处的深度内所处的流型,具体过程为:
步骤401,计算所述的一次增压井深度处的地热流体的测试变量vt、δ和vgD,所述的测试变量用来确定地热流体所处的边界条件:
式中:vt为一次增压井深位置处的地热流体速度,m/s;vgD为一次增压井深位置处的无量纲气体速度;σ为一次增压井深位置处的流体表面张力,N;Ap为地热井筒的井筒截面积;
边界限制条件为:
式中:(L)B为气泡流-段塞流边界的整体无量纲速度数,dh为地热流体流动的当量直径;(L)S为段塞流-过渡流边界的整体无量纲速度数;(L)M为过渡流-雾状流边界的整体无量纲速度数;
步骤402,分别比较δ,vgD,(L)S和(L)M;进行如下判断,若δ<(L)B,则该压力值所对应井深位置处该地热流体所处的流型为气泡流;若δ>(L)B,vgD<(L)S,则该压力值所对应井深位置处该地热流体所处的流型为弹状流;若(L)M>vgD>(L)S,则该压力值所对应井深位置处该地热流体所处的流型为过渡流;若vgD>(L)M,则该压力值所对应井深位置处该地热流体所处的流型为雾状流;若为δ<(L)B情况,则执行下一步;否则执行步骤201,在原有的压力梯度初始值的基础上将压力梯度值调小;然后执行步骤二-四;
步骤五、确定一次增压井深位置处的地热流体的雷诺数、地热流体中的气泡雷诺数和地热井筒的相对粗糙度;
其中,地热流体的雷诺数NRe:
其中μL为地热流体液体动力粘度系数;
通过对vb进行迭代求解地热流体中的气泡雷诺数Nb,公式如下:
地热流体的平均气体上升速度vbf:
地热流体在井筒中的相对粗糙度为:
其中:ξ为地热井筒的绝对粗糙度,d为井筒内径;
步骤六、计算与一次增压井深位置处对应的地热流体与井筒之间的摩擦力f和地热流体的平均密度具体公式为:
当NRe<2400,
当NRe≥2400,所以地热流体的平均密度/>为:
计算地热流体的摩擦损失梯度τf为:
其中:gc为引力常量;
步骤七、计算地热流体在井筒中的闪蒸深度D,即一次增压井深位置处:
步骤八,计算地热流体预测结垢位置的相对误差为:
其中:ψ1为预测值,即通过上述步骤计算出的地热流体在井筒中的闪蒸深度D;
ψ2为实际值,即通过实际测量得到的地热流体在井筒中的结垢位置;
ω为潜水泵深度,即在地热井筒中潜水泵下潜的深度。
本发明的有益效果:
(1)数学模型简捷,不用很复杂的数学推导,主要通过气液两相流在地热井筒中的不同流型来判断地热流体在井筒不同位置处所处的流型,进而计算出地热流体在井筒中的闪蒸点位置,来判断地热流体在井筒中的结垢位置。
(2)建立数学模型来预测地热流体的结垢位置,最终形成一套通过井口参数简要确定结垢位置的计算方法。根据该计算方法,可以建立从地热井井口到地热井井底不同深度处的地热流体的流型,最终确定地热流体发生闪蒸的位置,可以确定地热流体的结垢位置。
(3)该数学模型能够更快的预测地热流体在井筒中的闪蒸位置,而不需要考虑太复杂的影响因素。
(4)该数学模型简单易行,适合工业化应用。
(5)解决了以往软件计算中的黑箱问题,对地热流体在井筒中发生闪蒸,进而产生沉淀的理论模型的构建。
具体实施方式
下面将结合实施例对本发明作进一步的详细说明。此处所描述的具体实施例方式仅用于解释本发明,并不用于限定本发明的保护范围。
本发明的预测地热井筒中结垢位置的数学模型的构建方法,包括以下步骤:
步骤一、测量地热井中井口处的地热流体的流量、温度、压力以及流体特性;所述的流体特性包括井口处的地热流体的气体溶解度、标准临界压力、临界温度以及黏度,其中,
气体溶解度:Rs=20.481;
标准临界压力:
临界温度:
黏度:μ=1.0×10-3Pa·s;
步骤二、通过计算得到地热流体中的气体压缩系数,具体过程如下:
步骤201,根据井口处的地热流体的压力值设定井深范围内的压力梯度的初始值,然后通过公式(1)计算得到一次增压压力值p1,与一次增压压力值对应的距离井口一定距离的位置记为一次增压井深位置;井口处的地热流体的压力越大,通常压力梯度也越大,压力梯度的的初始值可以设置为0.05-2.0MPa,梯度过大或者过小都会影响最终的预测结果。
公式如下:
式中p0为井口处的压力值,Δp为假设的压力梯度的初始值;
步骤202,根据公式(2)和(3),利用对比温度和对比压力矫正井口处的地热流体的温度和压力的测量值,得到井口处地热流体的无量纲温度Tr和无量纲压降pr;
无量纲温度:
无量纲压降:
式中:为井口处的地热流体温度测量值,Tpe为地热流体的临界温度,/>为井口处地热流体压降测量值,ppe为地热流体的临界压力。
步骤203,根据井口处地热流体的无量纲温度Tr和无量纲压降pr,在无量纲气体压缩系数图(参见2009年化学工业出版社出版的《化工热力学》)中查找得到井口处的地热流体中的气体压缩系数Z;
步骤三、修正井口处的地热流体中的液体体积流量和气体体积流量;之后利用根据经验参数修正地热流体的液体质量流量和气体质量流量;最后再计算地热流体中修正后的液相密度和气相密度,具体公式如下:
地热流体在井口处根据经验参数修正后的液体体积流量为:
qL=6.49×10-5q0 (4)
其中,qL为地热流体在井口处根据经验参数修正后的液体体积流量,q0为地热流体井口处的实际体积流量。
地热流体在井口处根据经验参数修正后的气体质量流量为:
其中,qg为根据经验参数修正后的井口处地热流体中气体体积流量;R是气体溶解度,在此处可取为100-120,z为井口处的地热流体中的气体压缩系数。
公式(4)和(5)的修正参数分别为6.49×10-5和3.27×10-7,可以查找文献《Predicting two-phase pressure drops in vertical pipe在竖直井中的两相压力降预测》(期刊Journal of Petroleum Technology,1967年版,出版页数为829-838)得到。
qt=qL+qg (6)
式中:qt为根据经验参数修正后的地热流体在井口处的地热流体总的体积流量。
修正后的井口处的地热流体中的液相质量流量wL(kg/s)为:
wL=q0(4.05×10-8γ0+8.85×10-7γgRs) (7)
其中:γ0为地热流体中水对油的流体比重,γg为地热流体中的气含率。
修正后的井口处的地热流体中的气体质量流量wg和地热流体中总的液体质量流量wt为:
wg=8.85×10-7q0γg(R-Rs) (8)
wt=wL+wg (9)
修正参数8.85×10-7和4.05×10-8可以查找文献《Predicting two-phasepressure drops in vertical pipe》(期刊Journal of Petroleum Technology,1967年版,出版页数为829-838)得到。
修正后的井口处的地热流体中的液相密度ρL和气相密度ρg分别为:
步骤四、确定在计算出地热流体所处的深度内所处的流型,具体过程为:
步骤401,计算所述的一次增压井深度处的地热流体的测试变量,即公式(12),(13),(14)中所涉及到的变量vt、δ和vgD,所述的测试变量用来确定地热流体所处的边界条件:
式中:vt为一次增压井深位置处的地热流体速度(m/s);vgD为一次增压井深位置处的无量纲气体速度;σ为一次增压井深位置处的流体表面张力(N),从2017年高等教育出版社出版的《化工原理》查表可以得到。Ap为地热井筒的井筒截面积。
边界限制条件为:
式中:(L)B为气泡流-段塞流边界的整体无量纲速度数,dh为地热流体流动的当量直径。(L)S为段塞流-过渡流边界的整体无量纲速度数;(L)M为过渡流-雾状流边界的整体无量纲速度数;
步骤402,分别比较δ,vgD,(L)S和(L)M。进行如下判断,若δ<(L)B,则该压力值所对应井深位置处该地热流体所处的流型为气泡流;若δ>(L)B,vgD<(L)S,则该压力值所对应井深位置处该地热流体所处的流型为弹状流;若(L)M>vgD>(L)S,则该压力值所对应井深位置处该地热流体所处的流型为过渡流;若vgD>(L)M,则该压力值所对应井深位置处该地热流体所处的流型为雾状流;若为δ<(L)B情况,则执行下一步;否则执行步骤201,在原有的压力梯度初始值的基础上将压力梯度值调小。然后执行步骤二-四。
(L)B、(L)S和(L)M是Orkiszewski根据Griffith和Wallis以及Duns和Ross针对其余三种流态提出的气泡和段塞流之间的边界报告的参数,(L)S为段塞流-过渡流边界的整体无量纲速度数,(L)M为过渡流-雾状流边界的整体无量纲速度数;参见文献《Predictingtwo-phase pressure drops in vertical pipe》(期刊Journal of PetroleumTechnology,1967年版,出版页数为829-838)。
表2地热流体中流型的边界条件
表2内容来源于文献《Predicting two-phase pressure drops in verticalpipe》(期刊Journal of Petroleum Technology,1967年版,出版页数为829-838)。
表2中的x为地热流体的蒸汽质量,也叫蒸汽品质,或者蒸汽干度(现场取样测量得到)。
步骤五、确定一次增压井深位置处的地热流体的雷诺数、地热流体中的气泡雷诺数和地热井筒的相对粗糙度;
其中,地热流体的雷诺数NRe
其中μL为地热流体液体动力粘度系数,可以通过查表得到,表参见2017年高等教育出版社出版的《化工原理》;
通过对vb进行迭代求解地热流体中的气泡雷诺数Nb,公式如(19)、(20)、(21):
地热流体的平均气体上升速度vbf:
气泡雷诺数Nb求解过程中由于地热流体中的气泡上升速度是非线性相关,对气泡的上升速度进行迭代是必要的,迭代过程中可以假设地热流体中的气泡上升的初始速度可定为vb=1.7。
地热流体在井筒中的相对粗糙度为:
其中:ξ为地热井筒的绝对粗糙度,d为井筒内径。
所述的相对粗糙度的公式推导过程如下:
确定地热流体的液体分配系数Γ和摩擦系数f
所以根据计算出的分配系数可知地热流体中的气泡最终的上升速度为:
进而可知地热流体中水的雷诺数为
通过计算得到Nw,所以地热流体在井筒中的相对粗糙度为:
其中:ξ为地热井筒的绝对粗糙度,d为井筒内径。
步骤六、地热流体沿井筒的速度可以在短距离内发生显著变化,并导致可变的地热流体与井筒之间的摩擦损失。所以需要计算不同地热流体的雷诺数下的地热流体与井筒之间的摩擦力。根据公式(18)和(25)可知地热流体的雷诺数和地热井筒的相对粗糙度,计算与一次增压井深位置处对应的地热流体与井筒之间的摩擦力f和地热流体的平均密度具体公式为:
当NRe<2400,
当NRe≥2400,所以地热流体的平均密度/>为:
计算地热流体的摩擦损失梯度τf为:
其中:gc为引力常量;
步骤七、计算地热流体在井筒中的闪蒸深度(即一次增压井深位置处)D,一般地热流体在井筒中发生闪蒸之后还可以发生第二次闪蒸,但是本方法不考虑这种情况,
如:假设井口处的深度为0m,则此时计算出的结果符合所处的气泡流流型边界条件(即表2的条件)内时,此处的深度值即为发生气泡流的深度值,即为地热流体发生闪蒸的位置:
步骤八,最后考虑到地热井的结垢位置通常发生在潜水泵上方,采用潜水泵深度作为公式(31)的分母,计算地热流体预测结垢位置的相对误差为:
其中:ψ1为预测值,即通过上述步骤计算出的地热流体在井筒中的闪蒸深度D。
ψ2为实际值,即通过实际测量得到的地热流体在井筒中的结垢位置。
ω为潜水泵深度,即在地热井筒中潜水泵下潜的深度,可通过实地地热井测量得到。
实施例1
某县第一口地热井的井口参数是:
根据上表中某县第一口地热中的井口参数,将该地热井的井口参数代入具体实施方式中建立的数学模型中计算,可以预测得到地热流体在地热井筒中发生闪蒸的位置为6.3m,根据实际的结垢位置0m来计算其相对误差为5.7%。
实施例2
某县第二口地热井的井口参数是:
根据上表中某县第二口地热中的井口参数,将该地热井的井口参数代入具体实施方式中建立的数学模型中计算,可以预测得到地热流体在地热井筒中发生闪蒸的位置为34.6m,根据实际的结垢位置30m来计算其相对误差为1.5%。
值得注意的是,尽管上面对本发明的优选实施例进行了详细描述,但本发明并不仅仅局限于上述的具体实施方式,上述的实施方式仅仅是示意性的,并不是限制性的,本领域的相关技术人员受本发明的启示,在不脱离本发明宗旨和权利要求保护范围的情况下,还可以做出很多形式,这些均属于本发明的保护范围。
Claims (1)
1.预测地热井筒中结垢位置的数学模型的构建方法,其特征在于包括以下步骤:
步骤一、测量地热井中井口处的地热流体的流量、温度、压力以及流体特性;所述的流体特性包括井口处的地热流体的气体溶解度、标准临界压力、临界温度以及黏度,其中:
气体溶解度:Rs=20.481;
标准临界压力:ppe=22.055kPa;
临界温度:Tpe=374.15;
黏度:μ=1.0×10-3Pa·s;
步骤二、通过计算得到地热流体中的气体压缩系数,具体过程如下:
步骤201,根据井口处的地热流体的压力值设定井深范围内的压力梯度的初始值,然后计算得到一次增压压力值p1,与一次增压压力值对应的距离井口一定距离的位置记为一次增压井深位置;
式中p0为井口处的压力值,Δp为假设的压力梯度的初始值;
步骤202,利用对比温度和对比压力矫正井口处的地热流体的温度和压力的测量值,计算得到井口处地热流体的无量纲温度Tr和无量纲压降pr;
无量纲温度:
无量纲压降:
式中:为井口处的地热流体温度测量值,Tpe为地热流体的临界温度,/>为井口处地热流体压降测量值,ppe为地热流体的临界压力;
步骤203,根据井口处地热流体的无量纲温度Tr和无量纲压降pr,在无量纲气体压缩系数图中查找得到井口处的地热流体中的气体压缩系数Z;
步骤三、修正井口处的地热流体中的液体体积流量和气体体积流量;然后利用根据经验参数修正地热流体的液体质量流量和气体质量流量;最后再计算地热流体中修正后的液相密度和气相密度,具体公式如下:
地热流体在井口处根据经验参数修正后的液体体积流量为:
qL=6.49×10-5q0
其中,qL为地热流体在井口处根据经验参数修正后的液体体积流量,q0为地热流体井口处的实际体积流量;
地热流体在井口处根据经验参数修正后的气体质量流量为:
其中,qg为根据经验参数修正后的井口处地热流体中气体体积流量;R是气体溶解度,z为井口处的地热流体中的气体压缩系数;
qt=qL+qg
式中:qt为根据经验参数修正后的地热流体在井口处的地热流体总的体积流量;
修正后的井口处的地热流体中的液相质量流量wL为:
wL=q0(4.05×10-8γ0+8.85×10-7γgRs)
其中:γ0为地热流体中水对油的流体比重,γg为地热流体中的气含率;
修正后的井口处的地热流体中的气体质量流量wg和地热流体中总的液体质量流量wt为:
wg=8.85×10-7q0γg(R-Rs)
wt=wL+wg
修正后的井口处的地热流体中的液相密度ρL和气相密度ρg分别为:
步骤四、确定在计算出地热流体所处的深度内所处的流型,具体过程为:
步骤401,计算所述的一次增压井深度处的地热流体的测试变量vt、δ和vgD,所述的测试变量用来确定地热流体所处的边界条件:
式中:vt为一次增压井深位置处的地热流体速度,m/s;vgD为一次增压井深位置处的无量纲气体速度;σ为一次增压井深位置处的流体表面张力,N;Ap为地热井筒的井筒截面积;
边界限制条件为:
式中:(L)B为气泡流-段塞流边界的整体无量纲速度数,dh为地热流体流动的当量直径;(L)S为段塞流-过渡流边界的整体无量纲速度数;(L)M为过渡流-雾状流边界的整体无量纲速度数;
步骤402,分别比较δ,vgD,(L)S和(L)M;进行如下判断,若δ<(L)B,则该压力值所对应井深位置处该地热流体所处的流型为气泡流;若δ>(L)B,vgD<(L)S,则该压力值所对应井深位置处该地热流体所处的流型为弹状流;若(L)M>vgD>(L)S,则该压力值所对应井深位置处该地热流体所处的流型为过渡流;若vgD>(L)M,则该压力值所对应井深位置处该地热流体所处的流型为雾状流;若为δ<(L)B情况,则执行下一步;否则执行步骤201,在原有的压力梯度初始值的基础上将压力梯度值调小;然后执行步骤二-四;
步骤五、确定一次增压井深位置处的地热流体的雷诺数、地热流体中的气泡雷诺数和地热井筒的相对粗糙度;
其中,地热流体的雷诺数NRe:
其中μL为地热流体液体动力粘度系数;
通过对vb进行迭代求解地热流体中的气泡雷诺数Nb,公式如下:
地热流体的平均气体上升速度vbf:
地热流体在井筒中的相对粗糙度为:
其中:ξ为地热井筒的绝对粗糙度,d为井筒内径;
步骤六、计算与一次增压井深位置处对应的地热流体与井筒之间的摩擦力f和地热流体的平均密度具体公式为:
当NRe<2400,
当NRe≥2400,
所以地热流体的平均密度为:
计算地热流体的摩擦损失梯度τf为:
其中:gc为引力常量;
步骤七、计算地热流体在井筒中的闪蒸深度D,即一次增压井深位置处:
步骤八,计算地热流体预测结垢位置的相对误差为:
其中:ψ1为预测值,即通过上述步骤计算出的地热流体在井筒中的闪蒸深度D;
ψ2为实际值,即通过实际测量得到的地热流体在井筒中的结垢位置;
ω为潜水泵深度,即在地热井筒中潜水泵下潜的深度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310405346.9A CN116738113A (zh) | 2023-04-17 | 2023-04-17 | 预测地热井筒中结垢位置的数学模型的构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310405346.9A CN116738113A (zh) | 2023-04-17 | 2023-04-17 | 预测地热井筒中结垢位置的数学模型的构建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116738113A true CN116738113A (zh) | 2023-09-12 |
Family
ID=87914018
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310405346.9A Pending CN116738113A (zh) | 2023-04-17 | 2023-04-17 | 预测地热井筒中结垢位置的数学模型的构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116738113A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454063A (zh) * | 2023-12-26 | 2024-01-26 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
-
2023
- 2023-04-17 CN CN202310405346.9A patent/CN116738113A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117454063A (zh) * | 2023-12-26 | 2024-01-26 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
CN117454063B (zh) * | 2023-12-26 | 2024-03-12 | 西南石油大学 | 一种井筒油-气-水多相流流态判别及持水率计算方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104895560B (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
CN102226499B (zh) | 一种模拟深水油气水混输实验装置 | |
Xie et al. | The influence of backpressure on severe slugging in multiphase flow pipeline-riser systems | |
Sun et al. | Effects of phase transition on gas kick migration in deepwater horizontal drilling | |
Wang et al. | Numerical simulation of two phase flow in wellbores by means of drift flux model and pressure based method | |
CN116738113A (zh) | 预测地热井筒中结垢位置的数学模型的构建方法 | |
Carrizales et al. | Prediction of multiphase flow in pipelines: literature Review | |
Zhou et al. | Optimization of methods for liquid loading prediction in deep condensate gas wells | |
Spesivtsev et al. | Study of the Transient Terrain–Induced and Severe Slugging Problems by Use of the Drift–Flux Model | |
Pickering et al. | The prediction of flows in production risers-truth & myth | |
Mao et al. | Heat transfer mechanism of cold-water pipe in ocean thermal energy conversion system | |
Spesivtsev et al. | Comparison of drift-flux and multi-fluid approaches to modeling of multiphase flow in oil and gas wells | |
CN115831243B (zh) | 一种结垢临界温度压力预测方法 | |
Gudmundsdottir et al. | The wellbore simulator FloWell | |
Francke | Thermo-hydraulic model of the two-phase flow in the brine circuit of a geothermal power plant | |
Champel | Discrepancies in brine density databases at geothermal conditions | |
Gao et al. | Study on the migration of gas kicks in undulating sections of horizontal wells | |
Garber et al. | Internal corrosion rate prediction in pipelines and flowlines using a computer model | |
Liu et al. | Experimental study of flow patterns and pressure drops of heavy oil-water-gas vertical flow | |
Yao et al. | Case study on diagnosis and identify the degree of bottom hole liquid accumulation in double-branch horizontal wells in PCOC | |
Wang et al. | Experimental Study on Friction Coefficient of Gas–Oil Two-Phase Flow in a Large Annulus | |
Hegde et al. | CSMFlow: a multiphase flow modeling tool for hydrates in flow assurance | |
Zheng et al. | Research progress on the calculation model of critical liquid carrying flow of gas well | |
Ogai et al. | Research of pressure gradient in a gas well working with a foaming agent | |
Androulakis et al. | Impact of Liquid Pressure Losses and Solid- Phase, in the Performance of a Three- Phase Flow Air- Lift Pump |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |