賈鎖剛，董長(zhǎng)銀，武 龍，隆佳佳，馮勝利，陳 君，阿雪慶

(1.中國(guó)石油青海油田分公司鉆采院，甘肅敦煌736200;2.中國(guó)石油大學(xué)石油工程學(xué)院，山東青島266580;3.川慶鉆探工程有限公司鉆采工程技術(shù)研究院長(zhǎng)慶分院，陜西西安710018)

水平井礫石充填包括管外地層礫石充填和套管內(nèi)循環(huán)礫石充填，是一種效果好、有效期長(zhǎng)的防砂(完井)工藝。管外礫石充填是指高壓大排量條件下礫石顆粒由攜砂液帶入地層，在管外地層中形成高滲透帶。管外礫石充填技術(shù)不但具有良好的防砂效果，而且還可以最大程度降低對(duì)油氣井產(chǎn)能的傷害，甚至具有增產(chǎn)作用［1-3］。以上優(yōu)點(diǎn)使水平井管外礫石充填成為疏松砂巖水平井礫石充填防砂完井的發(fā)展方向之一，并已有少量實(shí)際應(yīng)用。文獻(xiàn)［4］中提出了疏松砂巖地層水平井管外礫石充填的兩種模式:裂縫開裂延伸和塑性擠壓破壞模式，并研究了管外地層巖石破壞模式的判別方法。對(duì)于水平井地層巖石的破壞研究，重點(diǎn)主要集中在低滲透水平井水力壓裂領(lǐng)域的裂縫開裂條件和裂縫延伸領(lǐng)域［5］，對(duì)疏松砂巖地層在高壓擠注條件下的塑性擠壓充填模擬研究幾乎處于空白。文獻(xiàn)［6］根據(jù)摩爾－庫(kù)侖準(zhǔn)則和平面應(yīng)變軸對(duì)稱問(wèn)題的基本方程，建立了垂直井近井地帶塑性變形區(qū)的應(yīng)力與應(yīng)變關(guān)系數(shù)學(xué)模型。文獻(xiàn)［7］依據(jù)離散元理論，使用商業(yè)模擬軟件PFC2D對(duì)地層巖石按照膠結(jié)強(qiáng)度強(qiáng)、膠結(jié)強(qiáng)度中等、膠結(jié)強(qiáng)度弱等3種情況進(jìn)行模擬，得到一些定性認(rèn)識(shí)。但上述研究均未對(duì)水平井開展針對(duì)性研究，也未形成系統(tǒng)的模擬理論與方法。筆者針對(duì)上述水平井管外地層礫石充填問(wèn)題，基于沿水平井井身軌跡變化的主應(yīng)力及地層巖石力學(xué)參數(shù)非均質(zhì)分布的基礎(chǔ)，建立疏松砂巖地層高壓擠注條件下的水平井管外地層塑性擠壓充填模擬數(shù)學(xué)模型，預(yù)測(cè)充填半徑隨施工時(shí)間的變化及其沿水平井井身軌跡的分布規(guī)律。

1 水平井近井地應(yīng)力分布

1.1 水平井近井地應(yīng)力

假設(shè):在地層深處，高壓擠壓井段圍壓起主要作用，并且井筒直徑遠(yuǎn)小于施工井段長(zhǎng)度，不考慮井眼軸線方向的應(yīng)變。

基于上述假設(shè)，疏松砂巖地層水平井管外地層塑性擠壓充填可簡(jiǎn)化為垂直于井筒軸線的平面應(yīng)變問(wèn)題。使用柱坐標(biāo)系，以具有普遍意義的斜井為研究對(duì)象，原始地層主應(yīng)力可分解為斜井直角坐標(biāo)系下的 σxx、σxy、σyy、σyz、σzz和 σxz6 個(gè)應(yīng)力分量。由于坐標(biāo)轉(zhuǎn)換后井周截面上的最大和最小主應(yīng)力不一定是σxx和σyy，因此，為了得到井周截面上的最大和最小主應(yīng)力，仍需將坐標(biāo)系繞z軸逆時(shí)針旋轉(zhuǎn)(圖1)。在繞井眼軸線旋轉(zhuǎn)過(guò)程中，各個(gè)面上的應(yīng)力將發(fā)生變化。根據(jù)旋轉(zhuǎn)變換矩陣可得到各個(gè)面上轉(zhuǎn)換后的應(yīng)力分量 σxx′、σyy′、σzz′、σxy′、σxz′和 σyz′。上述轉(zhuǎn)換是為了保證得到一個(gè)最大應(yīng)力σxx′和對(duì)應(yīng)的最小應(yīng)力 σyy′。

