劉飛，趙立強（油氣藏地質(zhì)及開發(fā)工程國家重點實驗室（西南石油大學(xué)），四川 成都 610500）

孟向麗（中海油能源發(fā)展股份有限公司工程技術(shù)分公司，天津 300452）

馬旭，趙振峰（中石油長慶油田分公司油氣工藝研究院低滲透油氣田國家工程實驗室，陜西 西安 710021）

復(fù)雜巖性、非常規(guī)、特殊油氣藏的勘探開發(fā)對常規(guī)增產(chǎn)技術(shù)提出了巨大挑戰(zhàn)，常規(guī)水力壓裂由于壓裂液破膠不徹底、壓裂液殘渣滯留、壓裂液濾液傷害、壓裂液濾餅傷害等不能充分發(fā)揮水力裂縫的增產(chǎn)改造效果①趙立強，等 .前置酸酸液體系優(yōu)化 .西南石油大學(xué)，2011.［1，2］，將酸化和壓裂工藝相結(jié)合是未來增產(chǎn)改造技術(shù)的重要發(fā)展方向［3］，前置酸加砂壓裂工藝［4～8］在這方面做了主要嘗試，并已在現(xiàn)場取得了較好的應(yīng)用效果。前期研究成果表明：前置酸通過溶蝕或穩(wěn)定濾失區(qū)黏土礦物，改善近縫基巖孔隙度和滲透率，減少后續(xù)壓裂液濾液造成的水鎖傷害、返排期間促使壓裂液破膠、溶解壓裂液殘渣、清洗支撐裂縫等，可在常規(guī)水力壓裂的基礎(chǔ)上增產(chǎn)10%～23%；而前置酸通常采用有機土酸、多氫酸、氟硼酸、復(fù)合酸等緩速酸液體系，確保壓裂施工結(jié)束后濾失區(qū)酸液仍具有一定反應(yīng)活性。目前對砂巖與緩速酸液體系的酸巖反應(yīng)模擬理論尚不完善，前置酸加砂壓裂工藝設(shè)計帶有極強的主觀經(jīng)驗性。筆者基于緩速酸液體系緩速機理、緩速酸液與砂巖礦物的反應(yīng)機理，結(jié)合前置酸加砂壓裂工藝，建立水力裂縫內(nèi)酸濃度分布模型和濾失區(qū)酸巖反應(yīng)的 “三酸三礦物”模型。以HCl＋HBF4＋HF的復(fù)合緩速酸液體系為例，通過模擬計算裂縫內(nèi)酸液濃度分布，作為計算濾失區(qū)酸液濃度的內(nèi)邊界條件，模擬濾失區(qū)酸巖反應(yīng)過程可確定出濾失區(qū)酸濃度分布、礦物濃度分布以及基巖孔隙度和滲透率的改善程度。本研究完善了前置酸加砂壓裂理論，模擬計算結(jié)果為優(yōu)選緩速酸液類型、優(yōu)化前置酸用量提供數(shù)據(jù)支撐，并能指導(dǎo)前置酸加砂壓裂工藝設(shè)計。

1 裂縫內(nèi)酸巖反應(yīng)模擬

在前置酸壓裂中，采用緩速酸液體系壓開儲層，形成并延伸裂縫，酸液在裂縫內(nèi)的濃度分布是計算酸蝕裂縫導(dǎo)流能力的基礎(chǔ)，同時也是計算酸液濾失區(qū)酸液濃度分布的內(nèi)邊界條件。

假設(shè)條件：①酸液不可壓縮；②HF／H2SiF6與縫壁巖石的酸巖反應(yīng)受表面反應(yīng)控制；③碳酸鹽巖礦物優(yōu)先與注入酸液中的HCl反應(yīng)；④裂縫內(nèi)酸巖反應(yīng)符合一次反應(yīng)動力學(xué)。

取一個裂縫微元段，如圖1所示。根據(jù)酸液質(zhì)量守恒可以建立酸液在裂縫內(nèi)的濃度分布偏微分方程組。

圖1 裂縫微元示意圖

式中：A 為距井筒x 處的裂縫截面積，m2；C0為緩速酸濃度，mol／m3；C1為 HF濃度，mol／m3；C2為 H2SiF6濃度，mol／m3；u為距井筒x處裂縫內(nèi)的酸液流速，m／s；R0為緩速酸反應(yīng)速度，mol／（m3·s）；R1為HF與巖石的表面反應(yīng)速度，mol／（m3·s）；R2為H2SiF6與巖石的表面反應(yīng)速度，mol／（m3·s）；H 為裂縫高度，m；ula為酸液濾失速度，m／s；λ1為1mol緩速酸反應(yīng)生成 HF的物質(zhì)量，mol／mol；λ2為1mol HF反應(yīng)生成 H2SiF6的物質(zhì)量，mol／mol。

由于酸液不可壓縮，裂縫內(nèi)流體的質(zhì)量守恒可以轉(zhuǎn)化為體積守恒，則裂縫內(nèi)流體的連續(xù)性方程為：

