国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

測(cè)井?dāng)?shù)據(jù)自適應(yīng)多項(xiàng)式擬合去噪

2013-12-03 10:29:46蔡志東張固瀾劉振亞
測(cè)井技術(shù) 2013年1期
關(guān)鍵詞:厚層階次薄層

蔡志東,張固瀾,劉振亞

(東方地球物理公司新興物探開發(fā)處,河北 涿州072751)

0 引 言

受到地層中隨機(jī)噪聲、井徑大小、泥漿擾動(dòng)和儀器自身等諸多因素的影響,測(cè)井資料中總是存在著不同程度的噪聲。這些噪聲在測(cè)井曲線上表現(xiàn)為許多與地層屬性無關(guān)的統(tǒng)計(jì)起伏變化或是毛刺干擾[1]。測(cè)井曲線去噪通常采用曲線平滑、中值濾波、頻率域低通濾波和小波變換等方法。

曲線平滑、中值濾波是直接、有效的高頻噪聲去除方法,但容易使測(cè)井曲線對(duì)地層屬性的反映變模糊,降低屬性界面的定位精度,影響了薄層和薄互層的識(shí)別能力;頻率域?yàn)V波也是一種常用的去噪方法,通常采用低通濾波濾除高頻噪聲,但不利于局部屬性細(xì)節(jié)分析[2];小波變換被用于抑制測(cè)井?dāng)?shù)據(jù)噪聲,可以有效地對(duì)白噪聲進(jìn)行去除,但是對(duì)于脈沖噪聲的衰減有限,當(dāng)噪聲水平比較高時(shí)會(huì)使原始數(shù)據(jù)的高頻端信號(hào)受到損傷[3]。以上幾種方法都普遍應(yīng)用于測(cè)井?dāng)?shù)據(jù)的去噪處理之中,但都是基于對(duì)測(cè)井曲線整體去噪處理的方法,不能有區(qū)別地對(duì)薄、厚層分別進(jìn)行處理。真實(shí)的地層屬性變化多種多樣,同一口井所在的地層可能同時(shí)存在著大套的同屬性介質(zhì)、屬性漸變的介質(zhì)、薄層與薄互層介質(zhì)等情況。本文以滑動(dòng)多項(xiàng)式擬合法為基礎(chǔ),提出一種自適應(yīng)改變多項(xiàng)式階次進(jìn)行測(cè)井曲線擬合的方法,通過統(tǒng)計(jì)窗口內(nèi)樣點(diǎn)值的變化規(guī)律,評(píng)估出最優(yōu)多項(xiàng)式階次,用該階次多項(xiàng)式擬合窗口內(nèi)的測(cè)井曲線,取擬合曲線的中點(diǎn)進(jìn)行樣點(diǎn)替換,然后向后移動(dòng)一個(gè)樣點(diǎn)重復(fù)上面過程,最終完成測(cè)井曲線擬合重構(gòu)。

1 方法原理

假設(shè)采集測(cè)井記錄用X(x)表示,噪聲成分用S(x)表示,實(shí)際測(cè)井記錄用Y(x)表示,則

式中,x為測(cè)井記錄中離散的數(shù)據(jù)樣點(diǎn)。如果利用多項(xiàng)式對(duì)全部測(cè)井曲線Y(x)進(jìn)行擬合,則

式中,n為多項(xiàng)式階次,當(dāng)n=0時(shí),Y(x)為常數(shù)a0,適合用來擬合屬性不變的厚層;當(dāng)n=1,Y(x)適合用來擬合屬性線性漸變的厚層;當(dāng)n≥2時(shí),Y(x)則可用來擬合屬性非線性漸變的厚層、薄層和薄互層等,式(2)可以通過最小二乘法求得多項(xiàng)式Y(jié)(x)的系數(shù)ai。

當(dāng)采用逐個(gè)樣點(diǎn)連續(xù)滑動(dòng)多項(xiàng)目式擬合時(shí),以m表示滑動(dòng)窗口內(nèi)的樣點(diǎn)數(shù)(為了使擬合數(shù)據(jù)中點(diǎn)對(duì)原數(shù)據(jù)中點(diǎn)的替換,約定m為奇數(shù);程序?qū)崿F(xiàn)時(shí),如果判斷m為偶數(shù)時(shí)將自動(dòng)增加1個(gè)樣點(diǎn)),則第j個(gè)m長(zhǎng)度的窗口內(nèi)擬合多項(xiàng)式為