應(yīng)力轉(zhuǎn)換關(guān)系［8］表達(dá)式為

其中繞z軸旋轉(zhuǎn)的變化矩陣為

轉(zhuǎn)換坐標(biāo)后的應(yīng)力分量可表示為

式中，σxx、σyy、σzz和 σxx′、σyy′、σzz′分別為坐標(biāo)系轉(zhuǎn)換前、后的主應(yīng)力，MPa;σxy、σyz、σzx和 σxy′、σxz′、σyz′分別為坐標(biāo)轉(zhuǎn)換前、后的剪切應(yīng)力，MPa;γ為繞井眼軸線z逆時(shí)針旋轉(zhuǎn)的角度，rad。

圖1 坐標(biāo)系旋轉(zhuǎn)示意圖Fig.1 Sketch map of coordinate system transformation around z axis

在繞井眼軸線z逆時(shí)針旋轉(zhuǎn)的角度γ從0增加到π的過(guò)程中，會(huì)出現(xiàn)一個(gè)最大的主應(yīng)力σxx′和對(duì)應(yīng)的最小主應(yīng)力σyy′。研究表明，3個(gè)剪切應(yīng)力相對(duì)于主應(yīng)力值小一個(gè)數(shù)量級(jí)，計(jì)算時(shí)忽略剪切應(yīng)力的作用，即認(rèn)為巖石受力變形后其應(yīng)力分布仍保持在垂直井筒軸線的平面上，且井筒截面上應(yīng)力σxx和σyy的差值要比應(yīng)力值小一個(gè)數(shù)量級(jí)。所以，可認(rèn)為井筒受均勻壓應(yīng)力 σ(R)=(σxx′+σyy′)/2 的作用，建立數(shù)學(xué)模型可得到沿井周的破壞半徑。

1.2 擠壓條件下水平井近井破壞帶劃分

假設(shè):①水平井裸眼井段井眼軸線為直線，無(wú)井筒迂曲;②較短的井眼長(zhǎng)度單元情況下地應(yīng)力沿井軸方向無(wú)變化，只隨角度和半徑變化;③近井地帶巖石各向同性;④剪應(yīng)力數(shù)量級(jí)比正應(yīng)力小，可忽略;⑤近井地層圍巖受均勻壓應(yīng)力作用;⑥破壞模式為塑性擠壓充填。

如圖2所示，將近井地層分為井筒、中間塑性壓實(shí)區(qū)和遠(yuǎn)處彈性變形區(qū)。對(duì)水平井實(shí)施高壓擠壓充填作業(yè)時(shí)，水平井筒內(nèi)壓力為pd。開始充填時(shí)，攜砂液對(duì)井筒周圍巖石進(jìn)行擠壓充填，巖石發(fā)生變形。如圖2所示的充填壓實(shí)過(guò)程中水平井筒半徑由r0增大為r1，在塑性壓實(shí)區(qū)域半徑r2以外的地層巖石仍保持原始彈性變形狀態(tài)。

圖2 水平井近井?dāng)D壓充填帶劃分示意圖Fig.2 Sketch map of compaction packing zone divide around horizontal well bore

2 水平井管外地層塑性擠壓充填模型

2.1 彈性變形區(qū)的應(yīng)力與位移

采用柱坐標(biāo)系，對(duì)于給定的深度h，由r、θ確定空間上一點(diǎn)。任意點(diǎn)(r、θ)處的巖石應(yīng)變張量的幾何方程［8］為

式中，εz、εr和εθ分別為軸向應(yīng)變、徑向應(yīng)變和切向應(yīng)變;u為巖石單位徑向位移，m;w為巖石單位軸向位移，m。

根據(jù)廣義胡克定律可以得到巖石變形的本構(gòu)方程為

其中

式中，σr、σθ和σz分別為徑向、切向和垂向主應(yīng)力，MPa;υ為泊松比;E為彈性模量，MPa。

