蔣星達(dá)， 張 偉， 王仔軒， 李 光， 陳曉非

(1.中國(guó)科學(xué)技術(shù)大學(xué) 地球和空間科學(xué)學(xué)院，合肥 230026；2.南方科技大學(xué) 地球與空間科學(xué)系，深圳 518055；3.遼河油田勘探開(kāi)發(fā)研究院，盤(pán)錦 124000；4.東方地球物理公司 遼河物探處開(kāi)發(fā)地震服務(wù)中心，盤(pán)錦 124000)

0 引言

在水力壓裂條件下，微地震監(jiān)測(cè)裂縫成像精度受檢波器的布設(shè)位置、拾取誤差、速度模型準(zhǔn)確度、定位算法以及解釋算法等因素的影響，其中速度模型的準(zhǔn)確度對(duì)裂縫的定位結(jié)果具有重大影響[1]。在井下監(jiān)測(cè)過(guò)程中，聲波測(cè)井?dāng)?shù)據(jù)經(jīng)常被用于構(gòu)建井中監(jiān)測(cè)初始速度模型[2]。在地下地質(zhì)條件較為簡(jiǎn)單的地區(qū)，可以利用垂向一維層狀速度模型近似研究區(qū)域速度模型作為初始速度模型，這種簡(jiǎn)化的速度模型在很多地區(qū)得到了成功地應(yīng)用[3-4]。然而，獲取的聲波測(cè)井資料頻率范圍為幾百赫茲到幾千赫茲，而微地震信號(hào)相對(duì)頻率較低，且獲取的聲波測(cè)井信號(hào)傳播方向主要沿井筒方向，不能反映地層橫向速度變化。因此，獲取的初始速度模型通常需要根據(jù)射孔點(diǎn)信號(hào)進(jìn)行校正。

針對(duì)速度模型的校正，許多學(xué)者根據(jù)初始速度模型確定層位信息，利用各種算法校正各層速度值的大小。Oye等[5]提出通過(guò)擬合P波或者S波的走時(shí)信息獲得對(duì)應(yīng)的速度模型和震源位置;Pei等[4]提出同時(shí)反演速度模型和射孔信號(hào)發(fā)震時(shí)刻的方法，通過(guò)反演起始時(shí)刻獲得了走時(shí)信息;Zhou[6]提出以某個(gè)地震臺(tái)站記錄到的信號(hào)作為參考信號(hào)，計(jì)算其他臺(tái)站信號(hào)與該信號(hào)記錄的到時(shí)差作為觀測(cè)數(shù)據(jù)，計(jì)算地震事件位置；Zhou等[7]利用同一地震事件不同檢波器的P波到時(shí)差和同一檢波器的P-S波到時(shí)差作為觀測(cè)數(shù)據(jù)，在目標(biāo)函數(shù)中去掉了發(fā)震時(shí)刻的影響。利用到時(shí)差作為觀測(cè)數(shù)據(jù)，可以有效解決起始時(shí)刻不準(zhǔn)確帶來(lái)的問(wèn)題。

目前，微地震速度模型校正的許多方法都是對(duì)速度模型的層位深度或者速度值大小進(jìn)行校正，這些方法沒(méi)有考慮到聲波建立的初始速度模型層位個(gè)數(shù)、層位深度以及速度值大小均存在誤差的情況。為了避免由于聲波測(cè)井曲線的誤差而導(dǎo)致的速度模型誤差問(wèn)題，筆者提出一種不依賴于聲波測(cè)井曲線、僅使用射孔信號(hào)數(shù)據(jù)校正速度模型的方法。