式中,nj為第j個(gè)擬合多項(xiàng)式的階次;l為測(cè)井?dāng)?shù)據(jù)記錄的總樣點(diǎn)數(shù)。該段擬合多項(xiàng)式的自由度vj為

為了階次nj可以自適應(yīng)判斷選擇,引入了卡方檢驗(yàn)法,式(5)用來計(jì)算χ2值。

式中,mj和m′j分別為第j個(gè)m窗口的第1個(gè)和最后一個(gè)樣點(diǎn)序號(hào),均值則卡方檢驗(yàn)函數(shù)為

式中,Γ(z)為Γ函數(shù),即用 Δj表示擬合的良好度,則

分別對(duì)第j個(gè)窗口內(nèi)0階到最高階的擬合的良好度Δj進(jìn)行對(duì)比判斷,即可優(yōu)選出nj,此時(shí)擬合多項(xiàng)式為Yj(x),當(dāng)x=xj時(shí)取到該多項(xiàng)式的中點(diǎn)Yj(xj),即可得到最終的擬合測(cè)井曲線為

式(8)中邊界部分采用首尾多項(xiàng)式替換的方法,通過這種自適應(yīng)多項(xiàng)式擬合的方法可以有區(qū)別地對(duì)薄層和厚層進(jìn)行擬合,使測(cè)井曲線更能反映實(shí)際地層的屬性。

2 仿真分析

2.1 仿真數(shù)據(jù)生成

根據(jù)上述理論,為了確定自適應(yīng)多項(xiàng)式擬合測(cè)井曲線的參數(shù)和分析去噪效果,首先利用模型數(shù)據(jù)進(jìn)行仿真試驗(yàn)。模型測(cè)井?dāng)?shù)據(jù)由隨機(jī)變化的模擬測(cè)井曲線和隨機(jī)噪聲2種組分合成,數(shù)據(jù)采樣間隔為0.1m、記錄長(zhǎng)度100m。圖1(a)為模型數(shù)據(jù)的第1種組分隨機(jī)噪聲;圖1(b)為模擬數(shù)據(jù)的第2種組分模擬測(cè)井曲線;圖1(c)為由2種組分所合成的模型數(shù)據(jù)。

圖1 模型數(shù)據(jù)的生成

2.2 擬合多項(xiàng)式的階次范圍

自適應(yīng)滑動(dòng)多項(xiàng)式擬合的最低擬合階次為0階,即允許利用常數(shù)直線擬合測(cè)井記錄,這是與記錄中同屬性厚層的測(cè)井曲線特征相符合的;同時(shí),也對(duì)最高擬合階次進(jìn)行了限制,如果階次過高,不僅將影響計(jì)算效率,又可能對(duì)有效數(shù)據(jù)薄層識(shí)別造成傷害。

圖2為隨機(jī)截取的一段10m長(zhǎng)的模型數(shù)據(jù),利用0~10階多項(xiàng)式分別進(jìn)行擬合。結(jié)果表明,隨著階次的增大,擬合曲線的波動(dòng)隨之增強(qiáng)。

通過不同時(shí)長(zhǎng)的數(shù)據(jù)進(jìn)行測(cè)試,并結(jié)合實(shí)際數(shù)據(jù)處理經(jīng)驗(yàn)和計(jì)算效率分析,認(rèn)為以6~8階為最高擬合階次比較合適。

圖2 模型數(shù)據(jù)(10m)不同階次的多次項(xiàng)式擬合

2.3 擬合優(yōu)良度門檻值的確定

擬合優(yōu)良度門檻值是用來判斷多項(xiàng)式擬合效果的一個(gè)臨界值,用來自適應(yīng)的優(yōu)選多項(xiàng)式階次,這個(gè)參數(shù)可以不進(jìn)行設(shè)置,默認(rèn)取優(yōu)良度最好的階次進(jìn)行擬合。

