A Logging Data-Based Calculation Method for the Horizontal TIV Formation In-Situ Stress
-
摘要:
低孔低渗地层常具有TIV各向异性特征,准确计算水平地应力对该类储层的射孔和压裂设计至关重要。为了更加准确地计算TIV地层的水平地应力,针对该类地层的低压特征,优选Bowers法求取地层孔隙压力,然后利用阵列声波测井资料求取TIV地层的岩石力学参数(垂直与水平方向上的杨氏弹性模量和泊松比);同时考虑层理面产状对水平地应力的影响,改进了传统的Sn模型,建立了TIV地层水平地应力新计算模型。用新模型计算了鄂尔多斯盆地合水地区长6、长7和长8段TIV地层的水平地应力,计算结果与实测最大、最小水平地应力的最大相对误差分别为8.70%和7.86%,低于Sn模型的相对误差。研究结果表明,新模型计算的水平地应力更符合实际地应力纵横向分布的变化规律,可为致密油储层的水力压裂设计提供更可靠的参考依据。
Abstract:Given the TIV anisotropy characteristics of low porosity/permeability formations, the accurate calculation of horizontal in-situ stress is essential for the perforation and fracturing designs of those reservoirs. In order to more accurately and effectively calculate the horizontal in-situ stress of a TIV formation, considering the low-pressure characteristics of such formation, the Bowers method was used to obtain the formation pore pressure. Then, the array acoustic logging data was used to obtain the anisotropic rock mechanical parameters (the vertical/ horizontal Young’s modulus of elasticity and Poisson's ratio) of TIV formation. Considering the influence of bedding plane occurrence on horizontal in-situ stress, the traditional Sn model was improved and a new calculation model for the horizontal in-situ stress of TIV formation was established. The horizontal in-situ stresses of the TIV formations in the Chang 6, Chang 7 and Chang 8 sections of the Heshui area of the Ordos Basin were calculated by this new model. The maximum relative errors between the calculated and the measured maximum/minimum horizontal in-situ stresses were 8.70% and 7.86%, respectively, which were smaller than those of in-situ stresses calculated by Sn model. The results showed that the horizontal in-situ stress calculated by this new model was more in line with the variation laws of the vertical and horizontal distributions of the actual in-situ stress, which could provide a more reliable reference for the hydraulic fracturing design of tight oil reservoirs.
-
随着石油开采技术的发展,对井下流量测量的要求越来越高。特别是,井下低流量(流量小于5 m3/d)的检测一直是难点。相关的测量方法中,热式质量流量计启动流量小、测量精度高,认为是测量井下微小流量的最佳方法。
目前,热式质量流量计通常用于检测气体流量[1-2],将其用于测量井下液相流量方面的研究国内尚不多见;而且,目相关研究主要集中在恒功率热式流量计方面[3-5]。与恒功率热式流量计相比,恒温差热式流量计功耗低、响应速度快,不仅适合于低流量测量,而且其加热器对被测流体所在环境温度影响小,特别适合阵列化检测结构。不过,由热式流量计测量原理可知,被测流体的物性参数会随流体温度变化而改变,从而影响流量计的输出(恒温差模式下,会影响加热功率;恒功率模式下,会影响反映温差的输出电压),导致计量误差较大[6]。实际应用中,不同深度的井段井温不同,井下环境各异,应用恒温差热式流量计时需考虑井下被测流体温度的变化对测量结果的影响,对其测量结果进行温度校正。
目前,仅有学者研究了热式气体流量计的温度校正算法,并取得了较好的补偿效果[6-7]。为此,笔者以水作为测量流体,分析了温度对液体物性参数的影响,研究了温度对恒温差热式流量计测量结果的影响规律,通过数值模拟分析了井下不同深度处温度、压力对恒温差热式流量计测量结果的影响,为井下恒温差热式流量计测量结果的温度、压力校正提供了理论和试验依据。
1. 恒温差热式流量计测量原理
热式流量计测量的物理基础是热传递。根据传热学知识,热式流量计测量时加热源和流体的热量交换以强迫对流传热为主。流体流动时,强迫对流传热从加热源表面带走的热流量可以表示为[5]:
Φf=hA(th−te) (1) 式中:
Φf 为强迫对流传热从加热源表面带走的热流量,W;h为强迫对流平均传热系数,W/(m2∙℃);A为加热源换热表面积,m2;th和te分别为加热源和流体所在环境的温度,℃。对于长度为l、直径为d的热线式加热源,表面积可表示为:
A=πld (2) 为了确定换热系数与流体物性参数的关系,引入了努塞尔数(Nu)、普朗特数(Pr)和雷诺数(Re)
等3个热力学参数,其表达式分别为: Nu=hdλf (3) Pr=ηCpλf (4) Re=ρvdη (5) 式中:
λf 为被测流体的热导率,W/(m∙℃);η 为流体的动力黏度,Pa∙s;Cp 为流体定压比热容,J/(kg∙℃);ρ 为流体密度,kg/m3;v为流体速度,m/s。当流体流动时,可忽略自然对流换热传热的影响,则加热源产生的热流量
H 等于强迫对流流体带走的热流量Φf 。此时,式(1)可写为:H=πlλfNu(th−te) (6) 式(6)中,努塞尔数
Nu 是表示对流换热强烈程度的参数。学者们对其进行过深入研究,并提出了一些对流换热公式,也在很多场合进行了应用。根据H. A. Kramers[8]给出的对流换热公式,在一定条件下Nu 可以表示为:Nu=0.42Pr0.2+0.57Pr0.33Re0.5 (7) 式(7)中,与流速成正比的雷诺数
Re 的指数会随流体流速变化而发生变化,0.5仅在一定条件下适用。一般情况下,Re 的指数用m代替。因此,综合式(2)—式(7)可以得到:H=πlλf(th−te)[0.42Pr0.2+0.57Pr0.33(ρvdη)m] (8) 令
Ac=0.42πlλfPr0.2,Bc=0.57πlλfPr0.33(ρd/η)m ,并为了便于表示,将式(8)中各物性参数统一为Ac 和Bc ,则式(8)可简化为:H=(Ac+Bcvm)(th−te) (9) 在常压(1标准大气压)下,对于给定的加热源和确定的流体温度,
Ac 和Bc 反映此时被测流体热导率、普朗特数等物性参数的综合计算结果,可以视为常量。根据式(9),若保持th−te 不变,加热源维持温差恒定所产生的热流量H 与流体流速v 存在唯一且单调的关系,用这种通过维持温差恒定测量加热源产生的热流量来计算流体流量的方法被称为恒温差法。根据上述原理,恒温差热式流量计采用了2个相同的温度传感器(一个作为测温探头,测量流体的环境温度;另一个与加热器集成在一起作为测速探头,测量加热器温度)。考虑导热传热对测量有影响,测速探头采用隔热陶瓷,以减少加热器沿探头的热传导。测量时,将测温探头放置在流体的上游,测量流体环境温度;将测速探头放置在流体的下游,测量加热器的温度(见图1)。
电路工作时,测速探头加热器接通电源加热,加热器产生的热流量
H 与加热器的功率P 相关[9]。在热平衡状态下,式(9)可改写为:P=k(Ac+Bcvm)(th−te) (10) 式中:
k 为加热器功率因数。2. 定压条件下温度对流量计的影响
通过数值模拟和室内试验方式,分析了常压(1标准大气压)条件下温度对恒温差热式流量计测量结果的影响情况。
2.1 模拟分析
2.1.1 流体温度
由式(8)和式(9)可见,
Ac 和Bc 由加热源结构和被测流体物性参数所决定,是流体的普朗特系数Pr 、密度ρ 、热导率λf 和动力黏度η 的函数。由于这些参数会随流体温度变化而变化,当流体温度发生变化时,式(10)可改写为:P(te)=k[Ac(te)+Bc(te)vm](th−te) (11) 因此,在不同流体温度下,同样流速的流体会对应不同的输出功率。
根据常压下水的物性参数[10],可得到水温从0 ℃升高至150 ℃时
Pr,ρ,λf 和η 的变化曲线(见图2)。从图2可以看出,当水温由0 ℃升高至150 ℃时,水的物性参数随之变化,在低温区域(<90 ℃)变化尤为显著。由式(11)可知,物性参数的改变会导致恒温差热式流量计测量结果变化。
2.1.2 环境温度
为分析环境温度对强迫对流换热功率的影响,设加热器为长度6.4 cm、直径0.8 cm的圆柱体,理想情况下加热器功率因数k=1.0,流速指数m=0.5,在温差维持0.5 ℃条件下,根据式(11)和图2所示物性参数,利用数值模拟法,分析了环境温度分别为25,30,35和40 ℃ 时水流量与换热功率的关系,结果如图3所示。
由图3可知:环境温度变化会对反映流速的换热功率产生显著影响;相同流速条件下,液体换热功率随着温度升高而增大,主要原因是空气的动力黏度随温度升高而增大,而水的动力黏度随温度升高而降低。模拟结果表明,相同流量下,恒温差热式流量计的测量结果会随环境温度升高而增大。
2.2 试验分析
2.2.1 试验平台
基于恒温差热式流量计的理论模型,为进行等梯度温度试验、验证不同环境温度下流体流速与加热器功率的关系,搭建了试验平台(见图4)。标准流量计用来测量流体的实际流量,并与热式流量计输出信号进行比较;水泵用于调节实际流量,最大扬程6 m,最大流量960 L/h;恒温加热系统由加热带、温度控制器和显示屏组成,将流体温度控制在设定温度,误差不超过±0.1 ℃;模拟井筒直径124.0 mm,进水口与恒温加热系统相连,出水口与蓄水箱相连,在水泵的驱动下,流体在蓄水箱和模拟井筒中循环,模拟井筒中流体自下而上的流动状况。
试验的温差采集和恒温差控制方案如图5所示。其中,测温探头用于测量流体环境温度,测速探头用于测量加热器的温度,均与电路系统相连,电路系统最终将测量结果上传至主机显示、保存。
测量控制系统的基本技术指标:1)恒温差数据采集模块采样间隔为500 ms,采集的温差电压精度为1 μV,温度分辨率可达0.01 ℃;2)数字电源模块最大输出电压为4.8 V,调整最小幅度单位为0.1 mV,输出最大功率为1 W。
2.2.2 试验方式及结果
等梯度温度试验中,进行每组试验时保持环境温度不变,流量以1 m3/d为增量,依次获得1~15 m3/d的标准流量。每个标准流量下,自动调节加热器功率,使测温探头与测速探头的温差维持稳定(两温度传感器温差电压保持在2 mV,偏差不超过±0.1 mV),记录温差稳定后加热器的功率。
按照上述试验方式,依次测量环境温度分别为25,30,35和40 ℃时,不同流量下加热器的功率,结果如图6所示。
由图6可知:1)流量相同条件下,环境温度升高,会导致恒温差热式流量计的输出功率显著提高,这与数值模拟结果基本一致;2)相对于数值模拟结果,试验结果受环境温度的影响更大,低流量(<5 m3/d)时受到的影响更大。分析认为,试验结果受环境温度影响较大的原因是:数值模拟仅考虑了强迫对流换热,忽略了其他形式的换热;实际情况是,其他形式换热也不同程度地存在;另外,流量低时井筒内的流体循环速度变慢,维持井筒流体恒温系统的温度调整会呈现较大程度的延时,导致井筒流体环境温度在设定温度上下发生较大波动,进而影响测量结果。
3. 实测条件下流量计的影响因素分析
实际测量时,恒温差热式流量计要下入到井中,随着垂深增深,井筒中的温度和压力均会升高,流体的物性参数也会受到温度和压力的综合影响,进而影响流量计的测量结果[11-12]。由于条件所限,暂时无法开展现场试验研究,因此只进行了理论分析。
井下特定深度的温度和压力与垂深之间存在线性递增关系,为了研究不同垂深时井筒温度和压力对恒温差热式流量计测量结果的影响,将式(11)改为如下形式:
P(D)=k[Ac(D)+Bc(D)vm](th−te) (12) 式中:D为油井的垂深,m。
取井深温度梯度为3.0 ℃/100m,压力梯度为784 kPa/100m,根据文献[10],参考温度和表压力,计算得到不同垂深下水的物性参数(见表1),并代入式(12),得到不同垂深下水的流量与加热功率的关系,结果如图7所示。
表 1 不同井深条件下水的物性参数Table 1. Physical parameters of water at different well depths井深/
m温度/
℃表压力/
MPa密度/
(kg∙m–3)导热系数/
(W∙(m∙K)–1)动力黏度/
(Pa∙s)普朗
特数0 25 0 997.05 0.607 2 0.890 6.130 22 500 40 3.92 993.93 0.632 4 0.653 4.306 72 1 000 55 7.84 989.09 0.653 0 0.506 3.223 93 1 500 70 11.76 982.89 0.668 9 0.407 2.532 24 2 000 85 15.68 975.57 0.680 9 0.338 2.065 45 2 500 100 19.60 967.30 0.689 8 0.287 1.736 63 3 000 115 23.52 958.21 0.696 1 0.249 1.496 99 3 500 130 27.44 948.39 0.700 4 0.220 1.317 57 4 000 145 31.36 937.90 0.703 0 0.197 1.180 40 由图7可知:1)相同流量条件下,垂深增加,温度和压力同步升高,恒温差热式流量计的输出功率也相应升高,表明井筒温度、压力同步升高导致测量误差增大;2)相同流量条件下,尽管随着温度和压力同步升高,恒温差热式流量计的输出功率升高,但升高的幅度逐步减小。垂深0~2 000 m井段,垂深变化对恒温差热式流量计的输出功率影响很大;而在垂深2 000~4 000 m井段,垂深变化对输出功率影响较小。换言之,若恒温差热式流量计的输出功率为0.5 W,在垂深2 000~4 000 m井段对应最大的流量误差约为0.6 m3/d,在垂深0~2 000 m井段对应最大的流量误差却达到2.2 m3/d。其原因在于:垂深较浅的井段,温度、压力都比较低,此垂深下温度和压力变化引起的物性参数变化率高,变化趋势快,对恒温差热式流量计输出功率影响大;而垂深较深的井段,温度一般在85 ℃以上,对应的物性参数
Pr 、λf 和η的变化较小,且ρ 会随温度升高降低、随压力升高而升高,因此总的变化趋势相对平缓,对流量计输出功率的影响减小。4. 结论与建议
1)井筒温度和压力的变化会影响恒温差热式流量计的输出功率。在相同流量下,恒温差热式流量计的输出功率会随垂深增加而升高;在垂深浅的井段,温度和压力对该流量计输出功率的影响较大,垂深深井段的影响相对减小。
2)为了获得相对准确的流量检测结果,要对恒温差热式流量计的测量结果进行温度和压力校正,特别是需要校对垂深较浅井段的测量结果。
3)建议通过现场试验获取更全面的数据,分析不同垂深条件下、恒温差热式流量计不同流量下的输出功率,建立可靠的校正图版或网络,以对实测数据进行深度校正,从而获得更准确的测量结果。
-
表 1 最大、最小水平地应力实测值与模型计算值的对比
Table 1 Comparison of the measured and calculated values of the maximum/minimum horizontal in-situ stresses
井深/m 最大水平地应力/MPa 最大水平地应力相对误差,% 最小水平地应力/MPa 最小水平地应力相对误差,% 实测 新模型 Sn模型 新模型 Sn模型 实测 新模型 Sn模型 新模型 Sn模型 1 637.00 30.14 32.64 32.61 8.28 9.09 23.88 25.10 20.02 5.11 16.16 1 664.00 36.04 37.06 36.20 2.84 0.44 29.12 28.51 22.31 2.09 23.39 1 676.40 37.89 38.68 32.78 2.08 13.49 31.02 29.76 19.61 4.08 36.78 1 708.90 37.23 40.17 37.87 7.91 1.72 31.01 30.90 23.75 0.35 23.41 1 719.50 38.74 41.35 35.32 6.75 8.83 32.70 33.08 20.62 1.16 36.94 1 720.60 34.10 37.06 35.78 8.70 4.93 27.10 29.23 21.54 7.86 20.52 1 780.90 40.48 41.28 36.72 1.97 9.29 30.58 31.75 20.21 3.83 33.91 1 787.00 33.89 35.78 44.78 5.55 32.13 25.72 27.52 27.49 7.00 6.88 -
[1] 黄鑫,董秀成,肖春跃,等. 非常规油气勘探开发现状及发展前景[J]. 天然气与石油, 2012, 30(6): 38–41. HUANG Xin, DONG Xiucheng, XIAO Chunyue, et al. Present situation and development prospect of unconventional oil and gas exploration and development[J]. Natural Gas and Oil, 2012, 30(6): 38–41.
[2] 王越之,李自俊. 横向各向同性地层地应力的研究[J]. 石油学报, 1999, 34(1): 34–37. doi: 10.3321/j.issn:0253-2697.1999.01.007 WANG Yuezhi, LI Zijun. A study of in-situ stress in transversely isotropic formations[J]. Acta Petrolei Sinica, 1999, 34(1): 34–37. doi: 10.3321/j.issn:0253-2697.1999.01.007
[3] WATERS G A, HEINZE J R, JACKSON R, et al. Use of horizontal well image tools to optimize Barnett shale reservoir exploitation[R]. SPE103202, 2006.
[4] 姜梦奇.横向各向同性介质克希霍夫叠前深度偏移方法研究[D].青岛: 中国石油大学(华东), 2016. JIANG Mengqi. Study on Kirchhoff pre-stack depth migration method for transversely isotropic media[D]. Qingdao: China University of Petroleum(East China), 2016.
[5] 赵向原,曾联波,王晓东,等. 鄂尔多斯盆地宁县—合水地区长6、长7、长8储层裂缝差异性及开发意义[J]. 地质科学, 2015, 50(1): 274–285. doi: 10.3969/j.issn.0563-5020.2015.01.018 ZHAO Xiangyuan, ZENG Lianbo, WANG Xiaodong, et al. Differences of natural fracture characteristics and their development significance in Chang 6, Chang 7 and Chang 8 reservoir, Ningxian-Heshui area, Ordos Basin[J]. Scientia Geological Sinica, 2015, 50(1): 274–285. doi: 10.3969/j.issn.0563-5020.2015.01.018
[6] 刘忠华,宋连腾,王长胜,等. 各向异性快地层最小水平主应力测井计算方法[J]. 石油勘探与开发, 2017, 44(5): 745–752. LIU Zhonghua, SONG Lianteng, WANG Changsheng, et al. Evaluation method of the least horizontal principal stress by logging data in anisotropic fast formations[J]. Petroleum Exploration and Development, 2017, 44(5): 745–752.
[7] 黄荣樽,庄锦江. 一种新的地层破裂压力预测方法[J]. 石油钻采工艺, 1986, 8(3): 1–14. HUANG Rongzun, ZHUANG Jinjiang. A new method for predicting formation fracture pressure[J]. Oil Drilling & Production Technology, 1986, 8(3): 1–14.
[8] HIGGINS S M, GOODWIN S A, DONALD A, et al. Anisotropic stress models improve completion design in the Baxter Shale[R]. SPE 115736, 2008.
[9] 刘方,陈剑军. 地质统计学反演预测地层压力的方法[J]. 工业技术创新, 2015, 2(3): 366–369. LIU Fang, CHEN Jianjun. Formation pressure prediction with geostatistical inversion[J]. Industrial Technology Innovation, 2015, 2(3): 366–369.
[10] 杨振平,吴波,王勇. Eaton法预测M油田地层孔隙压力[J]. 石油天然气学报, 2012, 34(9): 181–182. YANG Zhenping, WU Bo, WANG Yong. Application of Eaton method in predicting the formation pore pressure of M Oilfield[J]. Journal of Oil and Gas Technology, 2012, 34(9): 181–182.
[11] BOWERS G L. Pore pressure estimation from velocity data: accounting for overpressure mechanisms besides undercompaction[J]. SPE Drilling & Completion, 1995, 10(2): 89–95.
[12] 徐志星. 西湖凹陷异常地层压力特征及其与油气成藏的关系[D]. 成都: 成都理工大学, 2015. XU Zhixing. Abnormal formation pressure and its relationship with hydrocarbon accumulation in Xihu Sag[D]. Chengdu: Chengdu University of Technology, 2015.
[13] 夏宏泉,么勃卫,毕刚. 塔里木盆地致密砂岩气层异常高压的测井计算方法[J]. 测井技术, 2018, 42(4): 423–426. XIA Hongquan, YAO Bowei, BI Gang. Logging calculate method of abnormal high pressure in tight sandstone gas reservoirs in Tarim Basin[J]. Well Logging Technology, 2018, 42(4): 423–426.
[14] WANG Zhijing. Seismic anisotropy in sedimentary rocks, part 1: a single-plug laboratory method[J]. Geophysics, 2002, 67(5): 1415–1422. doi: 10.1190/1.1512787
[15] SAROUT J, MOLEZ L, GUEGUEN Y, et al. Shale dynamic pro-perties and anisotropy under triaxial loading: experimental and theoretical investigations[J]. Physics and Chemistry of the Earth, 2007, 32(8–14): 896–906. doi: 10.1016/j.pce.2006.01.007
[16] SUAREZ-RIVERA R, DEENADAYALU C, YANG Y K. Unlock-ing the unconventional oil and gas reservoirs: the effect of laminated heterogeneity in wellbore stability and completion of light gas shale reservoirs[R]. OTC 20269, 2009.
-
期刊类型引用(8)
1. 刘金辉,刘升虎,党瑞荣. 横掠单管热量扩散特性的研究. 西安文理学院学报(自然科学版). 2025(01): 79-84 . 百度学术
2. 冯爽,魏勇,杜雪梅,李冰,刘杰,陈强,林斯. 基于时域积分的温差流量传感器仿真与试验. 石油机械. 2024(04): 110-119 . 百度学术
3. 陈强,刘国权,王志杰,魏宝军,冯爽,魏勇,甘如饴. 井下流量监测系统软件设计. 仪器仪表与分析监测. 2024(04): 25-29 . 百度学术
4. 鲁义攀,魏勇,陈强,刘国权,刘杰. 基于热传导时域积分的井下流量测量方法. 石油钻探技术. 2023(01): 106-114 . 本站查看
5. 谈聪,杨旭辉,刘平,李冰,吕品歧,霍东凯,杨启聪. 基于SOA与模糊PID的恒温差热式流量计. 石油机械. 2023(07): 113-120 . 百度学术
6. 张镨,周理,张佩颖,罗勤,蒲长胜. 天然气管道掺氢对天然气分析计量的影响. 天然气工业. 2023(08): 135-145 . 百度学术
7. 林斯,魏勇,魏宝军,刘杰,魏磊. 基于时域积分的热式流量检测模型与FPGA实现. 传感器与微系统. 2023(09): 85-88 . 百度学术
8. 秦宏伟,党瑞荣,党博,曹峰. 基于物联网的油田井口数据检测传输方法研究. 仪表技术与传感器. 2022(11): 96-100+112 . 百度学术
其他类型引用(9)