根據(jù)假設(shè)條件，軸向應(yīng)力變化和剪切應(yīng)力可以忽略，則垂直于井筒軸線的截面為平面應(yīng)變，變形后仍保持為平面。這樣巖石的應(yīng)力平衡方程式可簡(jiǎn)化為

聯(lián)立可得

采用分離變量法求解得到方程(7)的通解為

式中，c1和c2為積分常數(shù)。

對(duì)于彈性區(qū)，邊界條件為

式中，R為外邊界半徑，m;r2為塑性壓實(shí)區(qū)最大半徑，m。

求解可得

彈性區(qū)相對(duì)位移為

根據(jù)式(11)可得彈塑性交界處(r=r2)的相對(duì)位移為

2.2 塑性變形區(qū)域的應(yīng)力與位移

在水平井靠近井筒的塑性變形區(qū)域內(nèi)，地層巖石突破彈性變形條件，產(chǎn)生塑性破壞，其應(yīng)力－應(yīng)變關(guān)系遵從Mohr-Coulun準(zhǔn)則，可表示為

式中，σ1為最大主應(yīng)力，MPa;σ3為最小主應(yīng)力，MPa;C為巖石內(nèi)聚強(qiáng)度，MPa;φ為巖石內(nèi)摩擦角，rad。

在彈性變形階段，井筒內(nèi)邊界處的徑向及切向應(yīng)力表示為

對(duì)于水平井管外地層塑性擠壓充填，井內(nèi)壓力較高，一般有pd＞σ(R)。此時(shí) σθ＜σr，徑向應(yīng)力 σr表現(xiàn)為壓應(yīng)力，切向應(yīng)力σθ表現(xiàn)為拉應(yīng)力，水平井眼直徑擴(kuò)大，即有σ3=σθ和σ1=σr。根據(jù)方程(13)有

將式(15)帶入平衡方程(6)，得

在充填區(qū)與塑性區(qū)交界處，即r=r1處σr=pd，由此可得

將式(17)帶入式(16)得

在彈性區(qū)和塑性區(qū)交界處，即r=r2處有

由式(15)和邊界條件式(19)可得

化簡(jiǎn)可得到塑性區(qū)擴(kuò)張比為

式中，r1為裸眼塑性擴(kuò)張半徑，m。

2.3 管外地層擠壓充填區(qū)和壓實(shí)帶半徑計(jì)算

根據(jù)水平井塑性擠壓段的巖石體積變化等于彈性變形區(qū)的體積變化與塑性壓實(shí)帶的體積變化之和的原理，可得

式中，Δ為塑性壓實(shí)帶的平均體積應(yīng)變;r0為水平井筒半徑，m。

疏松砂巖油藏巖石由巖石骨架顆粒和孔隙空間組成，且顆粒間膠結(jié)疏松。在進(jìn)行高壓擠壓充填作業(yè)時(shí)，隨著泵壓增高，井底壓力增大，地層巖石介質(zhì)所承載的內(nèi)部應(yīng)力也不斷升高，巖石基質(zhì)顆粒間發(fā)生崩塌和聚結(jié)。擠壓充填過(guò)程中，水平井筒周圍地層巖石介質(zhì)顆粒之間的原始膠結(jié)被破壞，顆粒間隙變小，巖石孔隙度也變小。根據(jù)上述分析，塑性變形區(qū)的體積應(yīng)變可按疏松砂巖施工前后的巖石介質(zhì)孔隙度的比值計(jì)算:

式中，φ1為疏松砂巖地層巖石的原始孔隙度，即施工前的地層孔隙度;φ2為高壓擠壓塑性壓實(shí)區(qū)的巖石介質(zhì)孔隙度。

對(duì)方程(22)進(jìn)行簡(jiǎn)化并忽略［u(r2)］2項(xiàng)可得

在高壓擠壓作業(yè)中水平井井底施工壓力高于壓應(yīng)力，即pd＞σ(R)。因此根據(jù)計(jì)算得到的彈性變形區(qū)位移方程(21)可以推導(dǎo)得到

由式(25)可得井徑擴(kuò)張率為

方程(26)可用來(lái)計(jì)算塑性變形帶的擴(kuò)張率，從而進(jìn)一步得到井徑擴(kuò)張率r1/r0，最后可求出塑性擠壓充填帶的半徑r1和塑性壓實(shí)帶半徑r2。

3 模擬計(jì)算結(jié)果分析