如果給定了擬合優(yōu)良度門檻值,當(dāng)所有階數(shù)的擬合優(yōu)良度均不能滿足該門檻值時(shí),以優(yōu)良度最高的階次進(jìn)行擬合,效果等同于未設(shè)置該參數(shù);當(dāng)只有1個(gè)階次的擬合優(yōu)良度高于這一門檻值時(shí),就以這個(gè)階次進(jìn)行多項(xiàng)式擬合;當(dāng)存在多個(gè)階次同時(shí)滿足門檻值時(shí),將選擇最低階數(shù)的多項(xiàng)式擬合,這樣既減少了計(jì)算量也更符合測(cè)井曲線的地球物理特征。

圖2所示模型數(shù)據(jù)的擬合優(yōu)良度如圖3所示,隨著擬合階次的增加,擬合的優(yōu)良度呈現(xiàn)遞增趨勢(shì),建議該門檻值設(shè)置不小于80%。

2.4 窗口長(zhǎng)度的確定

圖3 模型數(shù)據(jù)(10m)不同階次的擬合優(yōu)良度

模型數(shù)據(jù)采樣間隔為0.1m,如果要識(shí)別的薄層厚度為1m,假定擬合多項(xiàng)式的最高階次可以取到8次,則窗口樣點(diǎn)數(shù)最大可以取到70個(gè),超過這一數(shù)值將影響薄層的識(shí)別效果;同時(shí)要保證對(duì)高頻噪聲的去除效果,窗口樣點(diǎn)數(shù)也不能過??;通常選用與曲線平滑、中值濾波等計(jì)算方法相當(dāng)或略大的窗口長(zhǎng)度即可,建議窗口長(zhǎng)度對(duì)應(yīng)測(cè)井深度在1~5m之間,對(duì)于不同的測(cè)井?dāng)?shù)據(jù),要根據(jù)具體情況在這一范圍內(nèi)進(jìn)行多次實(shí)驗(yàn)確定最佳參數(shù)。

2.5 衰減結(jié)果分析

利用前面所分析確定的參數(shù),對(duì)模型數(shù)據(jù)進(jìn)行了去噪處理。圖4(a)為自適應(yīng)滑動(dòng)多項(xiàng)式測(cè)井曲線擬合,顯示了擬合過程中用到的不同階次的所有多項(xiàng)式;圖4(b)為原始測(cè)井?dāng)?shù)據(jù)和自適應(yīng)擬合后的測(cè)井?dāng)?shù)據(jù);圖4(c)為擬合數(shù)據(jù)與合成前不含噪聲數(shù)據(jù)的差值。處理結(jié)果表明,擬合效果較好,可以去噪還原得到較可靠的測(cè)井曲線,從殘差結(jié)果看處理誤差大多數(shù)在10%以內(nèi),可以驗(yàn)證該方法的有效性。

對(duì)于實(shí)際測(cè)井資料,為了獲得最好的去噪效果,擬合多項(xiàng)式階次范圍、卡方檢驗(yàn)優(yōu)良度門檻值和窗口長(zhǎng)度這3個(gè)參數(shù)在一定范圍內(nèi)可以適當(dāng)?shù)剡M(jìn)行調(diào)整。

圖4 模型數(shù)據(jù)的自適應(yīng)滑動(dòng)多項(xiàng)式擬合去噪

3 實(shí)際資料處理應(yīng)用

3.1 原始測(cè)井?dāng)?shù)據(jù)分析

以土庫曼斯坦巴格德雷地區(qū)1口井的測(cè)井資料為實(shí)例進(jìn)行分析。該井區(qū)地層平緩、屬性穩(wěn)定,從錄井資料分析,該井段主要以層狀灰色泥巖為主,夾著多套厚度不大于5m的淺灰色粉砂巖薄層。圖5(a)為原始自然伽馬測(cè)井曲線,數(shù)據(jù)的噪聲水平較高,需要經(jīng)過去噪處理才能反映較真實(shí)的地層屬性。

圖5 實(shí)際測(cè)井資料的各種去噪方法對(duì)比

3.2 常規(guī)去噪效果

圖5(b)、圖5(c)和圖5(d)是分別采用了移動(dòng)平均平滑、中值濾波、頻率域低通濾波等3種方法去噪得到的處理結(jié)果,3種方法均能達(dá)到一定的高頻噪聲去除效果。綜合對(duì)比3種方法,移動(dòng)平均光滑效果不如其他2種方法;中值濾波法的較厚地層屬性界面略清晰,但薄層反映不明顯;頻率域低通濾波法的對(duì)薄層有一定反映,但地層屬性界面卻較不清晰。

