基于多井壓裂數(shù)據(jù)的煤巖水平地應力反演模型建立及應用

袁 征 王杏尊 黃 杰 郭布民 邢云龍 王緒性

(中海油田服務(wù)股份有限公司油田生產(chǎn)研究院，天津 300450)

基于多井壓裂數(shù)據(jù)的煤巖水平地應力反演模型建立及應用

袁 征 王杏尊 黃 杰 郭布民 邢云龍 王緒性

(中海油田服務(wù)股份有限公司油田生產(chǎn)研究院，天津 300450)

通過分析不同坐標系下井壁應力分布狀態(tài)，結(jié)合巖石起裂準則推導出裂縫破裂壓力的具體形式，在此基礎(chǔ)上建立了破裂壓力與水平主應力的超定方程組；通過對方程組在矩陣形式下分離變量、求導等，最終建立了水平地應力反演模型。基于此模型對柿莊南某區(qū)塊的水平地應力值進行了計算，結(jié)果顯示：本區(qū)最大、最小水平主應力與上覆巖層壓力值相近，三者應力差不超過3MPa，且最小水平主應力最弱，反演結(jié)果與裂縫監(jiān)測數(shù)據(jù)相符。

水力壓裂 破裂壓力 水平地應力 反演模型 超定方程組

1 模型建立

建立煤層氣井井眼坐標系OXYZ，如圖1所示。為建立井壁應力坐標系與地層三向主應力之間的聯(lián)系，在沿井軸方向上的某點為原點建立最大、最小主應力及上覆巖層壓力的0123坐標系。若井眼以最大水平主應力標定的方位角為α,井斜角為β，則兩坐標系存在如下幾何關(guān)系：首先沿上覆巖層壓力旋轉(zhuǎn)方位角α，再將上覆巖層壓力向井軸方向旋轉(zhuǎn)β。

圖1 井眼應力坐標系轉(zhuǎn)換關(guān)系示意圖

坐標系0123與坐標系OXYZ之間的轉(zhuǎn)換關(guān)系可以通過下列矩陣變換得到：

(1)

進而可以得到最大水平主應力σH、最小水平主應力σh與X、Y方向上應力分量關(guān)系式如下：

(2)

式中，σx為X方向的應力分量；σy為Y方向的應力分量；σv為巖石上覆巖層壓力。

根據(jù)坐標系之間的轉(zhuǎn)換關(guān)系，井壁巖石柱坐標系下的周向應力在直角坐標系下的轉(zhuǎn)換關(guān)系如下：

σθθ=(σx+σy)-2(σx-σy)cos2θ-4τxysin2θ-Pb

(3)

式中，σθθ為沿井壁方向的周向應力；τxy為XY面的剪應力；Pb為井底靜液柱壓力。

考慮儲層內(nèi)部孔隙壓力的影響，同時忽略剪應力對周向應力作用時，水力裂縫起裂準則：當最大周向拉伸應力大于巖石抗拉強度時裂縫開始起裂，裂縫起裂公式如下：

-σt=σθθ-pp

(4)

結(jié)合式(3)、(4)，考慮到不同井眼井斜角和方位角，0123坐標系在OXYZ坐標系內(nèi)的應力分量顯然不同，分類整理，可得到水力壓裂井眼起裂方程的兩種形式，即Kirch簡化方程如下：

(5)

式中，pp為儲層孔隙壓力；pf為地層起裂最低井底施工壓力。

將式(2)帶入到式(5)。分離變量整理可得

(6)

分析上述公式可以看出，破裂壓裂、地層孔隙壓力、巖石抗拉強度、井斜角、方位角以及上覆巖層壓力參數(shù)均可以從測井、壓裂施工中獲取，僅最大、最小水平主應力為未知參數(shù)。上述公式可簡寫為：

p=a·σH+b·σh

(7)

顯然存在兩組數(shù)據(jù)可即可得出最大、最小主應力，受測井及施工參數(shù)的錄取存在一定誤差，計算值與實際值存在一定誤差。顯然可增加同地區(qū)同層位的施工數(shù)據(jù)來降低誤差值。當存在n組數(shù)據(jù)時便可得到如下矩陣形式的超定方程組：

簡寫為P=C·X

(8)

當方程的個數(shù)大于未知數(shù)的個數(shù)時將演變成超定方程組，并且現(xiàn)場施工數(shù)據(jù)很慢滿足方程組相容性，然而應力最優(yōu)解應使誤差最小。借鑒最小二乘法思想，定義矩陣方程的方差如下：