針對(duì)無(wú)法準(zhǔn)確知道射孔信號(hào)的初始時(shí)間問(wèn)題，我們利用不同檢波器的P波走時(shí)差、S波走時(shí)差和同一檢波器的P/S波走時(shí)差作為觀測(cè)數(shù)據(jù)。以均勻但是具有細(xì)層間隔的速度模型作為初始速度模型，避免了采用聲波測(cè)井?dāng)?shù)據(jù)進(jìn)行層位劃分。除此之外，引入總變分正則化約束，對(duì)一維速度模型進(jìn)行反演校正?？傋兎?TV)正則化約束方法廣泛地應(yīng)用于圖像的去噪領(lǐng)域[8]，在去除噪聲的基礎(chǔ)上能夠保留圖像的邊緣信息和不連續(xù)性，達(dá)到保留圖像紋理細(xì)節(jié)的目的。地下介質(zhì)具有明顯的成層性，在地球物理反演過(guò)程中利用TV正則化約束，相對(duì)于Tikhonov約束，可以更好地保護(hù)介質(zhì)的邊緣信息。井下微地震監(jiān)測(cè)中的一維速度模型具有明顯的成層性，利用TV正則化約束速度模型，在擬合射孔信號(hào)走時(shí)差的基礎(chǔ)之上保護(hù)了層位邊界，最終獲得的速度模型具有不連續(xù)的速度值大小。這樣根據(jù)速度值的不同可以分辨層位，同時(shí)確定了層位的個(gè)數(shù)、層位深度和各層速度值的大小。

1 方法原理

(1)

針對(duì)同一個(gè)射孔信號(hào)，第j個(gè)檢波器和第i個(gè)檢波器的P波到時(shí)差表示為式(2)。

(2)

同理，可以計(jì)算出第j個(gè)檢波器和第i個(gè)檢波器的S波到時(shí)差為式(3)。

(3)

第j個(gè)檢波器P波信號(hào)和S波信號(hào)的到時(shí)差式(4)。

(4)

根據(jù)式(4)可以明顯看出，檢波器的到時(shí)差即為走時(shí)差，與射孔信號(hào)的發(fā)震時(shí)刻沒(méi)有關(guān)系。利用以上得到時(shí)差信息，可以構(gòu)建如下觀測(cè)函數(shù)：

(5)

一維層狀速度模型在速度值上表現(xiàn)為“階梯性”，利用具有保持邊界屬性的總變分(TV)約束，能夠更加有利于一維層狀速度模型的反演[9]。一階TV正則化約束項(xiàng)用相鄰細(xì)層的速度差表示為：

(6)

(7)

其中：m表示反演的速度模型參數(shù)；n表示速度模型參數(shù)的個(gè)數(shù)；Lm近似為反演模型的一階偏微分近似[10]。最終的目標(biāo)函數(shù)表示為式(8)。

OF=‖G(m)-d‖1+α‖Lm‖1

(8)

其中：G(m)表示反演速度模型獲得的走時(shí)差；d表示觀測(cè)得到的到時(shí)差；α表示權(quán)重因子；‖Lm‖1表示TV正則化約束項(xiàng)。

筆者利用田玥等[11]提出的快速兩點(diǎn)間射線追蹤方法，計(jì)算層狀速度模型中源到檢波器的走時(shí)信息?；谀Ｐ蚅1范數(shù)的問(wèn)題可以采用內(nèi)點(diǎn)方法(interior points solver)求解[12]，這里借助于CVX軟件包[13]求解目標(biāo)函數(shù)。

2 人工合成數(shù)據(jù)測(cè)試

對(duì)四層速度模型進(jìn)行了測(cè)試。觀測(cè)系統(tǒng)由20級(jí)布設(shè)于垂直觀測(cè)井的檢波器構(gòu)成，檢波器之間的間距為10 m。射孔點(diǎn)與檢波器的水平距離為500 m。初始速度模型由10 m 間隔，每個(gè)細(xì)層中P波或S波速度值大小均相同的細(xì)層構(gòu)成。反演過(guò)程中，我們進(jìn)行了20次迭代優(yōu)化目標(biāo)函數(shù)。

圖1 四層速度模型反演結(jié)果Fig.1 The inverted result of a four-layer velocity model(a)P波：3 500 m/s、S 波：1 750 m/s; (b)P波：4 000 m/s、S 波：2 000 m/s;(c)P波：4 500 m/s、S波：2 250 m/s

圖2 三種初始速度模型目標(biāo)函數(shù)收斂速度Fig.2 The convergence rate of three initial velocity models