3.3 新方法去噪效果

圖5(e)采用的是自適應(yīng)滑動(dòng)多項(xiàng)式擬合測(cè)井?dāng)?shù)據(jù)去噪的方法,擬合階次范圍為0~8階次、窗口長(zhǎng)度為2m(21個(gè)樣點(diǎn))、擬合優(yōu)良度門檻值為80%,從去噪效果看,其地層屬性界面清晰、薄層反應(yīng)良好,總體優(yōu)于其它幾種常規(guī)的去噪效果,通過與圖5(f)錄井剖面的對(duì)比可知地層屬性反應(yīng)良好,層界面位置準(zhǔn)確。

4 結(jié) 論

(1)提出了一種全新的利用自適應(yīng)滑動(dòng)多項(xiàng)式擬合測(cè)井曲線去噪方法,不需要進(jìn)行頻率域轉(zhuǎn)換,可以對(duì)目前測(cè)井?dāng)?shù)據(jù)中所存在的噪聲進(jìn)行有效去除,測(cè)井記錄的數(shù)據(jù)量較小,多次循環(huán)運(yùn)算對(duì)處理效率的影響不大。

(2)自適應(yīng)的階數(shù)選擇可以自行判斷并區(qū)別對(duì)薄層和厚層進(jìn)行不同階次的測(cè)井曲線擬合,從而達(dá)到去噪目的。

(3)該方法已應(yīng)用于土庫曼斯坦巴格德雷地區(qū)多口測(cè)井的去噪處理試驗(yàn),在實(shí)際生產(chǎn)中應(yīng)用效果較好。

[1] 雍世和,張超模.測(cè)井?dāng)?shù)據(jù)處理與綜合解釋[M].東營(yíng):中國石油大學(xué)出版社,1996.

[2] 鄒長(zhǎng)春,楊欣德,潘令枝,等.一種基于小波變換的測(cè)井曲線去噪新方法[J].物探與化探,1999,6(23):462-466.

[3] 司清亮.基于四種改進(jìn)閾值函數(shù)的小波去噪[J].河南師范大學(xué)學(xué)報(bào):自然科學(xué)版,2009,5(37):140-142.

[4] 施振飛,許海萍,高宏祥,等.提高薄層測(cè)井資料解釋分辨率的新方法[J].測(cè)井技術(shù),2001,2(25):131-135.

[5] 劉詠梅,錢真,申麗然,等.一種新的非線性曲線平滑方法及在測(cè)井曲線識(shí)別中的應(yīng)用[J].模式識(shí)別與人工智能,2005,4(18):447-452.

猜你喜歡
厚層階次薄層
溫州地區(qū)厚層軟土UU三軸試驗(yàn)與FLAC3D數(shù)值模擬
階次分析在驅(qū)動(dòng)橋異響中的應(yīng)用
基于Vold-Kalman濾波的階次分析系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)*
Winter Walking in Wales威爾士的冬季徒步旅行
基于齒輪階次密度優(yōu)化的變速器降噪研究
維藥芹菜根的薄層鑒別
SiN_x:H膜沉積壓強(qiáng)與擴(kuò)散薄層電阻的匹配性研究
厚層塊狀特低滲礫巖油藏水平井壓裂參數(shù)優(yōu)化
參芪苓口服液的薄層色譜鑒別
哮天龍水庫大壩厚層碾壓混凝土技術(shù)研究
沙坪坝区| 南城县| 巩义市| 延寿县| 柘荣县| 九龙坡区| 湘潭市| 克拉玛依市| 嘉黎县| 丹阳市| 纳雍县| 长乐市| 铜川市| 奉贤区| 新津县| 陇西县| 辽阳市| 五华县| 桃园市| 礼泉县| 云龙县| 利津县| 将乐县| 北安市| 苗栗县| 佛冈县| 武夷山市| 肃北| 巴塘县| 武宁县| 凯里市| 乐亭县| 公安县| 博乐市| 全州县| 宜丰县| 新干县| 江永县| 汉寿县| 平邑县| 胶州市|