井筒流壓快速分析方法及其在CO₂井中的應用

井筒作為井筒內工質流體和地下儲層發(fā)生物質和能量交換的主要通道，其內部流體溫度場和壓力場的準確預測是深部流體注采工程中至關重要的內容，控制工程生產效率與安全性評估。通過總結前人研究工作，發(fā)現井筒內流體壓力隨深度近似線性變化的特征，在此基礎上，提出一種關于井筒內流壓的快速分析方法；建立了基于井筒內流體的質量守恒方程、動量守恒方程和能量守恒方程的井筒內流體溫度、壓力耦合計算模型，并利用該耦合模型就所提快速分析方法中待定參數與流體工況、井筒參數、地層物性等參數變化的敏感性進行分析；針對CO₂井，給出了井筒內流壓快速分析表達式，結合工程現場數據進行驗證，壓力解精度在工程應用許可的范圍內。所提井筒內流壓快速分析方法基于文獻統計得出，適用于氣井、液井、生產井、注入井等，覆蓋范圍較廣，具有計算參數少、方便快捷的特點，針對某指定工質流體的快速分析表達式中的待定參數可由室內研究人員給出。研究結果對現場工程師具有極大的便利性。

井筒流動模型；快速分析方法；一維穩(wěn)態(tài)流動；流體壓力；敏感性分析;CO₂井

TE33

A

014

Rapid Analysis Method for Fluid Pressure in

Wellbore and Its Application in CO₂ Wells

He Bangchao¹ Ci Huiling² Bai Bing² Liu Jianfeng³ Zou Yan² Lei Hongwu² Cui Yinxiang²

（1.Jiangsu Changjiang Geological Survey Institute；2.Hubei Key Laboratory of Environmental Geotechnology，Institute of Rock and Soil Mechanics，Chinese Academy of Sciences；3.State Key Laboratory of Hydraulics and Mountain River Engineering，Sichuan University）

Wellbore is the main channel for material and energy exchange between the working fluid in the wellbore and the subsurface reservoir，and the accurate prediction of fluid temperature and pressure fields in the wellbore is a crucial part in deep fluid injection and production engineering，and it conteols the operation efficiency and safety evaluation.Literature review reveals that the fluid pressure in the wellbore changes approximately linearly with depth.On this basis，a rapid analysis method for fluid pressure in the wellbore was proposed.Then，a coupled calculation model for fluid temperature and pressure in a wellbore was built based on the mass conservation equation，momentum conservation equation and energy conservation equation of the fluid in the wellbore，and the coupled model was used to analyze the variation sensitivity of the parameters such as the undetermined parameters，fluid conditions，wellbore parameters and physical properties of formation proposed in the rapid analysis method.Finally，for CO₂ wells，a rapid analysis expression for fluid pressure in the wellbore was presented，and verified with engineering field data，showing that the accuracy of the pressure solution is within the range permitted by engineering applications.The rapid analysis method for fluid pressure in the wellbore proposed in the paper is obtained from the statistics of literatures，suitable for gas wells，liquid wells，production wells and injection wells，having a wide coverage range，and has the characteristics of convenient and fast calculation with fewer parameters.The undetermined parameters in the rapid analysis expression for a specified working fluid can be provided by office researchers.Overall，the proposed method has great convenience for field engineers.

wellbore flow model；rapid analysis method；one-dimensional steady-state flow；fluid pressure；sensitivity analysis;CO₂ well

0 引 言

近年來，深部地下空間和深部能源資源的開發(fā)利用的研究日益受到重視并得到不斷推進，如深部地熱能開采、二氧化碳地質封存、地下鹽穴儲氣庫、油田增滲驅油等^［1-4］。在這類工程中，井筒發(fā)揮了關鍵作用。以地熱能開采為例，井筒作為地表和地下熱儲層發(fā)生物質和能量交換的主要通道，其結構設計不僅會影響生產效率，還會影響工程安全性，如井筒泄漏風險、地面震動響應等。一般來說，影響井筒結構設計的2個關鍵因素包括流體實際壓力和地層允許壓力的比值、流體質量流率（也稱為生產效率）^［5］。