圖1測(cè)試了三種不同初始速度模型的校正結(jié)果(實(shí)線表示真實(shí)速度模型，虛線表示初始速度模型，點(diǎn)線表示反演速度模型)。從圖1可以看出，最終反演結(jié)果受初始速度模型影響較小。反演的各細(xì)層速度值與真實(shí)速度值的偏差都不超過(guò)1%。從圖2可以看出(實(shí)心圓表示圖1(a)，五角星表示圖1(b),空心圓表示圖1(c)，經(jīng)過(guò)5次迭代以后，目標(biāo)函數(shù)已經(jīng)收斂接近于“0”。說(shuō)明利用本文發(fā)展的方法，可以快速地獲得層狀速度模型的層位個(gè)數(shù)、層位深度以及各層速度值的大小。

針對(duì)圖1中的四層速度模型，測(cè)試了在較少檢波器情況下速度模型的校正情況。圖3(a)顯示了四層模型在10個(gè)檢波器監(jiān)測(cè)情況下速度模型的校正情況。通過(guò)反演的速度模型可以看出，存在檢波器的層位(前10層)，反演的速度模型各細(xì)層偏離真實(shí)速度值的大小不超過(guò)2%，能夠通過(guò)反演的速度值大小劃分層位，獲得各層速度值。不存在檢波器的層位(后10層)，不能夠分辨出正確的層位深度，這是由于單個(gè)射孔信號(hào)的觀測(cè)數(shù)據(jù)包含的信息不足造成的。為了增加地震射線路徑在速度模型中的覆蓋程度，通過(guò)采用多個(gè)射孔點(diǎn)數(shù)據(jù)同時(shí)反演以增強(qiáng)對(duì)細(xì)層的分辨能力。圖3(b)表示采用三個(gè)射孔數(shù)據(jù),各射孔點(diǎn)水平距離大約為25 m，的速度模型校正結(jié)果。各細(xì)層偏離真實(shí)速度值的大小均不超過(guò)2%，表明增加射孔點(diǎn)的個(gè)數(shù)，增加了射線路徑在速度模型中的覆蓋程度，有利于反演的速度模型準(zhǔn)確性地提高。為了模擬獲得實(shí)際信號(hào)過(guò)程中的隨機(jī)誤差，在合成的到時(shí)信息中加入平均值為“0”，標(biāo)準(zhǔn)差為 1 ms的高斯噪聲，測(cè)試本文方法在含噪聲數(shù)據(jù)中反演的穩(wěn)定性。圖3(c)表明添加噪聲后反演的速度模型，含有噪聲情況下反演的P波和S波速度模型各細(xì)層偏離真實(shí)值均不超過(guò)2%，表明該方法具有良好的抗噪性，在噪聲情況下仍然能夠獲得正確的層位深度以及各層速度值的大小。

圖3 十個(gè)檢波器速度模型校正結(jié)果Fig.3 The inverted results of 10-receiver velocity model(a)單個(gè)射孔點(diǎn)校正結(jié)果；(b)三個(gè)射孔點(diǎn)沒(méi)有噪聲數(shù)據(jù)校正結(jié)果；(c)三個(gè)射孔點(diǎn)存在噪聲數(shù)據(jù)校正結(jié)果。

3 實(shí)際資料處理

將該方法應(yīng)用到遼河油田某水力壓裂監(jiān)測(cè)資料中，測(cè)試該方法的實(shí)際應(yīng)用效果。采用遼河油田某次水力壓裂單級(jí)射孔數(shù)據(jù)，進(jìn)行速度模型校正。該次壓裂具有三簇射孔段，射孔點(diǎn)位于檢波器的下方。三個(gè)射孔點(diǎn)布設(shè)于水平井距離地面大約1 550 m的深度左右，水平間隔大約為25 m。井中檢波器組由10級(jí)三分量檢波器組成，間隔為10 m，布設(shè)于距離地面1 425 m～1 515 m深度。射孔段與檢波器的水平距離大約為400 m。首先通過(guò)聲波測(cè)井曲線劃分速度模型(圖4(a))，通過(guò)該速度模型對(duì)射孔點(diǎn)進(jìn)行重定位，發(fā)現(xiàn)三個(gè)射孔點(diǎn)在不考慮方位角差異的二維垂直平面上的重定位誤差在350 m～430 m范圍，存在較大誤差(表1)。通過(guò)本文方法利用間隔為10 m、垂向深度為200 m的均勻細(xì)層初始速度模型，結(jié)合射孔數(shù)據(jù)對(duì)速度模型進(jìn)行校正，獲得多層速度模型(實(shí)線表示反演速度模型，圖4(b))。利用校正的速度模型對(duì)射孔點(diǎn)進(jìn)行重定位，發(fā)現(xiàn)重定位后三個(gè)射孔點(diǎn)，在不考慮方位角差異的二維垂直平面上的定位誤差分別為10.1 m、2.6 m以及8.1 m，大大縮小了射孔點(diǎn)的重定位誤差(表1)。說(shuō)明利用本方法可以有效地獲得地下速度模型的層位信息和各層速度值的大小，同時(shí)能夠有效地減少射孔點(diǎn)的重定位誤差。利用該速度模型對(duì)此次水力壓裂微地震事件進(jìn)行定位處理，獲得的定位結(jié)果如圖5所示(不同顏色代表不同壓裂段)。由圖5可以看出，微地震事件之間具有良好的連續(xù)性，能夠反映水力壓裂人工裂縫的主要延伸方向。由于水力壓裂產(chǎn)生的人工裂縫預(yù)期主要沿最大主應(yīng)力方向延伸，將微地震定位獲得的裂縫方位角與地質(zhì)資料提供的該地區(qū)最大主應(yīng)力方向進(jìn)行對(duì)比(圖6)，發(fā)現(xiàn)兩者具有高度的一致性，驗(yàn)證了此次微地震監(jiān)測(cè)成果的正確性，說(shuō)明本方法在真實(shí)微地震監(jiān)測(cè)中具有良好的應(yīng)用價(jià)值。

4 結(jié)論

為了獲得井中微地震監(jiān)測(cè)速度模型的層位個(gè)數(shù)、層位深度以及各層速度值大小，我們采用P波 走時(shí)差、S波走時(shí)差、P/S波走時(shí)差作為觀測(cè)數(shù)據(jù)， 在TV正則化約束條件下校正目標(biāo)層位速度模型的 層位信息和速度值的大小：

圖4 兩種速度模型對(duì)比Fig.4 The comparison between two velocity models(a)根據(jù)聲波測(cè)井曲線劃分的速度模型；(b)根據(jù)本文方法反演獲得的P波速度模型和S波速度模型

圖5 水力壓裂微地震監(jiān)測(cè)定位結(jié)果Fig.5 The fracture mapping of the microseismic monitoring in the hydraulic fracturing

射孔編號(hào)聲波測(cè)井速度模型校正速度模型ΔH/mΔZ/mΔL/mΔH/mΔZ/mΔL/m1-2.3350.85357.99.991.4510.12137.81350.82376.90.902.412.63-27.81422.24423.27.513.028.1

ΔH表示水平誤差，ΔZ表示垂直誤差，ΔL表示總誤差

1)該方法有效解決了射孔信號(hào)起始時(shí)間未知對(duì)速度模型反演的影響。同時(shí)用均勻細(xì)層速度模型作為初始速度模型，有效避免了聲波測(cè)井信息不準(zhǔn)對(duì)速度模型層位深度和速度值大小的影響。

2)利用射孔數(shù)據(jù)校正速度模型過(guò)程中，在TV正則化約束條件下，有效地保護(hù)了速度模型的層位信息，同時(shí)獲得了層位深度和各層速度值的大小。

3)合成數(shù)據(jù)和真實(shí)資料均顯示，利用本文方法反演獲得的速度模型與真實(shí)速度模型更加接近，可以獲得更好的射孔事件重定位結(jié)果，顯示該方法在微地震監(jiān)測(cè)資料處理中具有良好的實(shí)用性。

圖6 最大主應(yīng)力方向和裂縫方位角對(duì)比
Fig.6 The comparison between the maximum principal stress directions and the azimuths in the microseismic monitoring
(a)地質(zhì)資料提供的最大主應(yīng)力方向；(b)微地震監(jiān)測(cè)獲取的裂縫方位角