某水平井井斜角88.63°，垂深1.18382 km，井斜方位角248.5°。最大原始水平主應(yīng)力方位角120°，利用測(cè)井計(jì)算得到的平均地層巖石密度2002.3 kg/m3，最大水平應(yīng)力構(gòu)造系數(shù)0.902，最小水平應(yīng)力構(gòu)造系數(shù)0.257。已知該井2.1～2.2 km井段的測(cè)井?dāng)?shù)據(jù)，并根據(jù)其計(jì)算得到巖石力學(xué)參數(shù)沿井身軌跡的分布規(guī)律。

3.1 管外擠壓充填破壞剖面的動(dòng)態(tài)模擬結(jié)果

利用本文模型對(duì)該算例水平井水平段管外塑性擠壓破壞剖面進(jìn)行動(dòng)態(tài)模擬。不同井底壓力下的充填形態(tài)如圖3所示。

根據(jù)圖3的模擬結(jié)果可知，由于水平段較長(zhǎng)，導(dǎo)致巖石力學(xué)參數(shù)及地應(yīng)力的非均質(zhì)性，管外地層塑性擠壓充填也表現(xiàn)出較強(qiáng)的非均勻性，其充填剖面的分布規(guī)律與巖石強(qiáng)度參數(shù)的分布規(guī)律相似。對(duì)比分析結(jié)果表明，巖石強(qiáng)度相對(duì)較低的井段，最終擠壓破壞的半徑相對(duì)較大。算例中，在該井后半段巖石強(qiáng)度相對(duì)較高，模擬計(jì)算得到的擠壓充填破壞半徑也相對(duì)較小。

圖3 水平井管外地層破壞剖面動(dòng)態(tài)模擬Fig.3 Dynamic simulation of failure radius distribution along horizontal wellbore

當(dāng)井底施工壓力較低時(shí)，水平段管外塑性擠壓破壞形態(tài)相對(duì)比較均勻，充填破壞半徑沿井身軌跡的變化差異不十分明顯。當(dāng)井底施工壓力達(dá)到較高值時(shí)，水平段管外塑性擠壓破壞形態(tài)及尺寸逐漸沿井身軌跡方向呈現(xiàn)出一定的非均勻性。

3.2 水平井?dāng)D壓充填橫向破壞形態(tài)靜態(tài)模擬結(jié)果

使用上述算例水平井井身2.15 km處的巖石力學(xué)數(shù)據(jù)進(jìn)行管外塑性擠壓充填破壞形態(tài)的模擬計(jì)算與分析。該位置處地層巖石的內(nèi)摩擦角為39.1°，巖石內(nèi)聚力為3.85 MPa，泊松比為0.314，巖石彈性模量為11952 MPa。不同井底壓力下的管外破壞形態(tài)如圖4所示。

由圖4可以看出:當(dāng)井底施工壓力較低時(shí)，水平井管外地層塑性擠壓破壞形態(tài)相對(duì)比較均勻，其沿井周角的非均質(zhì)性不十分明顯;而當(dāng)井底施工壓力逐步達(dá)到較高值后，管外地層巖石的塑性擠壓破壞形態(tài)也表現(xiàn)出一定的非均勻性。這表明，疏松砂巖油氣藏水平井在高壓擠注情況下，近井地帶地層巖石主要受井眼內(nèi)流體壓力的影響較大，而受外圍地層地應(yīng)力的影響相對(duì)較小。但隨著管外地層塑性擠壓破壞半徑的增加，外圍地層地應(yīng)力的作用和影響逐漸體現(xiàn)出來(lái)，表現(xiàn)為橫向塑性擠壓剖面的非均勻性。

圖4 水平井管外地層徑向破壞剖面動(dòng)態(tài)模擬Fig.4 Dynamic simulation of failure radius around horizontal wellbore

3.3 擠壓充填擴(kuò)張比敏感性分析

3.3.1 巖石內(nèi)摩擦角對(duì)充填結(jié)果的影響

設(shè)定如下計(jì)算參數(shù):內(nèi)聚力2.0 MPa，彈性模量800 MPa，巖石泊松比0.3，體積應(yīng)變0.1，平均井筒截面壓應(yīng)力σ(R)=20 MPa。不同疏松砂巖內(nèi)摩擦角條件下，高壓擠壓過(guò)程中，井徑擴(kuò)張率r1/r0及壓實(shí)區(qū)擴(kuò)張率r2/r1隨井底壓力pwf變化曲線如圖5所示。