自1962年H.J.RAMEY^［6］提出第一個半穩(wěn)態(tài)井筒傳熱模型，廣大學者已針對井筒內溫度與壓力分布開展了廣泛研究。H.J.RAMEY所提模型雖然是顯式的，但不能用于瞬態(tài)行為^［6-8］，且忽略了摩擦力和重力做功，以及流體相態(tài)的變化；BAI B.等^［9］利用H.J.RAMEY解求解井筒溫度場，進而提出一種快速顯式有限差分方法（Fast Explicit Finite Difference Method，FEFDM）來確定井口或井底的壓力值，其本質是壓力場和溫度場的解耦計算；之后，由LIU M.Z.等^［10］給出了注采井的統一模型，可同時計算注采井的井口和井底壓力；WU H.Q.等^［11］推導并建立了井筒壓力與溫度的雙向耦合微分方程，通過雙FEFDM方法進行求解；李小江等^［8］分析研究了不同熱量源匯對流體流動與傳熱特性的影響，雖然考慮CO₂在工質流體情況下作為浮升力的影響而采用LIAO S.M.等^［12］提出的強迫對流傳熱系數計算公式，但該公式在微小直徑井筒（0.70、1.40、2.16 mm）、低壓力（7.4～12 MPa）情況下得到，其在工程實際中的應用有待考量。近年來，隨著研究的不斷深入，井筒內流體溫壓分析中考慮的因素越來越多，如不穩(wěn)定流動、工質流體相態(tài)變化、瞬態(tài)效應等^{［13-16］。}

雖然前人開展了諸多研究，但對于現場工程師而言，上述方法或公式中計算參數多，且通常涉及迭代計算或耦合計算時，求解井筒內溫度或壓力復雜、耗時，不具備快速分析的特征。因此，為保證結果精度，一種形式簡單的井筒內流體溫度和壓力的計算方法有待提出。通過閱讀文獻發(fā)現，不同工況下的井筒內流體溫度剖面各異，而井筒內流壓多呈現出隨深度線性變化的特征。基于此規(guī)律，本文提出一種關于井筒內流壓的快速分析方法，并具體針對CO₂井，通過流體溫壓耦合分析模型，對快速分析方法中待定參數進行敏感性分析，且進一步給出井筒內流壓的快速計算表達式，再結合工程實例對其進行驗證，即所得到的壓力解的精度是否在工程允許范圍內，從而說明所建立快速分析方法及針對CO₂井的計算表達式的可行性與有效性。

1 井筒流壓快速分析模型

綜合對比分析前人通過數值計算得到的井筒內流體溫度場與壓力場分布以及現場數據^{［5，10，11，13，16-27］}，可得到井筒內流體溫度與流體壓力典型分布形式，如圖1所示。由圖1a可以發(fā)現，在溫度場中，溫度剖面受流體初始溫度、排量、地層特性等諸多因素的影響而各異，但大致可分為3種類型、2個階段。其中流體與地層之間的換熱行為包括如下3種類型：①被注入井筒的低溫流體被迅速加熱，升溫速度逐漸變緩，當達到某臨界深度時，流體溫度隨深度線性增加，即為類型ⅰ；②當被注入的流體為高溫時，首先將發(fā)生高溫流體的降溫過程，當達到某臨界深度時，流體溫度回升，并呈現隨深度線性增長趨勢，即為類型ⅱ；③當流體溫度與地表溫度相差不大時，流體與地層之間的換熱行為不劇烈，流體溫度剖面不存在曲線階段或不明顯，而呈近似線性變化趨勢，即為類型ⅲ。由此引出“臨界深度”與“臨界溫度”的概念（具體見本文2.2節(jié)）。由圖1b可知，壓力場一般隨深度線性變化。基于此分布規(guī)律，提出一種井筒內流體壓力p快速分析方法的一般表達式：

p=Az+B（1）

式中：z為深度，m；A為壓力梯度，Pa/m，A是一待定系數，與流體初始壓力、流體初始溫度等工況信息相關；B在注入井中為流體注入壓力，在生產井中為地層破裂壓力，Pa。

2 參數敏感性分析

為進一步確定待定系數表達式，首先采用高精度的、不忽略摩擦力與重力做功的雙向耦合算法，探討式（1）中待定系數對注入參數、井筒結構參數、地層物性等參數變化的敏感性。

2.1 井筒流體溫壓耦合模型

2.1.1 基本假設

