張潔，陳文，康曉娜，劉松青，趙志良，王雙威

（1.中國(guó)石油鉆井工程技術(shù)研究院，北京100083；2.北京師范大學(xué)，北京100875；3.中國(guó)石油長(zhǎng)城鉆探工程公司井下作業(yè)公司，北京100101）

儲(chǔ)層段漏失易造成儲(chǔ)層污染嚴(yán)重等諸多問(wèn)題［1-3］，發(fā)展儲(chǔ)層段防漏堵漏技術(shù)迫在眉睫，但一直沒(méi)有配套的評(píng)價(jià)方法。采用常規(guī)堵漏儀，很難與儲(chǔ)層參數(shù)（滲透率、孔隙度、孔隙直徑等）相匹配［4-7］；而采用常規(guī)的儲(chǔ)層保護(hù)評(píng)價(jià)方法，不僅很難真實(shí)地反映儲(chǔ)層段堵漏材料的堵漏和儲(chǔ)層保護(hù)效果，而且會(huì)損害儀器。國(guó)外學(xué)者發(fā)明的可視式中壓砂床濾失儀，可以直觀定量地評(píng)價(jià)儲(chǔ)層段滲漏及中小漏的堵漏效果。但是，如何根據(jù)儲(chǔ)層孔隙度、孔隙直徑等參數(shù)確定砂床規(guī)格，一直沒(méi)有規(guī)范可以參考。為此，筆者借助掃描電鏡，利用回歸分析模型研究了玻璃微珠目數(shù)與砂床特征的關(guān)系。

1 實(shí)驗(yàn)方法

1.1 孔隙直徑測(cè)定

將6 種不同規(guī)格的玻璃微珠烘干24 h 以上，分別放進(jìn)玻璃圓筒中，輕拍上口多次，使得玻璃微珠緊實(shí)穩(wěn)定。然后，迅速將圓筒倒置于導(dǎo)電膠布之上，并施加一定的壓力，使得玻璃微珠不至于變形而影響孔隙直徑的測(cè)定。最后，將6 種樣品分別放入掃描電鏡樣品池中，通過(guò)掃描電鏡自帶的圖片分析軟件測(cè)定孔隙直徑。

1.2 孔隙體積及孔隙度測(cè)定

將玻璃微珠烘干至少24 h，分別將50 mL 玻璃微珠裝入內(nèi)徑為2.5 cm 的具塞量筒中，緩緩加入蒸餾水，使得砂床液面保持在50 mL。速度要慢，時(shí)間要充足，使得水充分滲透，將孔隙中的空氣排干凈。記錄此時(shí)滴入具塞量筒中的水的體積，即為孔隙體積Vp，而砂床的體積為50 mL，則孔隙度為φ=Vp/50×100%。

2 實(shí)驗(yàn)結(jié)果與數(shù)據(jù)處理

2.1 孔隙直徑測(cè)定結(jié)果

整理不同目數(shù)的玻璃微珠對(duì)應(yīng)砂床的孔隙直徑數(shù)值，結(jié)果如表1所示。

表1 不同目數(shù)玻璃微珠對(duì)應(yīng)砂床的孔隙直徑數(shù)值

2.2 回歸分析模型對(duì)比

回歸分析是研究因變量與自變量之間變動(dòng)比例關(guān)系的一種方法，分為線性回歸方法和非線性回歸方法2 種［8-10］。本文以線性回歸模型、線性模型EM 算法及非線性模型EM 算法分別建立了相應(yīng)的3 個(gè)回歸分析模型。

2.2.1 線性回歸模型

線性回歸模型為Y=β0+β1X+e，e～N（0，σ2）。其中，X 為玻璃微珠的目數(shù)，Y 為砂床的孔隙直徑，e 為殘差，σ 為標(biāo)準(zhǔn)估計(jì)誤差，N（0，σ2）為正態(tài)分布。應(yīng)用最小二乘估計(jì)法，得到的估計(jì)值是＝153.276，＝-1.309，＝33.523 10。線性回歸模型的優(yōu)點(diǎn)是簡(jiǎn)單方便，但只能用到確定目數(shù)的數(shù)據(jù)，20～40 目及40～60 目玻璃微珠對(duì)應(yīng)砂床的孔隙直徑無(wú)法用上，此模型不一定合理。