在內(nèi)摩擦角相同的情況下，井徑擴(kuò)張率r1/r0隨井底壓力的增大而增大。在井底壓力略高于井筒截面壓應(yīng)力時(shí)井徑擴(kuò)張率r1/r0接近于1，說(shuō)明較低的井底壓力不會(huì)引起較大的井徑擴(kuò)張，但其增大速率隨著井底壓力的升高有加大的趨勢(shì)。在井底壓力相同的情況下，內(nèi)摩擦角減小會(huì)導(dǎo)致井徑擴(kuò)張率r1/r0增大。內(nèi)摩擦角數(shù)值越小，井徑擴(kuò)張率r1/r0增加的幅度越大;反之，增大幅度越小。

在內(nèi)摩擦角相同的情況下，井底壓力越高，塑性壓實(shí)區(qū)域擴(kuò)張率r2/r1越高，兩者成正比關(guān)系。當(dāng)井底壓力相同時(shí)，內(nèi)摩擦角越小，塑性壓實(shí)區(qū)擴(kuò)張率r2/r1越大。當(dāng)內(nèi)摩擦角數(shù)值較小時(shí)，塑性壓實(shí)區(qū)的擴(kuò)張率r2/r1增幅隨井底壓力的增大而增大。

3.3.2 內(nèi)聚力的影響

取疏松砂巖的內(nèi)摩擦角 20°，彈性模量 800 MPa，泊松比0.3，體積應(yīng)變?nèi)?.1，平均井筒截面壓應(yīng)力σ(R)=20 MPa。不同疏松砂巖內(nèi)聚力條件下，高壓擠壓過(guò)程中，井徑擴(kuò)張率r1/r0以及壓實(shí)區(qū)擴(kuò)張率r2/r1隨井底壓力的變化如圖6所示。

由圖6可以看出:在井底壓力較低的情況下，內(nèi)聚力對(duì)井徑擴(kuò)張率及壓實(shí)區(qū)擴(kuò)張率影響較小;在較高井底壓力下，隨內(nèi)聚力減小，內(nèi)聚力對(duì)井徑擴(kuò)張率及壓實(shí)區(qū)擴(kuò)張率的影響增大。當(dāng)內(nèi)聚力較低時(shí)，井底壓力增加到一定程度后，較小的壓力變化就可以引起較大的井徑擴(kuò)張率變化;當(dāng)內(nèi)聚力較高時(shí)，井徑擴(kuò)張率隨井底壓力的變化較小。

圖5 內(nèi)摩擦角對(duì)充填擴(kuò)張率的影響Fig.5 Influence of rock internal friction angel on extension ratio of packing

圖6 內(nèi)聚力對(duì)充填擴(kuò)張率的的影響Fig.6 Influence of cohesive strength on extension ratio of packing

3.3.3 彈性模量的影響

取疏松砂巖的內(nèi)摩擦角20°，內(nèi)聚力2 MPa，泊松比0.3，體積應(yīng)變0.1，平均井筒截面壓應(yīng)力 σ(R)=20 MPa。圖7為水平井高壓擠壓過(guò)程中不同彈性模量條件下的井徑擴(kuò)張率r1/r0隨井底施工壓力的變化關(guān)系曲線。

圖7 彈性模量對(duì)井徑擴(kuò)張率的影響Fig.7 Influence of elastic modulus on extension ratio of borehole radius

由圖7可看出，在壓力較低的情況下彈性模量對(duì)井徑擴(kuò)張率影響較小。隨井底壓力升高，彈性模量對(duì)井徑擴(kuò)張率影響變大。在較高的井底壓力下較小的壓力變化就可以引起較大的井徑擴(kuò)張率變化。

3.3.4 泊松比的影響

取疏松砂巖的內(nèi)摩擦角20°，內(nèi)聚力2 MPa，彈性模量800 MPa，體積應(yīng)變0.1，平均井筒截面壓應(yīng)力 σ(R)=20 MPa。

圖8 泊松比對(duì)井徑擴(kuò)張率的影響Fig.8 Influence of Poisson ratio on extension ratio of borehole radius