圖2為井筒沿深度方向的示意圖及井筒橫截面的幾何示意圖。圖2中θ為井筒與地面之間夾角，r_ti為油管內半徑，r_to為油管外半徑，r_ci為套管內半徑，r_co為套管外半徑，r_wb為井筒半徑。井筒內流體流動模型基于以下基本假設：①考慮井筒內流體沿深度方向的一維流動；②油管、絕熱層、環(huán)空、套管、水泥環(huán)等與井筒同心；③考慮徑向傳熱，忽略垂直傳熱，井筒內為穩(wěn)態(tài)傳熱，地層為非穩(wěn)態(tài)傳熱。

2.1.2 壓力模型

以注入井為例，任意截取一井筒內流體微元段，其質量守恒方程和運動方程分別為：

式中：z為微元段深度，m；ρ為流體密度，kg/m³；v為流體流速，m/s；g為重力加速度，m/s²；τ為油管壁處的剪切應力，Pa；A_ti為油管內截面積，A_ti=πr_ti²，m²；r_ti為油管內半徑，m。

剪切應力與流速之間的關系式為：

式中：λ為管道摩阻系數，無量綱。

根據式（3）可知ρv=常數，定義參變量C=ρv，進一步可得：

式中：C為質量流速，kg/s。

將式（2）與式（4）代入式（3）中，得到注入井中流體壓力計算式為：

2.1.3 傳熱模型

根據能量守恒原理，可得微元段（上截面為1-1，下截面為2-2）的能量方程：

其中：

式中：E₁、E₂分別為微元段在截面1-1和2-2處的內能，J；p₁、p₂分別為微元所在截面1-1和2-2處的壓力，Pa;V₁、V₂分別為微元段在截面1-1和2-2處流體體積，m³；U、E_k、E_p分別為流體的內能、動能和勢能，J；H為流體比焓，J；V為流體體積，m³；ΔQ_r、ΔW分別為地層向井筒內流體的徑向傳熱量與摩阻損失能量，J;m為流體質量，kg。

流體比焓與流體比熱容之間的關系為：

式中：c_p為流體定壓比熱容，J/（kg·K）。

當假設所截取微元段足夠小而假定壓力均勻不變時，式（8）中等式右側第2項為0。

地層與井筒之間傳遞的熱量為：

式中：U_t為表征地層換熱的總傳熱系數，W/（m²·K）；f_t為地層非穩(wěn)態(tài)傳熱的無因次時間函數；k_e為地層導熱系數，W/（m·K）；T_e為z深度處的地層溫度，K；T為z深度處井筒內流體溫度，K。U_t、f_t的詳細表達式可見文獻［11］。

摩阻損失能量的計算式為：

將式（4）代入式（11），之后將式（9）、式（10）、式（11）代入到式（8）、式（7），可得到：

最終可建立注入井內流體溫度場計算式為：

類似地，可推導建立采出井內流體溫度與壓力的計算式。

2.1.4 模型參數獲取

井筒內壓力場與溫度場預測結果的精確度，一方面受計算方法（如解耦計算、半耦合計算、全耦合計算）的影響，另一方面受流體物性參數精度的影響。目前，流體物性參數的計算方法有很多，但對于不同流體的計算精度有所差異，因此，針對不同的工質流體，應選擇最為適用的流體物性參數計算模型。以CO₂為例，可利用Span-Wagner狀態(tài)方程^［28］計算CO₂流體的密度、比熱容等熱力學參數，采用Fenghour和Vesovic輸運定理^［29-30］計算CO₂流體黏度等輸運性質，所采用公式在計算過程中均考慮了流體相態(tài)的變化。此外，式（4）中管道摩阻系數、式（10）中總傳熱系數、無因次時間參數的計算均參照文獻［11］。

2.2 敏感參數

2.2.1 工程作業(yè)參數

圖3為CO₂流體在不同注入溫度工況下由雙向耦合方法計算得到的流體溫度與壓力，詳細的井底壓力值見表1。流體溫度與壓力剖面符合圖1中典型分布形式。此外，還注意到如下情況：①當超過某一深度（即“臨界深度”）時，流體注入溫度對溫度剖面的影響可忽略不計，且溫度斜率接近地溫梯度；②存在一“臨界注入溫度”，當以該溫度向井筒內注入工質流體時，流體溫度變化平行于地溫呈線性增加趨勢。

圖4為CO₂井中壓力梯度A隨流體注入溫度、注入壓力等參數變化的發(fā)展趨勢圖。從圖4可見：隨注入溫度升高，壓力梯度降低，下降幅度約為6 Pa/（m·K），說明流體溫度對壓力梯度的影響不明顯；流體壓力梯度隨流體注入壓力的增加而增加；流體壓力梯度隨注入流體的質量流量先增加后降低，說明存在最佳注入流量，合理的注入流量設計有助于提高作業(yè)效率。相較于低溫流體，高溫流體注入下的流體壓力梯度對質量流量更為敏感；流體壓力梯度與井筒傾角呈正相關關系，與式（12）具有一致性。在所研究的參數區(qū)間內，流體壓力梯度均在10 kPa/m附近，與重力加速度數值接近。

2.2.2 井筒結構參數與地層參數

除注入工況相關參數，井筒內流體溫度、壓力分布特征還受井筒結構參數（包括油管尺寸、套管尺寸、井筒尺寸等）、地層參數（包括地溫梯度、地表溫度、地層導熱系數）等因素的影響，這些參數主要影響井筒內流體與地層間的換熱進程，從而影響壓力剖面。基于以上考量，研究了不同井筒結構參數、地層參數下的壓力剖面，結果仍呈近似線性變化特征（見圖5），具體的流體壓力值列于表2中。從圖5可以看出：井筒結構參數對壓力分布的影響較小，可忽略不計；在其他參數相同的情況下，井底壓力隨地溫梯度和地表溫度的升高而稍有降低，主要是由井底流體溫度差異引起，但變化量很小；地層導數系數對流體壓力的影響微小。總體看來，相關參數對井筒內流體壓力敏感性影響順序為：

地溫梯度≈地表溫度gt;地層導數系數。

3 CO₂井流壓預測方法與應用

3.1 CO₂井流壓預測表達式

經研究，當以CO₂作為注入井的工質流體時，待定參數的表達式為：

式中：S、X、Y均為參變量，S=T_f1+g_T²T_e^1+g_T²，X=max0，T_f-T_e12，Y=maxsignT_f-T_e4，-signT_f-T_e1.5；G_t為流體注入流量，t/d。

式（14）中系數均為試算得到。同樣地，針對采出井等自下而上的流體流動，可得快速壓力計算表達式：

式中：g_T為地溫梯度，K/m；T_f、T_e分別為流體溫度與地層溫度，K；對注入井來說，流體溫度=注入溫度，地層溫度=地表溫度，即T_f=T_inj，T_e=T_surf；ρ由注入參數注入溫度與注入壓力共同決定，針對采出井，因為其井底流體溫度一般未知，且流體溫度小于等于地層溫度，當流體溫度與地層溫度相差較大時，在采出過程中流體依然與高溫巖層發(fā)生換熱，使溫度升高，此時式（15）中ρ由地層溫度T_e和地層破裂壓力p_frac共同決定。

3.2 工程應用實例

3.2.1 草舍油田草8井試驗

該井注入工況、地層熱力學參數及井筒結構參數等如表3所示。井筒內流體溫度壓力沿深度變化如圖6所示。

從圖6可看出，現場實測壓力亦呈現近似線性分布特征，當采用本文所述快速分析模型時，預測井底3 100 m處的壓力值為56.84 MPa，與實測值（56.35 MPa）的相對誤差為0.87%，處于工程可接受范圍內。從圖6還觀察到，雙向耦合方法得到的流體溫度場和H.J.RAMEY解^［6］高度重合，說明H.J.RAMEY方法在此工程中是高度適用的。

3.2.2 紅河油田長8油藏紅河156井CO₂注入先導性試驗

注入井溫壓參數、地層熱力學參數及井筒結構參數等如表4所示。由于缺少部分資料，表4括號中內容為估計值。井筒內流體溫度、壓力沿深度變化如圖7所示。再次證明，通過雙向耦合方法計算得到的流體溫度與文獻［6］中溫度基本重合。根據文獻［31］，當注入壓力為16.5 MPa時，井底壓力將達到地層破裂壓力（32.8 MPa），而經本文提出CO₂井流壓預測表達式及WU H.Q.等^［11］提出的雙FEFDM方法計算得到的井底處壓力值分別為31.47、33.76 MPa，與實際地層破裂壓力的相對誤差為4.05%、2.93%，均在工程可接受范圍（lt;5%）內，但本文方法具有計算參數少、不需要迭代計算的優(yōu)點。