e2=(P-C·X)T(P-C·X)

(9)

根據(jù)矩陣知識上式對XT、X求導具有等價性，對上式求導，方差在最小點時導數(shù)值為0，結(jié)果如下：

(10)

進而得：

X=(CTC)-1·CT·P

考慮到深度的影響，在計算過程中應力用當量密度表示。同時考慮到實際方位角與本文出現(xiàn)的方位角存在一定的差別，在計算過程在首先要給出最大或者最小水平地應力的方向，并結(jié)合實際井斜角φ得到上述公式中標定的井斜角α。

通過上述模型可以看出，在計算之前無法判斷應力分量σx、σy的大小。因此每組參數(shù)都會得出兩個形式的方程，當存在n組數(shù)據(jù)時，顯然可以得出2n個形式相同的方程組，進而可以得出2n組解，其中可使方差最小的解可認為最優(yōu)解。

2 沁水盆地南某區(qū)塊地應力計算

沁水盆地為我國煤層氣勘探開發(fā)的主力區(qū)，水裂壓裂施工為認識和分析煤層氣產(chǎn)層提供了豐富的數(shù)據(jù)。本文選取沁水盆地南部柿莊南某區(qū)塊的6口井為例，計算此地區(qū)的地應力值。6口煤層氣井的靶點間距在300m左右，埋深在800～900m之間，屬于二疊系下統(tǒng)山西組。裂縫監(jiān)測結(jié)果顯示最大主應力方向與正北方向的夾角為75°，煤巖抗拉強度1MPa，具體參數(shù)見表1。

表1 柿莊南6口井壓裂基本參數(shù)表

表2 最大、最小水平主應力計算結(jié)果表

3 地應力計算結(jié)果

計算得此地區(qū)最大水平主應力梯度2.14g/cm3、最小水平主應力梯度1.87g/cm3、上覆巖層壓力梯度2.10g/cm3(見表2)。結(jié)果顯示：最小水平主應力小于上覆巖層壓力，最大、最小水平主應力接近，各應力差小于3MPa。對于上述應力狀態(tài)，水力壓裂時易產(chǎn)生垂直裂縫和多裂縫現(xiàn)象，這與裂縫監(jiān)測結(jié)果相一致。

[1] 李志明,張金珠. 地應力與油氣勘探開發(fā)[M].北京：石油工業(yè)出版社，1997:113-121.

[2] 景鋒.我國原位地應力測量與地應力場分析研究進展[J].巖土力學,2011,2(32):51-58.

[3] 張重遠,吳滿路等.地應力測量方法綜述[J]. 河南理工大學學報,2012,3(31):305-311.

[4] 廖永遠,羅東坤,李婉棣等.中國煤層氣開發(fā)戰(zhàn)略[J]. 石油學報,2012,6(33):1098-1102.

[5] 羅陶濤.沁水盆地煤巖儲層特征及壓裂增產(chǎn)措施研究[D].成都理工大學,2010:50-74.

[6] 曹代勇,張守仁,穆宣社等.中國含煤巖系構(gòu)造變形控制因素討論[J]. 中國礦業(yè)大學學報,1999,1(28):26-29.

(責任編輯 劉 馨)

Establishment and Application of Horizontal In-situ Stress Inversion Model in Coal Rock Based on Multiple Fracturing Data

YUAN Zheng, WANG Xingzun, HUANG Jie, GUO Bumin, XING Yunlong, WANG Xuxing

(Production Optimization Department Research Institute, China Oilfield Services Ltd., Tianjin 300450)

An overdetermined equation about fracture pressure and horizontal in-situ stress is established based on analyzing stress distribution on borehole wall and deducing definite formation of fracture pressure. In-situ stress inversion model is given with some settle on the equation set, such as variables separation, derivation. Horizontal in-situ stresses of south Shizhuang is calculated based on the inversion model, and the result shows that the overburden pressure is similar with horizontal in-situ stresses in this area, and difference value under 3MPa among the three stresses. Minimum horizontal in-situ stress value is the smallest, in accordance with crack monitoring data. Keywords:Hydraulic fracture; fracture pressure; horizontal in-situ stress; inversion model; overdetermined equation

袁征，男，壓裂工藝工程師，現(xiàn)從事非常規(guī)致密氣、煤層氣水力壓裂相關(guān)工作。