2.2.2 線性模型EM 算法

該方法是保留確定目數(shù)及范圍目數(shù)的數(shù)據(jù)。線性EM 算法模型為Y=β0+β1X+e，e～N（0，σ2）。應(yīng)用最小二乘估計(jì)法，得到的估計(jì)值是＝178.827 0，＝－1.712 352，＝33.476 21。相較于線性回歸模型，線性模型EM 算法的優(yōu)點(diǎn)是用到了所有的數(shù)據(jù)。

2.2.3 非線性模型EM 算法

非線性EM 算法模型為Y=β0+β1X-1+e，e～N（0，σ2）。應(yīng)用最小二乘估計(jì)法，得到的估計(jì)值是＝17.263 34，＝3 132.142，＝29.306 08。非線性模型EM算法與線性回歸模型對(duì)比，優(yōu)點(diǎn)是用到了所有數(shù)據(jù)；與線性模型EM 算法對(duì)比，優(yōu)點(diǎn)是考慮了實(shí)際情況，即砂床孔隙直徑與玻璃微珠目數(shù)的關(guān)系呈反比例趨勢(shì)。

這3 種模型的顯著性均小于0.05，說(shuō)明模型都通過(guò)了參數(shù)顯著性檢驗(yàn)，建立的回歸方程具有統(tǒng)計(jì)學(xué)意義。3 種模型的擬合程度對(duì)比結(jié)果，如表2所示。其中，R2表示模型的擬合程度，其值越大，說(shuō)明模型擬合越好。由表2可看出，非線性模型EM 算法的R2最大，說(shuō)明其擬合程度相對(duì)其他2 種方法更好。所以，決定采用非線性模型EM 算法，來(lái)描述玻璃微珠目數(shù)與緊密排列時(shí)孔隙直徑之間的關(guān)系。

表2 模型擬合程度對(duì)比

2.3 模型預(yù)測(cè)關(guān)系

非線性模型EM 算法的結(jié)果為

Y=17.263 34+3 132.142/X+e，e～N（0，29.306 082）

用橫軸表示殘差e 的假設(shè)分布的分位數(shù)，縱軸表示真實(shí)分布的分位數(shù)，繪制該模型的常規(guī)Q-Q 圖，表征殘差e 的真實(shí)分布（真實(shí)模型）和假設(shè)分布（假設(shè)模型）的匹配程度（見(jiàn)圖1）。

圖1 常規(guī)Q-Q

從圖1可看出，這些散點(diǎn)恰好落在第一象限的角分線上，說(shuō)明真實(shí)分布與假設(shè)分布對(duì)應(yīng)良好，即殘差e的真實(shí)分布可認(rèn)為呈假設(shè)的正態(tài)分布。以殘差e 為縱坐標(biāo)、以Y 值為橫坐標(biāo)，繪制該模型殘差圖，以驗(yàn)證殘差e 的真實(shí)分布是否為假設(shè)的正態(tài)分布（見(jiàn)圖2）。

從圖2可看出，殘差e=Y-17.263 34-3 132.142X-1，這些點(diǎn)均勻散落在以0 軸為中心的帶狀區(qū)域內(nèi)，說(shuō)明殘差e 的分布是均值為0 的正態(tài)分布。將20～40 目及40～60 目玻璃微珠的真實(shí)目數(shù)補(bǔ)出來(lái)，繪制Y 與X-1的關(guān)系圖（見(jiàn)圖3）。結(jié)果表明，用模型預(yù)測(cè)補(bǔ)出的散點(diǎn)與真實(shí)的情況大概落在一條直線上，說(shuō)明建立Y 與X-1為直線關(guān)系是合理的。

圖2 非線性EM 模型殘差

圖3 Y-X-1 關(guān)系

下面根據(jù)非線性模型EM 算法的結(jié)果，來(lái)預(yù)測(cè)砂床孔隙直徑最大概率下的波動(dòng)范圍。當(dāng)X 取20 目時(shí)，孔隙半徑應(yīng)在173.870 μm 附近波動(dòng)，在置信度為95%的情況下，波動(dòng)范圍為116.431 60～231.309 30 μm。當(dāng)X 取其他值時(shí)，用上述方法進(jìn)行估計(jì)，結(jié)果如表3所示。