4 研究結果討論

（1）當假設井筒內流體密度不變時，由式（5）亦可推導出如下線性公式：

圖 8為某工況下本文所提快速分析方法（式（13））、式（15）與雙向耦合壓力解法的對比。從圖8可以看出，因式（16）不包含流體流量、地層特性等因素，將密度過度簡化為常數的方法會引起較大的誤差，于工程設計無益。

（2）當高流速情況下，通過本文快速分析方法計算得到的壓力偏大，對工程設計安全，但存在較大的相對誤差，該快速分析方法在高速流動中的應用有待進一步完善。

（3）當工質流體為水、氮氣或其他鉆井液時，井筒內流體壓力亦呈現隨深度近似線性的關系^{［17，23，32］}，故本文所提快速分析方法具有普遍適用性，可進一步向其他工質流體方面擴展，但式（1）中待定系數需室內研究人員試算決定；對于含節(jié)流裝置的地下井筒，井筒內流體壓力分段滿足本文所提的流體壓力隨深度線性變化的特征^［33-34］，故可采用分段描述的方法建立快速分析表達式。不同的是，式（1）中系數B應為分段函數，節(jié)流前后的流體壓力差應是節(jié)流裝置的結構參數與所處深度的函數。

（4）考慮地層的非均質特征，針對復雜多儲層情況下的井筒內流壓分布也開展了案例研究，結果如圖9所示。雙向耦合方法和本文所提快速分析方法計算得到的井底壓力分別為43.90和44.82 MPa。該案例表明：即使是復雜儲層，本文所提關于流體壓力的快速分析方法仍然可以得到滿意的解；在溫度場中，H.J.RAMEY方法與通過雙向耦合方法得到的解高度重合，H.J.RAMEY方法的高精度特征再次被證明。

5 結 論

（1）考慮當前井筒內流體壓力計算復雜、參數多、耗時的研究現狀，本文總結發(fā)現井筒內流體壓力隨深度近似線性變化的規(guī)律。

（2）基于所發(fā)現規(guī)律，本文提出一種井筒內流壓快速分析方法，并利用經典的、高精度的溫壓雙向耦合方法對比分析了流壓快速分析方法中待定參數對工程作業(yè)參數、井筒結構參數、地層參數等參數變化的敏感性，發(fā)現井筒結構參數與地層參數對待定壓力梯度的影響甚微。在工程作業(yè)參數中，流體注入壓力和質量流量以及井筒傾角對其在井筒內壓力變化率有明顯影響。

（3）以工質流體CO₂為例，給出了具體的井內流體壓力快速預測分析表達式，結合工程應用，驗證了該表達式的可用性，其精度在工程應用許可的范圍內。

（4）本文所提井筒內流壓快速分析方法是根據文獻統計得出，適用于氣井、液井、生產井、注入井等，并可擴展至含節(jié)流裝置的井筒，覆蓋范圍較廣，具有形式簡單、易于理解、計算參數少、方便快捷的特點，針對某指定工質流體的快速分析表達式中的待定參數可由室內研究人員給出，對現場工程師的工作具有極大的便利性。

［1］ RANJITH P G，ZHAO J，JU M H，et al.Opportunities and challenges in deep mining：a brief review［J］.Engineering，2017，3（4）：546-551.

［2］ LU S M.A global review of enhanced geothermal system （EGS）［J］.Renewable and Sustainable Energy Reviews，2018，81，Part 2：2902-2921.

［3］ WANG T T，MA H L，SHI X L，et al.Salt cavern gas storage in an ultra-deep formation in Hubei，China［J］.International Journal of Rock Mechanics and Mining Sciences，2018，102：57-70.

［4］ JIANG K，ASHWORTH P.The development of carbon capture utilization and storage （CCUS） research in China：a bibliometric perspective［J］.Renewable and Sustainable Energy Reviews，2021，138：110521.

［5］ BAI B，LI X C，WU H Q，et al.A methodology for designing maximum allowable wellhead pressure for CO₂ injection：application to the Shenhua CCS demonstration project，China［J］.Greenhouse Gases：Science and Technology，2017，7（1）：158-181.

［6］ RAMEY H J，Jr.Wellbore heat transmission［J］.Journal of Petroleum Technology，1962，14（4）：427-435.