圖8為水平井高壓擠壓過(guò)程中不同泊松比下的井徑擴(kuò)張率r1/r0隨井底施工壓力的變化曲線。由圖8可看出，在低井底壓力情況下曲線基本重合，在較高的井底壓力下各曲線差異也很小，說(shuō)明泊松比對(duì)充填擴(kuò)張率的影響非常小，可忽略不計(jì)。

4 結(jié)論

(1)由于沿水平井井身軌跡巖石力學(xué)參數(shù)及地應(yīng)力的非均質(zhì)性，管外地層塑性擠壓充填剖面也表現(xiàn)出較強(qiáng)的非均勻性，并與巖石參數(shù)分布規(guī)律相關(guān)。巖石強(qiáng)度相對(duì)較低的井段，最終擠壓破壞的半徑相對(duì)較大。當(dāng)井底壓力較低時(shí)，水平段管外塑性擠壓破壞形態(tài)沿井身軌跡的變化不明顯;當(dāng)井底壓力較高時(shí)，擠壓破壞形態(tài)非均質(zhì)性逐步變得明顯。

(2)隨井底壓力增高，充填擴(kuò)張率對(duì)各影響因素的敏感性急劇增加。疏松砂巖巖石內(nèi)摩擦角、內(nèi)聚力、彈性模量隨著井底壓力的增加對(duì)充填擴(kuò)張率影響增大。泊松比的影響最小，可以忽略。

［1］ 董長(zhǎng)銀.油氣井防砂技術(shù)［M］.北京:中國(guó)石化出版社，2009:48-60.

［2］ 董長(zhǎng)銀，張琪.水平井礫石充填過(guò)程實(shí)時(shí)數(shù)值模擬研究［J］.石油學(xué)報(bào)，2004，25(6):96-100.

DONG Chang-yin，ZHANG Qi.Real-time numerical simulation of gravel-packing process in horizontal wells［J］.Acta Petrolei Sinica，2004，25(6):96-100.

［3］ 李明忠，李彥超，王衛(wèi)陽(yáng)，等.考慮井筒變質(zhì)量流動(dòng)的礫石充填水平井產(chǎn)能預(yù)測(cè)［J］.中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版，2011，35(3):89-93.

LI Ming-zhong，LI Yan-chao，WANG Wei-yang，et al.Productivity prediction for gravel-packed horizontal well sonsidering variable mass flow in wellbore［J］.Journal of China University of Petroleum(Edition of Natural Science)，2011，35(3):89-93.

［4］ 董長(zhǎng)銀，劉春苗，武龍，等.高壓擠注條件下水平井管外地層破壞模式及其判別方法［J］.中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版，2012，36(5):45-49.

DONG Chang-yin，LIU Chun-miao，WU Long，et al.Formation rock failure mode and its recognition method under high pressure squeezing condition in horizontal wells［J］.Journal of China University of Petroleum(Edition of Natural Science)，2012，36(5):45-49.

［5］ 陳勉，陳治喜，黃榮樽.大斜度井水壓裂縫起裂研究［J］.石油大學(xué)學(xué)報(bào):自然科學(xué)版，1995，19(2):30-35.

CHEN Mian，CHEN Zhi-xi，HUANG Rong-zun.Hydraulic fracturing of highly deviated wells［J］.Journal of the University of Petroleum，China(Edition of Natural Science)，1995，19(2):30-35.

［6］ 閆相禎，王偉章，楊秀娟，等.近井地帶高壓擠壓?jiǎn)栴}的解析［J］.中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版，2008，32(3):103-107.

YAN Xiang-zhen，WANG Wei-zhang，YANG Xiu-juan，et al.Analytical solution of high-pressure squeezing problems at stratum near well［J］.Journal of China University of Petroleum(Edition of Natural Science)，2008，32(3):103-107.

［7］ 王偉章，閆相禎.高壓擠壓防砂充填帶形態(tài)仿真模擬［J］.濰坊學(xué)院學(xué)報(bào)，2009，9(2):1-5.

WANG Wei-zhang，YAN Xiang-zhen.Numerical simulation of packing zone morphology with high pressure gravelpacking sand-prevention technology in unconsolidated sandstone［J］.Journal of Weifang University，2009，9(2):1-5.

［8］ 閆鐵，李士斌.深部井眼巖石力學(xué)理論與實(shí)踐［M］.北京:石油工業(yè)出版社，2002:11-33.