緩速酸反應(yīng)生成HF的反應(yīng)速度R0可表示為［9～13］：

式中：k0為反應(yīng)速度常數(shù)，（mol／m3）－1·s－1；CH＋ 為酸液中H＋濃度，mol／m3；Ea為水解活化能，J／mol；R為通用氣體常數(shù)，8.314J／（mol·K）；T 為溫度，K。

由于酸液與縫壁巖石的酸巖反應(yīng)受表面反應(yīng)控制，故R1和R2可以表示為：

式中：Ef，i，j為酸液i與礦物j的反應(yīng)速度常數(shù)，mol·（m2·s·mol／m3）；Ea，i，j為酸液i與礦物j的反應(yīng)活化能，J／mol；n為與酸液i發(fā)生反應(yīng)的礦物數(shù)；為礦物j的比表面積，m2／m3；Vj為礦物j的體積分數(shù)，m3／m3；Ci，eqm為酸液i與礦物反應(yīng)的平衡濃度（常取0），mol／m3。

2 濾失區(qū)酸巖反應(yīng)模擬

在酸液濾失侵入?yún)^(qū)，假設(shè)酸液流動符合一維線性單相流動。緩速酸濃度分布可根據(jù)水解反應(yīng)動力學(xué)得到；HF和H2SiF6的濃度分布可以采用兩酸三礦物模型模擬［14］，其區(qū)別在于HF的來源除了濾液本身之外，還存在緩速酸液水解反應(yīng)生成的HF。采用酸液質(zhì)量守恒可以建立相應(yīng)的 “三酸三礦物模型”：

3 計算實例分析

在前置酸加砂壓裂中，裂縫內(nèi)酸液濃度分布是計算濾失區(qū)內(nèi)酸液濃度分布的內(nèi)邊界條件，而裂縫內(nèi)酸液濃度分布又需要裂縫幾何尺寸作為基礎(chǔ)參數(shù)，同時，裂縫內(nèi)及濾失區(qū)溫度場又決定了緩速酸水解反應(yīng)和酸巖反應(yīng)速度的快慢［15，16］。因此，模型的求解需要耦合多個方程組。具體求解步驟如下：①計算井筒溫度場；②計算某一注酸時刻形成裂縫的幾何尺寸參數(shù)；③計算裂縫溫度場；④計算裂縫內(nèi)酸液濃度分布；⑤計算濾失區(qū)內(nèi)的酸液及礦物濃度分布；⑥計算濾失區(qū)內(nèi)由于酸巖反應(yīng)引起的孔隙度和滲透率變化。

儲層及酸巖反應(yīng)基礎(chǔ)參數(shù)見表1，緩速酸液體系采用10%HCl＋10%HBF4＋1%HF，模擬計算結(jié)果如圖2～6所示，其中HBF4和H2SiF6的無量綱濃度為其實際濃度與HBF4的注入濃度的比值，而HF的無量綱濃度為其實際濃度與HF注入濃度之比。

表1 儲層及酸巖反應(yīng)基礎(chǔ)參數(shù)

裂縫內(nèi)酸液濃度分布如圖2所示。從圖2可以看出：裂縫內(nèi)HBF4無量綱濃度因水解反應(yīng)生成HF而緩慢降低；而HF無量綱濃度先升高后緩慢降低，雖然縫口段HBF4無量綱濃度高但溫度較低，水解反應(yīng)相對較慢，HF無量綱濃度緩慢增加，裂縫中前段的HBF4濃度和流體溫度均較高，故HF濃度增加速度較快，裂縫尖端的溫度較高但HBF4無量綱濃度低，HF與巖石的反應(yīng)速度大于HBF4水解生成HF的速度，故HF無量綱濃度開始緩慢降低。

圖2 不同注酸量時裂縫內(nèi)酸液濃度分布

濾失區(qū)內(nèi)酸液濃度分布如圖3所示?？傮w而言，隨濾失深度的增加，HBF4和HF的無量綱濃度整體呈現(xiàn)降低趨勢，HBF4僅在縫口較小的濾失區(qū)內(nèi)維持較高的無量綱濃度，并隨濾失深度增加快速降低；對于HF而言，由于存在HBF4水解作用，HF可維持較高無量綱濃度，整個濾失區(qū)HF無量綱濃度維持在0.5以上，縫口段濾失區(qū)HF濃度幾乎維持為注入濃度，而裂縫中前段濾失區(qū)濃度均大于HF的注入濃度；H2SiF6僅在裂縫尖端濾失區(qū)出現(xiàn)較高濃度，無量綱濃度可達到0.2左右，能溶解一定量的快反應(yīng)礦物，對改善基巖孔隙度和滲透率具有一定積極作用。

圖3 前置酸150m3時濾失區(qū)酸液無量綱濃度

圖4 前置酸150m3時濾失區(qū)礦物體積分數(shù)