［7］ L X R，ZHANG S C，YU B H，et al.The wellbore temperature test and simulation analysis for liquid carbon dioxide fracturing management［C］∥The SPE/IATMI Asia Pacific Oil amp; Gas Conference and Exhibition.Nusa Dua，Bali，Indonesia：SPE，2015：SPE 176187-MS.

［8］ 李小江，李根生，王海柱，等．超臨界CO₂壓裂井筒流動模型及耦合求解［J］．中國石油大學學報（自然科學版），2018，42（2）：87-94．

LI X J，LI G S，WANG H Z，et al.A wellbore flow model and coupling solution for supercritical CO₂ fracturing［J］.Journal of China University of Petroleum （Edition of Natural Science），2018，42（2）：87-94.

［9］ BAI B，LI X C，LIU M Z，et al.A fast explicit finite difference method for determination of wellhead injection pressure［J］.Journal of Central South University，2012，19（11）：3266-3272.

［10］ LIU M Z，BAI B，LI X C.A unified formula for determination of wellhead pressure and bottom-hole pressure［J］.Energy Procedia，2013，37：3291-3298.

［11］ WU H Q，BAI B，LI X C，et al.An explicit finite difference model for prediction of wellbore pressure and temperature distribution in CO₂ geological sequestration［J］.Greenhouse Gases：Science and Technology，2016，7（2）：353-369.

［12］ LIAO S M，ZHAO T S.An experimental investigation of convection heat transfer to supercritical carbon dioxide in miniature tubes［J］.International Journal of Heat and Mass Transfer，2002，45（25）：5025-5034.

［13］ WANG J T，SUN B J，CHEN W Q，et al.Calculation model of unsteady temperature-pressure fields in wellbores and fractures of supercritical CO₂ fracturing［J］.Fuel，2019，253：1168-1183.

［14］ BALE D S，SATTI R P.A fast computational model for wellbore pressure transients while perforating with dynamic underbalance［J］.Journal of Petroleum Science and Engineering，2020，185：106495.

［15］ 馮劍，許博越，付建紅，等．深水壓井井筒瞬態(tài)傳熱及復雜流動行為研究［J］．石油機械，2021，49（2）：1-6．

FENG J，XU B Y，FU J H，et al.Study on wellbore transient heat transfer and complex flow behavior during deepwater well killing［J］.China Petroleum Machinery，2021，49（2）：1-6.

［16］ 張銳堯，肖平，朱忠喜，等．CML雙梯度鉆井瞬態(tài)井筒溫壓耦合場［J］．石油機械，2023，51（8）：57-63．

ZHANG R Y，XIAO P，ZHU Z X，et al.Transient wellbore temperature-pressure coupling field in CML dual gradient drilling［J］.China Petroleum Machinery，2023，51（8）：57-63.

［17］ LIU X K，LIU L M，YU Z C，et al.Study on the coupling model of wellbore temperature and pressure during the production of high temperature and high pressure gas well［J］.Energy Reports，2022，8（3）：1249-1257.

［18］ WEI M H，WU C H，ZHOU Y X.Study on wellbore temperature and pressure distribution in process of gas hydrate mined by polymer additive CO₂ Jet［J］.Advanced Polymers and Related Composite Materials for Energy Technology，2020，2020：2914375.

［19］ GONG Q，XU Z G，WANG M Q，et al.Numerical investigation on wellbore temperature and pressure during carbon dioxide fracturing［J］.Applied Thermal Engineering，2019，157：113675.

［20］ YANG Z Z，YI L P，LI X G，et al.Model for calculating the wellbore temperature and pressure during supercritical carbon dioxide fracturing in a coalbed methane well［J］.Journal of CO₂ Utilization，2018，26：602-611.

［21］ LYU X，ZHANG S C，MA X F，et al.Numerical investigation of wellbore temperature and pressure fields in CO₂ fracturing［J］.Applied Thermal Engineering，2018，132：760-768.

［22］ 王宴濱，石小磊，高德利，等．深層高溫高壓氣井完井測試管柱失效分析：以順南地區(qū)某井為例［J］．石油鉆采工藝，2022，44（3）：302-308．

WANG Y B，SHI X L，GAO D L，et al.Failure analysis of completion test string for deep high-temperature and high-pressure gas well：a case study on a well in Shunnan area［J］.Oil Drilling amp; Production Technology，2022，44（3）：302-308.