表3 玻璃微珠目數(shù)與砂床孔隙直徑的關(guān)系

2.4 孔隙體積及孔隙度測(cè)定結(jié)果

不同目數(shù)玻璃微珠對(duì)應(yīng)砂床的孔隙體積及孔隙度值，如表4所示。

表4 玻璃微珠對(duì)應(yīng)砂床的孔隙體積及孔隙度

根據(jù)表3及表4中不同目數(shù)玻璃微珠對(duì)應(yīng)砂床的孔隙直徑及孔隙度的范圍，結(jié)合儲(chǔ)層段巖心的壓汞數(shù)據(jù)分析及常規(guī)孔滲分析數(shù)據(jù)，選擇最匹配目數(shù)的玻璃微珠作為砂床。

3 結(jié)論

1）測(cè)定了玻璃微珠對(duì)應(yīng)砂床的孔隙直徑及孔隙度值，并通過(guò)3 種模型對(duì)比，采用擬合度最高的非線性模型EM 算法來(lái)建立玻璃微珠目數(shù)與緊密排列時(shí)孔隙直徑的對(duì)應(yīng)關(guān)系，并根據(jù)該模型預(yù)測(cè)了95%置信度下砂床孔隙直徑的波動(dòng)范圍。

2）通過(guò)預(yù)測(cè)不同目數(shù)玻璃微珠對(duì)應(yīng)砂床的孔隙直徑及孔隙度的范圍，建立了玻璃微珠目數(shù)與砂床特征的對(duì)應(yīng)關(guān)系，從而改進(jìn)了可視式中壓砂床濾失評(píng)價(jià)方法，建立了一種可以與儲(chǔ)層參數(shù)相匹配的儲(chǔ)層段防漏堵漏評(píng)價(jià)方法。

［1］林英松，蔣金寶，秦濤.井漏處理技術(shù)的研究及發(fā)展［J］.斷塊油氣田，2005，12（2）：4-8.

［2］王業(yè)眾，康毅力，游利軍，等.裂縫性?xún)?chǔ)層漏失機(jī)理及控制技術(shù)進(jìn)展［J］.鉆井液與完井液，2007，24（4）：74-77.

［3］閆豐明，康毅力，孫凱，等.縫洞型碳酸鹽巖儲(chǔ)層漏失模型及控制對(duì)策［J］.鉆井液與完井液，2012，29（3）：78-80.

［4］竇斌，舒尚文，郭建華，等.高保真模擬漏失地層堵漏評(píng)價(jià)試驗(yàn)裝置設(shè)計(jì)［J］.石油機(jī)械，2009，37（11）：5-7.

［5］楊振杰，張玉強(qiáng)，吳偉，等.XAN-D 新型堵漏評(píng)價(jià)試驗(yàn)儀研究與應(yīng)用［J］.鉆井液與完井液，2010，27（1）：44-46.

［6］余維初，蘇長(zhǎng)明，鄢捷年，等.高溫高壓動(dòng)態(tài)堵漏評(píng)價(jià)系統(tǒng)［J］.鉆井液與完井液，2009，26（1）：20-22.

［7］余維初，蘇長(zhǎng)明，黃新然，等.高溫高壓動(dòng)態(tài)堵漏評(píng)價(jià)系統(tǒng)的研制［J］.石油天然氣學(xué)報(bào)，2008，30（5）：116-117.

［8］李華，劉帥，李茂，等.數(shù)據(jù)挖掘理論及應(yīng)用研究［J］.斷塊油氣田，2010，17（1）：88-91.

［9］杭義萍，田淑云，游革新.線性回歸方法處理熱釋汞分析數(shù)據(jù)［J］.石油實(shí)驗(yàn)地質(zhì)，1995，17（1）：96-98.

［10］胡慧芳.油藏?cái)?shù)值模擬歷史擬合質(zhì)量評(píng)價(jià)方法［J］.斷塊油氣田，2012，19（3）：354-358.