濾失區(qū)礦物體積分數(shù)分布如圖4所示。隨濾失深度增加，快反應(yīng)礦物和慢反應(yīng)礦物體積分數(shù)先減小后增大，直至等于原始體積分數(shù)，且快反應(yīng)礦物變化幅度遠大于慢反應(yīng)礦物。說明酸液主要溶蝕黏土類快反應(yīng)礦物，對石英等慢反應(yīng)礦物骨架溶蝕能力較??；硅膠沉淀僅出現(xiàn)在裂縫中前段的濾失區(qū)內(nèi)，即HF高濃度區(qū)域，說明HF濃度增加會導(dǎo)致二次沉淀的急劇增大，但總體濃度極小，體積分數(shù)數(shù)量級為10－4，產(chǎn)生的二次沉淀傷害微不足道。

施工結(jié)束時濾失區(qū)孔隙度和滲透率分布如圖5、6所示?？紫抖群蜐B透率均隨濾失深度增加而減小，最終降低至原始值，近裂縫壁面濾失區(qū)孔隙度從初始15%增加至32%，而滲透率從初始4mD增加至40mD，濾失區(qū)孔滲物性得到顯著改善。

圖5 施工結(jié)束時濾失區(qū)孔隙度分布

圖6 施工結(jié)束時濾失區(qū)滲透率分布

4 結(jié)論

1）用HBF4與HF的復(fù)合緩速酸液體系進行前置酸加砂壓裂時，裂縫內(nèi)HF濃度先增加后降低，酸液緩速性能較好。

2）緩速酸液體系可維持酸液濾失區(qū)在施工結(jié)束后保持相對較高的HF濃度，返排階段能促使壓裂液破膠、溶解壓裂液殘渣和清洗支撐裂縫，改善常規(guī)水力壓裂效果。

3）濾失區(qū)酸巖反應(yīng)可顯著改善基巖孔隙度和滲透率，對降低油氣流體向支撐裂縫匯聚的滲流阻力、提高壓裂井產(chǎn)能具有積極作用。

［1］Liu Pingli，Liu Fei，Zhao Liqiang，et al.Predicting productivity of hydraulic fracturing with pre－pad acid ［A］.ARMA 13－500Presented at 47thUS Rock Mechanics／Geomechanics Symposium ［C］.San Francisco，2013.

［2］劉飛 .前置酸壓裂工藝模擬研究 ［D］.成都：西南石油大學(xué)，2012.

［3］Kalfayan L J.Fracture acidizing：history，present state and future ［J］.SPE106371－MS，2007.

［4］劉紅現(xiàn)，趙立強，劉平禮，等 .前置酸改善水力壓裂效果機理分析 ［J］.石油與天然氣學(xué)報（江漢石油學(xué)院學(xué)報），2008，30（1）：332～334.

［5］王國峰，張立貴，李邊生，等 .前置酸壓裂工藝技術(shù)研究與應(yīng)用 ［J］.石油鉆采工藝.1997，19（增）：102～104.

［6］羅云，劉愛華，王俊明，等 .交聯(lián)酸加砂酸化壓裂技術(shù)在復(fù)雜巖性油藏的應(yīng)用 ［J］.石油學(xué)報，2008，29（2）：266～269.

［7］王文東，趙立強，劉平禮，等 .酸化加砂壓裂復(fù)合增產(chǎn)技術(shù) ［J］.斷塊油氣田，2009，16（5）：113～115.

［8］Rozo R，Paez J，Mendoza A，et al.Combining acid and hydraulic－fracturing technologies is the key to successfully stimulating the Orito Formation ［J］.SPE104610－MS，2007.

［9］王寶峰，趙忠揚，薛芳渝 .砂巖基質(zhì)酸化中的化學(xué)平衡研究 ［J］.西南石油學(xué)院學(xué)報，2000，22（2）：57～61.

［10］王寶峰，趙忠揚，薛芳渝 .砂巖基質(zhì)酸化中HF與硅鋁酸鹽的二次、三次反應(yīng)研究 ［J］.西南石油學(xué)院學(xué)報，2002，24（5）：61～64.

［11］李年銀，趙立強，劉平禮，等 .多氫酸酸化技術(shù)及其應(yīng)用 ［J］.西南石油大學(xué)學(xué)報，2009，31（6）：131～134.

［12］邢希金，唐洪明，趙峰，等 .伊利石與土酸／氟硼酸反應(yīng)實驗研究 ［J］.西南石油大學(xué)學(xué)報，2007，29（6）：29～31.

［13］曾凡輝，郭建春，趙金洲 .酸損傷降低砂巖儲層破裂壓力實驗研究 ［J］.西南石油大學(xué)學(xué)報，2009，31（6）：93～96.

［14］Li Chunlou.Fine scale sandstone acidizing core flood simulation ［D］ .Austin：University of Texas，2004.

［15］馬新仿，黃少云，張士誠，等.全三維水力壓裂過程中裂縫及近縫地層的溫度計算模型 ［J］.石油大學(xué)學(xué)報，2001，25（5）：38～41.

［16］王鴻勛，張士誠 .水力壓裂設(shè)計數(shù)值計算方法 ［M］.北京：石油工業(yè)出版社，1998.