［23］ 王建軍，曾祥俊，賈善坡，等．儲氣庫注采井井筒溫度場預測與影響因素分析［J］．科學技術與工程，2022，22（18）：7890-7902．

WANG J J，ZENG X J，JIA S P，et al.Wellbore temperature field prediction and influencing factors analysis of gas storage injection production wells［J］.Science Technology and Engineering，2022，22（18）：7890-7902.

［24］ 張銳堯，李軍，柳貢慧，等．雙梯度鉆井關鍵工具及井筒壓力動態(tài)變化規(guī)律［J］．石油機械，2022，50（1）：1-9．

ZHANG R Y，LI J，LIU G H，et al.Key tools and dynamic variation of wellbore pressure in dual-gradient drilling［J］.China Petroleum Machinery，2022，50（1）：1-9.

［25］ 吳春方，竇亮彬，劉建坤．CO₂干法壓裂井筒壓力與相態(tài)控制研究［J］．石油機械，2019，47（7）：71-79．

WU C F，DOU L B，LIU J K.Study on wellbore pressure and phase control for CO₂ fracturing［J］.China Petroleum Machinery，2019，47（7）：71-79.

［26］ 蘭建平，龔群，徐治國．CO₂壓裂參數對井內溫度和壓力的影響［J］．石油機械，2018，46（11）：97-103．

LAN J P，GONG Q，XU Z G.Effect of CO₂ fracturing parameters on well temperature and pressure［J］.China Petroleum Machinery，2018，46（11）：97-103.

［27］ 王軻，劉彪，張俊，等．高溫高壓氣井井筒溫度場計算與分析［J］．石油機械，2019，47（1）：8-13．

WANG K，LIU B，ZHANG J，et al.Calculation and analysis of wellbore temperature field in HTHP gas wells［J］.China Petroleum Machinery，2019，47（1）：8-13.

［28］ SPAN R，WAGNER W.A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1 100 K at pressures up to 800 MPa［J］.Journal of Physical and Chemical Reference Data，1996，25（6）：1509-1596.

［29］ FENGHOUR A，WAKEHAM W A，VESOVIC V.The viscosity of carbon dioxide［J］.Journal of Physical and Chemical Reference Data，1998，27（1）：31-39.

［30］ VESOVIC V，WAKEHAM W A，OLCHOWY G A，et al.The transport properties of carbon dioxide［J］.Journal of Physical and Chemical Reference Data，1990，19（3）：763-808.

［31］ 張永剛，羅懿，劉岳龍，等．注CO₂井井筒溫度和壓力分布模型研究及現場應用［J］．巖性油氣藏，2014，26（2）：108-113．

ZHANG Y G，LUO Y，LIU Y L，et al.Research and application of wellbore temperature and pressure distribution models for CO₂ injection well［J］.Lithologic Reservoirs，2014，26（2）：108-113.

［32］ 金永進，林伯韜，王如燕，等．注N₂井井筒溫度壓力耦合下的井底流壓計算［J］．石油鉆采工藝，2018，40（4）：489-493．

JIN Y J，LIN B T，WANG R Y，et al.A calculation model for the bottom hole flowing pressure of nitrogen injector with coupled well temperature and pressure［J］.Oil Drilling amp; Production Technology，2018，40（4）：489-493.

［33］ 張凌筱．井下節(jié)流氣井井底壓力計算方法應用研究［J］．承德石油高等專科學校學報，2022，24（4）：34-39．

ZHANG L X.Application of Bottom hole pressure calculation method for gas well with downhole throttling［J］.Journal of Chengde Petroleum College，2022，24（4）：34-39.

［34］ 朱禮濤．注采工況井筒溫度壓力耦合及井下節(jié)流研究［D］．西安：西安石油大學，2019．

ZHU L T.Research on temperature and pressure coupling and downhole throttling of wellbore under injection and production conditions［D］.Xi’an：Xi’an Shiyou University，2019.第一

何邦超，高級工程師，生于1985年，2008年畢業(yè)于中國礦業(yè)大學地質工程專業(yè)，現從事水文地質、地熱與干熱巖勘察開發(fā)工作。地址：（210046）江蘇省南京市。email：15062262272@163.com。

2024-02-08

楊曉峰