CN103558364B - 一种水泥浆塑态体积收缩率的确定方法 - Google Patents
一种水泥浆塑态体积收缩率的确定方法 Download PDFInfo
- Publication number
- CN103558364B CN103558364B CN201310583289.XA CN201310583289A CN103558364B CN 103558364 B CN103558364 B CN 103558364B CN 201310583289 A CN201310583289 A CN 201310583289A CN 103558364 B CN103558364 B CN 103558364B
- Authority
- CN
- China
- Prior art keywords
- cement paste
- grout
- volume shrinkage
- sleeve pipe
- pore
- 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.)
- Active
Links
Landscapes
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
本发明公开了一种水泥浆塑态体积收缩率的确定方法,依次包括以下步骤:(1)获取套管、井眼相关尺寸参数;(2)获取套管、地层岩石相关力学参数;(3)确定油气井相关工况参数;(4)计算储层压力、水泥浆初始液柱压力;(5)计算水泥浆的等温压缩系数;(6)配制水泥浆进行失重实验,并确定水泥浆终凝时间;(7)按相似性对比原则确定实际工况下的孔隙压力下降曲线;(8)计算水泥浆塑态体积收缩率。本发明水泥浆塑态体积收缩率数学模型综合考虑了井下套管、地层的实际受力工况和水泥浆的失重性能,原理可靠,操作简单,能有效确定水泥浆在凝结过程中的塑态体积收缩率变化规律,为水泥浆塑态体积收缩率的确定提供了一种科学的方法。
Description
技术领域
本发明涉及油气田开发领域固井过程中水泥浆塑态体积收缩率的确定方法。
背景技术
水泥浆体积收缩指固井后水泥浆在井底温度、压力条件下,从液态逐渐变成固态过程中、以及后期生产全过程中体积总的改变量。研究表明,水泥浆体积收缩主要是由水化反应引起的,原因在于水化反应所生成的水化物体积小于水和水泥颗粒的混合体积。随着水化反应的进行,液态水泥浆逐渐凝固成固态水泥石需经历液相、胶凝、凝固和硬化等四个阶段。因此,可将水泥浆体积收缩分为终凝前塑态体积收缩和后期水泥石体积收缩两个阶段。水泥浆塑态体积收缩会导致水泥浆柱的孔隙压力下降(即胶凝失重),造成水泥浆气侵;后期水泥石体积收缩则使得水泥环与地层和套管间胶结不良,形成微间隙,为气窜提供通道,造成油气资源的大量流失。已有研究表明,水泥浆体积收缩率最大可达6%左右,其中初凝前收缩很小,小于0.5%;90%以上的收缩主要发生在初凝至终凝这段水化反应最激烈的过渡期内;凝固成水泥石后的收缩仍会缓慢进行,但收缩量很小。因此,研究水泥浆塑态体积收缩对于早期气窜的防止、确保油气井勘探开发安全更有意义。
实验模拟测试是目前获得水泥浆体积收缩率的唯一手段,测试水泥浆体积收缩率的实验模拟装置主要有水泥自收缩测试法、“开环”法、比长仪法、ASTM膨胀测试法、高温高压体积变化率仪等。目前实验模拟测试装置主要用于测试水泥浆凝固后的体积收缩率,即终凝后的水泥石体积收缩率,用于测试水泥浆凝结全过程体积收缩率的仪器较少。另外,由于现有测试仪器存在无法模拟井底温度、压力工况,无法有效隔离加压介质和水泥浆等实验原理、方法的限制,目前有关水泥浆体积收缩的实验测试仍处于探索发展阶段,尚不能准确测量水泥浆的收缩与膨胀,无法满足防气窜研究需求。
发明内容
本发明的目的在于提供一种水泥浆塑态体积收缩率的确定方法,该方法原理可靠,实际操作方便,能有效确定水泥浆在凝固过程的塑态体积收缩率变化规律,为水泥浆防气窜能力评价提供了科学依据,具有广阔的应用前景。
为达到以上技术目的,本发明提供以下技术方案。
基于终凝前套管、地层在井下的实际受力工况,建立水泥浆凝固过程中液塑态“套管-水泥浆-地层”物理模型,并求解水泥浆塑态体积收缩对孔隙压力下降影响的解析关系。利用该关系模型,通过室内失重实验,测试水泥浆的失重压差,结合地层岩石、套管的弹性模量、泊松比等力学参数,计算获得终凝前水泥浆塑态体积收缩率大小。
本发明采用水泥浆塑态体积收缩率数学模型,依据弹性力学基础理论(徐芝纶,弹塑性力学[M],北京:高等教育出版社,2006),来推导出水泥浆塑态体积收缩率与孔隙压力下降的关系,具体推导过程如下:
(1)假设条件
为了建立水泥浆凝固前塑态体积收缩数学模型,首先对水泥浆所处的复杂井下环境及受力状态作如下简化:①套管、地层轴向尺寸远大于径向尺寸,可按照平面应变问题分析二者径向、周向应力应变;②以主力气层附近的水泥浆柱为研究对象,主要分析气层位置处水泥浆孔隙压力下降,近似为恒温问题,忽略候凝热传导和水泥水化反应放热引起的温度变化;③研究井段井眼为规则圆形井眼,套管居中度100%;④水泥浆孔隙压力下降所引起的套管、井壁围岩变形均为微小变形,故将套管、井壁围岩考虑为均质、连续、各向同性的理想弹性材料,候凝期间不发生塑性变形;⑤终凝前水泥浆尚未失去流动性,因此套管、水泥浆、地层彼此紧密相连,无微环隙。
(2)模型建立与求解
“套管-水泥浆-地层”物理模型见附图1,设ri、ro、rw、rf分别表示套管内半径、套管外半径、井眼半径和地层半径,显然地层半径rf→∞。由模型假设可知,套管、地层受力均为轴对称应力,为了简化计算,选择极坐标进行分析。由弹性力学理论知,轴对称应力状态下的应力分量半定解和位移分量分别见如下。
上式中A、B、C、H、I、K均是任意常数。
①套管位移变化分析
套管受力状态见附图2,其中qi为候凝时由套管内钻井液或候凝液产生的均布载荷,q0为候凝时由环空水泥浆产生的液柱压力。
显然,应力分布应当是轴对称的。套管受力边界条件为:
将边界条件带入应力分量半定解得:
考虑位移单值条件,套管上同一点不可能有不同的周向位移,因此必须满足B=0,求解上式得:
将上式带入应力半定解得拉梅解答:
当水泥浆孔隙压力下降后,套管所受的外挤压力qo下降,内压qi不变,径向应力、周向应力变化值计算如下:
将A、B、C带入位移分量,并考虑到孔隙压力下降,并令r=ro,可得固井第一界面处套管径向位移变化:
②井壁位移变化分析
“地层——井眼”模型见附图3。由假设条件①可知,地层可看作是半径无穷大的圆,因此地层无穷远处可认为应力为零。由于水泥浆未凝固,故井壁仅受水泥浆液柱压力qo的内压作用。
同理,在套管位移分析的基础上,取ri=rw、ro=rf→∞、qi=qo、qo=0可得:
井壁应力为:
将A、B、C带入位移分量,并考虑到孔隙压力下降,并令r=rw,可得固井第二界面处井壁岩石径向位移变化:
③孔隙压降下的体积膨胀分析
水泥浆具备一定的可压缩性,因孔隙压力变化引起的水泥浆自身体积膨胀量为:
计算水泥浆等温压缩系数Ccem时,认为水泥浆中固相颗粒不可压缩,而水相可以压缩。首先根据Brill-Begges提出的经典模型,计算高温高压下水的等温压缩系数:
然后根据等温压缩系数定义,得:
忽略水泥水化对固相含量的影响,令水固比等于w/s,根据水相体积Vw与水泥浆体积对应关系,得到水泥浆等温压缩系数:
④计算水泥浆体积收缩率
根据模型假设⑤套管、水泥浆、地层紧密相连无间隙,且根据假设②忽略热膨胀效应,因此可得:
Δvsh=Δvcase+Δvf-Δvc
即得:
故水泥浆塑态体积收缩率γ为:
式中:
ri、ro、rw、rf——分别表示套管内半径、套管外半径、井眼半径和地层半径,mm;
A、B、C、H、I、K——均为任意常数,无量纲;
σr、——分别表示径向应力、周向应力,Pa;
——分别表示σr、对应面的切应力,Pa;
qi、qo——分别表示候凝时由套管内钻井液或候凝液产生的均布载荷、候凝时由环空水泥浆产生的液柱压力,Pa;
Δσr、——分别表示径向应力、周向应力变化值,Pa;
Δppore——水泥浆孔隙压力变化值,Pa。
μf、μcase——分别表示地层、套管的泊松比,无量纲;
Ef、Ecase——分别表示地层、套管的弹性模量,Pa;
Δuf、Δucase——分别表示地层、套管的径向位移,mm;
Ccem、Cw——分别表示水泥浆、水的等温压缩系数,MPa-1;
dw、do——分别表示井眼直径、套管外径,mm;
Δh——水泥浆柱长度,m。
P——地层压力,MPa;
T——地层温度,℃。
w/s——水固比,无量纲;
Vw——水相体积,m3;
Δvsh、Δvcase、Δvf、Δvc——分别表示水泥浆塑态体积收缩量、套管膨胀恢复体积、井壁地层岩石膨胀恢复体积、水泥浆自身体积膨胀量,m3;
γ——水泥浆塑态体积收缩率,%。
一种水泥浆塑态体积收缩率的确定方法,依次包括以下步骤:
(1)获取套管、井眼相关尺寸参数:包括套管内半径ri、外半径ro(di=2ri),井眼半径rw(dw=2rw);
(2)获取套管、地层岩石相关力学参数:包括套管的弹性模量Ecase、泊松比μcase,地层岩石的弹性模量Ef、泊松比μf;
(3)确定油气井相关工况参数:包括井深H,地层压力系数PT,电测井温T,钻井液密度ρd,水泥浆密度ρs、水固比w/s;
(4)结合上述实际工况,按公式Pf=PT·H计算储层压力Pf;按公式Ps0=ρs·g·H计算水泥浆初始液柱压力Ps0;
(5)根据Brill-Begges经典模型,按下列公式计算水泥浆的等温压缩系数Ccem:
(6)配制水泥浆进行失重实验,实验模拟环空尺寸为rw-ro,取模拟浆柱为1m,实测模拟井筒环空底部浆柱压力Ps′随时间的变化规律:Ps′(t)~t;失重实验的同时,同步用维卡仪测试确定水泥浆的终凝时间t终凝;
(7)计算相似比按相似性对比原则:Ps(t)=Ps′(t)·H,确定实际工况下井底浆柱压力随时间的变化规律:Ps(t)~t,进而确定实际工况下孔隙压力下降曲线:ΔPpore(t)=Ps0-Ps(t);
(8)最后,按如下水泥浆塑态体积收缩率解析计算模型计算水泥浆体积收缩率随时间的变化关系:γ(t)~t;并由步骤6所测定的终凝时间t终凝确定终凝时的水泥浆塑态体积收缩率γ(t终凝):
与现有技术相比,本发明具有以下有益效果:水泥浆塑态体积收缩率数学模型综合考虑了井下套管、地层的实际受力工况和水泥浆的失重性能,原理可靠,操作简单,能有效确定水泥浆在凝结过程中的塑态体积收缩率变化规律,为水泥浆塑态体积收缩率的确定提供了一种科学的方法,为固井早期气窜评价提供了依据,可有效指导水泥浆防气窜性能评价。
附图说明
图1是“套管-水泥浆-地层”物理模型;
图2是套管受力状态;
图3是井壁受力状态;
图4是实测水泥浆(2.35g/cm3)孔隙压力下降曲线;
图5是实际工况下水泥浆(2.35g/cm3)孔隙压力下降曲线;
图6是水泥浆(2.35g/cm3)体积收缩率发展曲线。
具体实施方式
下面结合附图进一步说明本发明。
一种水泥浆塑态体积收缩率的确定方法,依次包括以下步骤:
(1)获取套管、井眼相关尺寸参数:套管内半径ri=78.55mm、外半径ro=88.9mm(di=177.8mm),井眼半径rw=107.95mm(dw=215.9mm);
(2)获取套管、地层岩石相关力学参数:套管的弹性模量Ecase=206.9GPa、泊松比μcase=0.3,地层岩石的弹性模量Ef=13.84GPa、泊松比μf=0.287;
(3)确定油气井相关工况参数:包括井深H=5085m,地层压力系数PT=0.0225MPa/m,电测井温T=122℃,钻井液密度ρd=2.30g/cm3,水泥浆密度ρs=2.35g/cm3、水固比w/s=0.263;
(4)结合上述实际工况,计算储层压力Pf=114.41MPa、水泥浆初始液柱压力Ps0=114.95MPa;
(5)计算水泥浆的等温压缩系数Ccem=3.68*10-4MPa-1;
(6)配制水泥浆进行失重实验,实验模拟环空尺寸为2cm,取模拟浆柱为1m,实测模拟井筒环空底部浆柱压力Ps′随时间的变化规律:Ps′(t)~t,见图4;失重实验的同时,同步用维卡仪测试确定水泥浆的终凝时间t终凝=487min;
(7)计算相似比H=5085,按相似性对比确定实际工况下井底浆柱压力随时间的变化规律:Ps(t)~t,见图5,进而确定实际工况下孔隙压力下降曲线:ΔPpore(t)=Ps0-Ps(t);
(8)计算水泥浆体积收缩率随时间的变化关系:γ(t)~t,见图6;并由步骤6所测定的终凝时间t终凝=487确定终凝时的水泥浆塑态体积收缩率γ(t终凝)=3.64%。
参看图4、图5、图6。
水泥浆注入环空后,随着水化反应的进行加剧,水泥浆孔隙水被不断地消耗,孔隙压力不断降低,对应的水泥浆塑态体积收缩率随时间不断增大。终凝时,水泥浆塑态体积收缩率高达3.64%,远大于1%,因此该井固井气窜风险特别大。施工过程顺利,未发生井下复杂事故,固井质量解释成果却令人失望。封固段4400~5080m固井质量非常差,综合解释结果合格的井段仅占8%,并在固井后数天检测到了套压异常,表明发生较严重气窜现象。
Claims (1)
1.一种水泥浆塑态体积收缩率的确定方法,依次包括以下步骤:
(1)获取套管、井眼相关尺寸参数:包括套管内半径ri、外半径ro,井眼半径rw;
(2)获取套管、地层岩石相关力学参数:包括套管的弹性模量Ecase、泊松比μcase,地层岩石的弹性模量Ef、泊松比μf;
(3)确定油气井相关工况参数:包括井深H,地层压力系数PT,电测井温T,钻井液密度ρd,水泥浆密度ρs、水固比w/s;
(4)按公式Pf=PT·H计算储层压力Pf;按公式Ps0=ρs·g·H计算水泥浆初始液柱压力Ps0;
(5)根据Brill-Begges经典模型,按下列公式计算水泥浆的等温压缩系数Ccem:
(6)配制水泥浆进行失重实验,实验模拟环空尺寸为rw-ro,取模拟浆柱为1m,实测模拟井筒环空底部浆柱压力Ps′随时间的变化规律:Ps′(t)~t;失重实验的同时,同步用维卡仪测试确定水泥浆的终凝时间t终凝;
(7)计算相似比按相似性对比原则:Ps(t)=Ps′(t)·H,确定实际工况下井底浆柱压力随时间的变化规律:Ps(t)~t,进而确定实际工况下孔隙压力下降曲线:ΔPpore(t)=Ps0-Ps(t);
(8)按如下水泥浆塑态体积收缩率解析计算模型计算水泥浆体积收缩率随时间的变化关系:γ(t)~t;并由步骤6所测定的终凝时间t终凝确定终凝时的水泥浆塑态体积收缩率:
式中:
ri、ro、rw——分别表示套管内半径、套管外半径、井眼半径,mm;
A、B、C——均为任意常数,无量纲;
Δppore——水泥浆孔隙压力变化值,Pa;
μf、μcase——分别表示地层、套管的泊松比,无量纲;
Ef、Ecase——分别表示地层、套管的弹性模量,Pa;
Ccem、Cw——分别表示水泥浆、水的等温压缩系数,MPa-1;
dw、do——分别表示井眼直径、套管外径,mm;
P——地层压力,MPa;
T——地层温度,℃;
w/s——水固比,无量纲;
γ——水泥浆塑态体积收缩率,%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310583289.XA CN103558364B (zh) | 2013-11-19 | 2013-11-19 | 一种水泥浆塑态体积收缩率的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310583289.XA CN103558364B (zh) | 2013-11-19 | 2013-11-19 | 一种水泥浆塑态体积收缩率的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103558364A CN103558364A (zh) | 2014-02-05 |
CN103558364B true CN103558364B (zh) | 2015-03-18 |
Family
ID=50012671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310583289.XA Active CN103558364B (zh) | 2013-11-19 | 2013-11-19 | 一种水泥浆塑态体积收缩率的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103558364B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108533209B (zh) * | 2018-04-02 | 2020-07-10 | 中国石油天然气股份有限公司 | 固井水泥的体积稳定性的确定装置和方法 |
CN110500082B (zh) * | 2018-05-18 | 2021-11-02 | 中国石油化工股份有限公司 | 一种固井水泥浆气侵危险时间的确定方法 |
CN110608967A (zh) * | 2019-09-20 | 2019-12-24 | 黔南民族师范学院 | 一种陶瓷伸缩率的测定方法 |
CN111208281B (zh) * | 2020-01-17 | 2022-06-24 | 西南石油大学 | 水泥水化过程中液塑固三态体积变化的测试装置及方法 |
CN111859701A (zh) * | 2020-07-29 | 2020-10-30 | 西南石油大学 | 一种固井水泥浆塑性体积收缩与孔隙压降评价方法 |
CN112253086B (zh) * | 2020-10-15 | 2022-04-12 | 中国石油大学(华东) | 一种固井初始作用力测量装置及方法 |
CN112782002B (zh) * | 2021-03-01 | 2022-04-22 | 西南石油大学 | 一种固井水泥浆侯凝过程中金属套管变形量测试装置 |
CN113960294B (zh) * | 2021-10-11 | 2023-12-29 | 中冶武汉冶金建筑研究院有限公司 | 一种测试耐火泥浆在砖缝中收缩率的装置及方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5246113B2 (ja) * | 2009-03-19 | 2013-07-24 | 住友大阪セメント株式会社 | コンクリート乾燥収縮ひずみの早期推定方法 |
CN102928578B (zh) * | 2011-08-11 | 2015-03-18 | 中国石油化工股份有限公司 | 一种油井水泥高温高压体积膨胀收缩测试仪 |
CN102590483B (zh) * | 2011-12-29 | 2015-12-16 | 江苏苏博特新材料股份有限公司 | 水泥基材料凝结时间的测试方法 |
JP2015507184A (ja) * | 2011-12-30 | 2015-03-05 | 江蘇博特新材料有限公司Jiangsu Bote New Materials Co.,Ltd | セメント基材料凝結時間のテスト方法 |
-
2013
- 2013-11-19 CN CN201310583289.XA patent/CN103558364B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN103558364A (zh) | 2014-02-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103558364B (zh) | 一种水泥浆塑态体积收缩率的确定方法 | |
CN108868748B (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
Cui et al. | A new approach to DEM simulation of sand production | |
Sampath et al. | Theoretical overview of hydraulic fracturing break-down pressure | |
Liu et al. | Effect of local loads on shale gas well integrity during hydraulic fracturing process | |
CN103206206B (zh) | 用于计算超深井地层破裂压力的方法 | |
Wang et al. | Adaptive finite element–discrete element analysis for the multistage supercritical CO2 fracturing and microseismic modelling of horizontal wells in tight reservoirs considering pre-existing fractures and thermal-hydro-mechanical coupling | |
Cheng et al. | The influence of hydraulic fractures on oil recovery by water flooding processes in tight oil reservoirs: An experimental and numerical approach | |
CN108590601A (zh) | 一种注水扩容膨胀施工参数优选的实验方法 | |
Jimenez et al. | Innovation of annular sealants during the past decades and their direct relationship with on/offshore wellbore economics | |
Li et al. | Deformation and damage of cement sheath in gas storage wells under cyclic loading | |
CN109522588A (zh) | 一种油井水泥环的水-热-力学多场耦合模型的建立及计算方法 | |
Han et al. | Study on the whole life cycle integrity of cement interface in heavy oil thermal recovery well under circulating high temperature condition | |
CN105403444B (zh) | 一种制备高温高压环境下固井水泥石试样的装置及方法 | |
Nie et al. | Continuous-discontinuous element method for three-dimensional thermal cracking of rocks | |
Yang et al. | The stress sensitivity of coal bed methane wells and impact on production | |
Arias | Use of finite-element analysis to improve well cementing in HTHP Conditions | |
CN110671088B (zh) | 一种考虑固相封堵主裂缝的新裂缝起裂压力预测方法 | |
Cai et al. | Numerical study on the evolution of mesoscopic properties and permeability in sandstone under hydromechanical coupling conditions involving industrial internet of things | |
CN111859701A (zh) | 一种固井水泥浆塑性体积收缩与孔隙压降评价方法 | |
Zhao et al. | Salt loading on casing in cased wellbore sections | |
Ma et al. | Investigation on cracking propagation patterns in hydraulic fracturing under different vug distributions: Experimental and numerical case | |
Yang et al. | Kinematical analysis of highway tunnel collapse using nonlinear failure criterion | |
Molina et al. | A computational fluid dynamics approach to predict pressure drop and flow behavior in the near wellbore region of a frac-packed gas well | |
Zhang et al. | Numerical performances of invariable and moving boundary methods during fluid penetration into anisotropic porous media |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |