油藏電纜地層測試壓力信號的小波分析

油藏電纜地層測試壓力信號的小波分析

湯小燕1，劉之的2，高曦3

(1.西安科技大學(xué) 地質(zhì)與環(huán)境學(xué)院，陜西 西安 710054；2.西安石油大學(xué) 地球科學(xué)與工程學(xué)院，陜西 西安 710065；3.中國石油長慶油田第九采油廠，寧夏 銀川 750006)

摘要：為了抑制非地層壓力等噪聲信號，有效地提取與產(chǎn)能評價密切相關(guān)的油藏壓力信息，詳細(xì)地探討了電纜地層測試的流動方式及物理方程之后，基于小波分析基本原理，深入剖析了電纜地層測試壓力信號的小波分析方法。利用小波分析高分辨率的特點，采用小波變換分析法從不同尺度上提取了電纜地層測試壓力的小波系數(shù)。結(jié)果表明，經(jīng)小波變換處理后的電纜地層壓力測試曲線能夠有效地反映地層中的流體信息；不同尺度上提取的電纜地層壓力的小波系數(shù)中，高頻小波系數(shù)d1和d2包含的產(chǎn)能及滲透率等信息更強；采用小波變換后降低了干擾信號的影響，突出了油藏產(chǎn)能信號，進(jìn)而提高了采樣數(shù)據(jù)的可靠性。這也充分說明利用小波分析方法對油藏電纜測試壓力信號處理后，油藏產(chǎn)能信息更為明晰，進(jìn)而提高了油藏電纜地層測試壓力有用信號的提取精度，且該法處理方式更為簡單，為今后油藏電纜地層測試壓力數(shù)據(jù)處理方法提供了一種新手段。

關(guān)鍵詞：壓力；電纜地層測試；小波分析

DOI:10.13800/j.cnki.xakjdxxb.2015.0216

文章編號：1672-9315(2015)02-0229-07

收稿日期:*2014-12-25責(zé)任編輯:李克永

基金項目：陜西省自然科學(xué)基礎(chǔ)研究計劃項目(2013JQ5008)

通訊作者：湯小燕(1977-)，女，四川廣安人，博士，講師，E-mail:lzdtxy2004@sina.com

中圖分類號：TE 122文獻(xiàn)標(biāo)志碼： A

Waveletanalysisforreservoirpressuresignalofwirelineformationtest

TANGXiao-yan1,LIU Zhi-di2,GAO Xi3

(1.College of Geology and Environment,Xi’an University of Science and Technology,Xi’an 710054,China；

2.School of Earth Sciences and Engineering，Xi’an Shiyou University，Xi’an 710065，China；

3.Oil Production Plant 9 of Petro China Changqing Oilfield Company,Yinchuan 750006,China)

Abstract：In order to restrain noise such as the formation pressure，and effectively abstract pressure information which is closely related to the production evaluation，the flow pattern and the physical equations of the wireline formation testing are derived.The wavelet analysis method for wireline formation test pressure is analyzed using the method of wavelet transform.In accordance with high-resolution features of wavelet analysis，the wavelet coefficients of the wireline formation pressure are extracted from different scales.The results show that the test curve of wireline formation pressure reflects the information flow of formation.The wavelet coefficients of high frequency wavelet coefficients of d1 and d2 in different scales from wireline formation pressure contains information such as the capacity and permeability stronger；The effects of interfering signal is reduced using wavelet transform，and the reservoir capacity signal is highlighted，thus the reliability of the sampling data is improved.Through wavelet decomposition for reservoir wireline formation test pressure data，the reservoir capacity information are more clear.Thereby the extraction accuracy of useful signal for reservoir wireline formation pressure test is improved，and the method is simple，it provides a new means for the data processing of wireline test pressure of the reservoir.

Key words：pressure；wireline formation test；wavelet analysis

0引言

在石油勘探開發(fā)中，電纜地層測試壓力是油藏產(chǎn)能評價的重要參數(shù)之一，但由于非地層壓力因素的存在會對油藏產(chǎn)能評價帶來誤差[1]。因此，在對壓力信號進(jìn)行相關(guān)性分析時，如何有效地提取有用的壓力信號是關(guān)鍵。油藏電纜地層測試是在地下復(fù)雜的地質(zhì)環(huán)境中測量的，測試的是壓力隨時間變化的序列，其壓力變化反映了不同的地質(zhì)條件[2]。該壓力變化不僅包含有用的地質(zhì)信息，而且還混雜著很多干擾信息。壓力曲線的明顯波動反映了壓力測試過程中的壓力波動或者壓力干擾，壓力波動有時即表現(xiàn)為干擾信息[3]。

地層壓力信號和噪聲混雜在一起，常用的濾波方法所得到的信號嚴(yán)重失真，失去了原信號所攜帶的地層壓力信息，從而影響信號的處理與分析。小波濾波能較好地解決該問題，小波分析作為一種新的數(shù)學(xué)工具，目前在信號處理、語音圖像分析、數(shù)據(jù)壓縮等眾多方面有顯著發(fā)展[4-6]，在油藏工程領(lǐng)域也得到了廣泛的應(yīng)用，如采用小波分析來完成濾波去噪等[1]。因此，有效抑制地層測試壓力過程中干擾信號、濾除混雜的噪聲、還原地層壓力信號的工作顯得尤為重要。

在信號分析領(lǐng)域，統(tǒng)計分析方法主要揭示信號在時域方面的信息，無法分析信號頻域方面的信息[7-8]；傅立葉分析方法則主要從頻域方面揭示信息特征，于是會丟失時域信息[9]。然而，小波分析方法克服上述方法缺陷的同時，可以提供原始信號的時域和頻域信息[10-12]。由于小波函數(shù)的時頻分辨能力較強，于是它成為信號分析與處理的一種良好工具[13-14]。小波分析的突出優(yōu)點是可以通過小取樣步長刻畫信號高頻成分的任意微小細(xì)節(jié)，該特性也被稱為時頻局部化性質(zhì)，鑒于此，小波分析被譽為數(shù)學(xué)“顯微鏡”[15-16]。油藏電纜地層測試的原始曲線是測試過程中壓力變化與時間的關(guān)系，可以把這個壓力變化看成是信號與時間之間的關(guān)系。電纜地層測試得到的壓力曲線包含了地質(zhì)信息，對該壓力曲線進(jìn)行小波變換分解，可得到不同分辨率下的小波系數(shù)，這些信息將反映電纜地層測試過程中壓力的變化特征，進(jìn)而反映了地層流體在流動過程中壓力及產(chǎn)量變化的信息。為此，文中利用小波分析高分辨率的特點，提取電纜地層壓力測試資料的細(xì)節(jié)特征，并使用小波系數(shù)刻畫壓力測試曲線的細(xì)節(jié)特征，以期為油藏產(chǎn)能評價提供可靠的地層壓力資料。

1電纜地層測試器的流動方式及物理方程

1.1電纜地層測試器的流動方式

地層中的流體在向井筒流動過程中，由于油藏條件、流速和流量的變化，流體的流動壓力與時間的關(guān)系可以分為穩(wěn)態(tài)流、半穩(wěn)態(tài)流和非穩(wěn)態(tài)流，用公式(1)表示這3種流動狀態(tài)[17]

(1)

式中p為流體流動壓力，psi；t為流動時間，d.

電纜地層測試器主要用于套管井，在產(chǎn)層厚度與封隔器間隔差距較小的情況下，地層流體主要以理想徑向流的方式流入井筒，然后進(jìn)入電纜地層測試儀器的取樣筒[2](圖1)。

圖1　理想徑向流的流動狀態(tài) Fig.1　Flow state of the ideal radial flow

1.2理想徑向流的物理方程

電纜地層測試的過程，實際上就是地層流體向取樣筒的流動過程[2]。聯(lián)立物質(zhì)連續(xù)方程、流動方程及壓縮方程，可以得到流體流動的一般描述方程[11]

(2)

式中r為距井軸的半徑，ft；K為單相流體在地層中的滲透率，mD；μ為流體粘度，cP；ρ為流體的密度，lb/ft3；φ為孔隙度，Decimal；Cf為地層的壓縮系數(shù)，1/psi.

考慮地層流體為微壓縮性和可壓縮性，分別對應(yīng)于原油和天然氣，方程(2)可以改寫為

(3)

式中Ct為地層與孔隙流體的綜合壓縮系數(shù)，1/psi；C為孔隙流體的壓縮系數(shù)，1/psi.

方程(3)適用于測試過程中原油樣品取樣的流動過程。

對于天然氣樣品測試采樣的流動過程

(4)

方程(4)與方程(2)聯(lián)立，得到天然氣采集樣品流動過程的描述方程

(5)

式中z為氣體壓縮因子，無量綱。

根據(jù)測試過程中測試條件對應(yīng)的邊界條件和初始條件，求解方程(3)和方程(5)，就可分別得到原油采樣和天然氣采樣流動過程的壓力變化表達(dá)式。

由于測試過程很短，只有幾分鐘到幾個小時，上述條件基本滿足，于是方程(3)的恒定邊界流量解為

(6)

式中p(r,t)為在t小時后，半徑為r處的壓力，psi；Q0為產(chǎn)量，bbl/d.

E(-x)稱為冪積分函數(shù)，在特定的x數(shù)值區(qū)間，可以使用近似方法計算，也可以使用方程(6)中的級數(shù)方法計算。

盡管方程(6)是在天然氣的前提下得到的。但對方程(6)進(jìn)行數(shù)值模擬發(fā)現(xiàn)，該方程可以適用于油氣兩相態(tài)。該方程中的參數(shù)，如μ,B0,Ct，z，都是溫度和壓力的函數(shù)，應(yīng)用已有的關(guān)聯(lián)方程可計算這些參數(shù)[2，17]，得到Q0-Δp關(guān)系圖(圖2)，當(dāng)參數(shù)ng(兩相態(tài)中氣態(tài)的摩爾分?jǐn)?shù)，完全氣態(tài)時，ng=1)比較小時(ng=0.15)，即以原油為主時，方程(6)表明，流體產(chǎn)量與生產(chǎn)壓差是線性關(guān)系，這一點與原油流動方程一致；在生產(chǎn)壓差大于35MPa以后，流體產(chǎn)量與生產(chǎn)壓差是下彎曲的非線性關(guān)系，即當(dāng)生產(chǎn)增大時，流體產(chǎn)量的增加趨勢逐漸變緩，這一點與天然氣的流動方程一致。圖2還顯示了凝析氣藏的特點，即當(dāng)生產(chǎn)增加時，流體中的部分物質(zhì)轉(zhuǎn)化成氣態(tài)，所以，流動過程呈現(xiàn)氣態(tài)的特征。當(dāng)ng增加為0.25時，方程(6)在30MPa壓差的時候，就開始顯示出氣態(tài)的流動特征了。但是，當(dāng)ng接近于1(ng=0.85)時，即接近于純氣態(tài)時，當(dāng)方程(6)計算的Q0-Δp關(guān)系無規(guī)律可循。其原因是，現(xiàn)有的狀態(tài)參數(shù)估算算法可能不完全適用于純氣態(tài)的流動方程。

圖2　方程(6)的數(shù)值模擬結(jié)果 Fig.2　Results of numerical simulation of equation(6)

對于可壓縮流體(如凝析天然氣藏)，需要求解方程(5)。其假設(shè)條件同原油油藏。但是，方程(5)的求解過程比較復(fù)雜，首先需要將方程(5)線性化。假設(shè)真實氣體位能m(p)為

(7)

對方程(7)微分，得到

(8)

聯(lián)立方程(8)和方程(5)，得到

(9)

同樣在恒流量邊界條件下，求解方程(9)，可以得到

(10)

式中pwf為井底的流動壓力，psi；pi為原始?xì)獠貕毫?，psi；Q為氣體產(chǎn)量，103ft3/d；pSC為標(biāo)準(zhǔn)壓力，psi；TSC為標(biāo)準(zhǔn)溫度，℃；h為油藏溫度，℃；rw為井眼半徑，ft；h為油藏厚度，ft；μi為原始壓力天然氣的粘度，cP；Ct為綜合壓縮系數(shù)，psi-1.

方程(10)可以計算出m(pwf)與測試產(chǎn)量之間的關(guān)系。使用方程(7)計算出m(p)-p曲線，pwf是測試過程中的流動壓力，是已知的。由此可以內(nèi)插確定m(pwf)，即m(pwf)也是已知的。同樣，m(pi)也是已知的。

在電纜地層測試資料的解釋中，根據(jù)測試的pwf，pi數(shù)值，可以得到m(pwf)，m(pi)值。當(dāng)氣體粘度和壓縮因子取不同的數(shù)值(μ，z)時，可以把方程(7)編成程序自動計算。然后，把m(pwf)，m(pi)分別帶入方程(10)，可以計算儲層滲透率、儲層產(chǎn)能等參數(shù)。

方程(10)還可以近似的壓力顯式方程表示。

(11)

(12)

至此為止，方程(5)、方程(9)和方程(10)分別表示了原油測試過程和天然氣測試過程中，在穩(wěn)定徑向流條件下，流體流動壓力，油氣藏原始壓力和流體采樣產(chǎn)量之間的關(guān)系。這些構(gòu)成電纜地層測試方法評價產(chǎn)能的基礎(chǔ)[18]。

2小波分析原理

小波變換是將基本小波[也叫母小波(mother wavelet)]的函數(shù)ψ(t)作位移τ后，在不同尺度α下與待分析的信號x(t)作內(nèi)積

(13)

等效的頻率表示是

(14)

式中X(ω)，ψ(ω)分別是x(t),ψ(t)的傅立葉變換。

用鏡頭觀測目標(biāo)x(t)(也就是待分析信號)，ψ(t)代表鏡頭所起的作用(例如：濾波或卷積)。τ相當(dāng)于使鏡頭相對于目標(biāo)平行移動，a的作用相當(dāng)于鏡頭向目標(biāo)推進(jìn)或遠(yuǎn)離。由此可見，小波變換具有以下特點：具有多分辨率，也叫多尺度的特點，可以由粗略到精細(xì)地逐步觀察信號；可以看到基本頻率特性為ψ(ω)的帶通濾波器在不同尺度a下對信號作濾波；適當(dāng)?shù)剡x擇基本小波，使ψ(t)在時域上為有限支撐，ψ(ω)在頻域上也比較集中，便可以使WT在時、頻兩域都具有表征信號局部特征的能力，因此，有利于檢測信號的瞬態(tài)或奇異性。

小波基具有尺度因子a和平移因子τ兩個重要參數(shù)。a增大，則時窗伸展，頻寬收縮，帶寬變窄，中心頻率降低，而頻率分辨率增高；a減小則帶寬增加，中心頻率升高，時間分辨率增高而頻率分辨率降低。這與油藏電纜地層測試壓力分析問題中高頻信號持續(xù)時間短、低頻信號持續(xù)時間長的特性相吻合。因此，小波變換能夠有效地對信號的時頻特性進(jìn)行分析。

將任意空間中的函數(shù)z(t)在小波基下進(jìn)行展開，稱作函數(shù)z(t)的連續(xù)小波變換，其表達(dá)式為

(15)

式中ψ(t)為基本小波或者稱母小波，也叫小波基；a為尺度因子，a>0；τ為位移因子。

小波是函數(shù)空間L2(R)中滿足“容許性條件”的一個函數(shù)或者信號

(16)

式中ψ(w)為小波母函數(shù)ψ(t)的傅立葉變換；R*=R-{0}表示非零實數(shù)全體。

而對于任意的實數(shù)對(a，τ)，稱如下形式的函數(shù)

(17)

式(17)是由小波母函數(shù)ψ(t)生成的依賴于參數(shù)對(a，τ)的連續(xù)小波函數(shù)，簡稱小波。小波具有很多種形式，所以也就有很多種的小波變換。

3電纜地層測試壓力曲線分析

電纜地層測試得到的壓力曲線包含了地層參數(shù)信息，對該壓力曲線進(jìn)行小波變換分解，可得到不同分辨率下的小波系數(shù)，這些信息將反映電纜地層測試過程中，壓力的變化特征，進(jìn)而反映地下流體信息和地層特征參數(shù)信息。但是，電纜地層測試是在地下一個復(fù)雜的地質(zhì)環(huán)境中進(jìn)行測量的，其得到的壓力曲線不僅包含有用的地質(zhì)信息，而且還混雜著很多干擾信息。壓力曲線中的壓力波動，有時表現(xiàn)為干擾信息。如圖4所示，在19s時刻附近存在的壓力擾動即表現(xiàn)為壓力干擾項。在流動段，壓力曲線變化率的快慢直接反映了地下流體的產(chǎn)能和滲透率參數(shù)。

為了對電纜地層測試壓力信號進(jìn)行小波變換，采用式(18)所示的數(shù)值近似積分來代替小波變換中的積分式

(18)

式中ΔT為采用間隔；x(n)為信號序列。

給定尺度因子a，依次計算不同位移因子τ值下上式的乘積和，就可求得a值下的一組小波系數(shù)。給定不同的尺度因子，就得到電纜地層測試壓力信號在不同變換尺度下的小波系數(shù)。

圖3～圖8是通過Db5小波[2-3]，進(jìn)行不同尺度因子的小波變換分解得到的小波系數(shù)，其尺度因子分別為a=2,4,8,16,32.從圖4～圖8可見，隨著尺度因子的增大，不管那種小波，其小波系數(shù)曲線都將變得平滑而沒有特征，只有低尺度因子的曲線顯示了原始壓力曲線流動段壓力有瞬時劇烈變化的特征，在圖上的表現(xiàn)形式為小波系數(shù)曲線快速跳動，并且在圖上也表現(xiàn)出壓力擾動也會產(chǎn)生明顯的小波高頻系數(shù)。隨著尺度因子的加大，壓力擾動產(chǎn)生的高頻消息系數(shù)會很快衰竭，但是流動段產(chǎn)生的壓力高頻小波系數(shù)依然存在。由此可見，不同分辨率的尺度因子對電纜地層測試壓力數(shù)據(jù)進(jìn)行變換處理結(jié)果表明，小波分析方法能夠放大數(shù)據(jù)的細(xì)節(jié)。

圖3　電纜地層測試壓力曲線分析 Fig.3　Curve analysis of cable formation testing pressure

圖4　Db5小波的小波系數(shù)d 1(尺度因子a=2) Fig.4　Wavelet coefficients d 1 of wavelet Db5(Scale factor a=2)

圖5　Db5小波的小波系數(shù)d 2(尺度因子a=4) Fig.5　Wavelet coefficients d 2 of wavelet Db5(Scale factor a=4)

圖6　Db5小波的小波系數(shù)d 3(尺度因子a=8) Fig.6　Wavelet coefficients d 3 of wavelet Db5(Scale factor a=8)

圖7　Db5小波的小波系數(shù)d 4(尺度因子 a=16) Fig.7　Wavelet coefficients d 4 of wavelet Db5(Scale factor a=16)

圖8　Db5小波的小波系數(shù)d 5(尺度因子a=32) Fig.8　Wavelet coefficients d 5 of wavelet Db5(Scale factor a=32)

根據(jù)電纜地層測試的理論知識[2-3]，壓力曲線的下降段反映了測試點的產(chǎn)量信息，恢復(fù)段反映滲透率參數(shù)，說明小波變換得到的小波系數(shù)包含有這些參數(shù)信息。因此，通過小波分析得到的高頻小波系數(shù)能夠反映電纜地層測試點的地質(zhì)信息，也就是說可能從中提取產(chǎn)能和滲透率等參數(shù)信息，特別是高頻小波系數(shù)d1和d2，包含的能量信息更強一些。

電纜地層測試的是壓力隨時間變化的序列，壓力曲線隨時間的變化反映了不同的地質(zhì)條件。壓力線的陡度反映壓力變化率，而流動段的壓力曲線變化更直接反映了儲層的相應(yīng)參數(shù)變化。

典型的電纜地層測試壓力曲線包含有鉆井液柱靜壓力、井底流壓、原始地層壓力等信息，壓力曲線的變化過程反映了電纜地層測試過程中儀器的實施和流體流動的過程。但是，實測的壓力曲線不會像標(biāo)準(zhǔn)曲線那樣完整和理想，電纜地層測試過程中會有一些壓力擾動或者干擾，而壓力曲線的明顯波動反映了壓力測試過程中的壓力波動或者壓力干擾。電纜地層測試的整個測試過程的壓力曲線包含了壓力擾動、流動段和地層壓力信息。以定性分析為基礎(chǔ)，對比各種實際的非理想的壓力記錄曲線與理想曲線的差別，也可以初步判定地層滲透性的高低等問題。從小波分析圖上可見，不管是油層還是氣層，其壓力曲線的小波變換均在壓力曲線流動段和壓力恢復(fù)段有高頻信息成分，反映了壓力變化過程中不同的地質(zhì)信息。

4結(jié)論

1)電纜地層壓力測試曲線反映了地層中流體信息。不同尺度上提取的電纜地層壓力的小波系數(shù)中，高頻小波系數(shù)d1和d2包含的產(chǎn)能及滲透率等信息更強；

2)油藏電纜地層測試壓力信號小波分析結(jié)果表明，采用小波變換后降低了干擾信號的影響，突出了油藏產(chǎn)能信號，進(jìn)而提高了采樣數(shù)據(jù)的可靠性，這也充分說明利用小波分析方法對油藏電纜測試壓力信號進(jìn)行分析是合理的。

參考文獻(xiàn)References

[1] 劉之的，趙靖舟，高秋濤.利用MDT資料預(yù)測油氣產(chǎn)能[J].地質(zhì)科技情報，2013，32(1)：163-166.

LIU Zhi-di，ZHAO Jing-zhou，GAO Qiu-tao.Predicting oil and gas deliverability using MDT data[J].Geological Science and Technology Information，2013，32(1)：163-166.

[2]匡立春.電纜地層測試資料應(yīng)用導(dǎo)論[M].北京：石油工業(yè)出版社，2005.

KUANG Li-chun.An introduction to application of wireline formation test data[M].Beijing：Petroleum Industry Press，2005.

[3]張聰慧，劉樹鞏，李義.利用電纜地層測試資料進(jìn)行低滲儲層流度計算和產(chǎn)能預(yù)測[J].中國海上油氣，2013(1)：43-45.

ZHANG Cong-hui，LIU Shu-gong，LI Yi.A discussion on methods to determine reservoir mohukunbility and productivity by using wireline formation test data[J].China Offshore Oil and Gas，2013(1)：43-45.

[4]周筱媛.一種使用小波分解的多尺度盲圖像復(fù)原方法[J].西安科技大學(xué)學(xué)報，2011，31(3)：363-365.

ZHOU Xiao-yuan.Blind image restoration using wavelet decomposition[J].Journal of Xi’an University of Science and Technology，2011，31(3)：363-365.

[5]齊愛玲，白煥莉.一種基于ICA的機(jī)械缺陷超聲信號提取方法[J].西安科技大學(xué)學(xué)報，2012，32(3)：394-397.

QI Ai-ling，BAI Huan-li.A method of mechanical defect ultrasonic signal extraction Based on ICA[J].Journal of Xi’an University of Science and Technology，2012，32(3)：394-397.

[6]姜友誼.小波神經(jīng)網(wǎng)絡(luò)在人臉識別中的應(yīng)用[J].西安科技大學(xué)學(xué)報，2012，32(5)：652-657.

JIANG You-yi.Face recognition Based on wavelet neural network[J].Journal of Xi’an University of Science and Technology，2012，32(5)：652-657.

[7]張賢達(dá).現(xiàn)代信號處理[M].北京：清華大學(xué)出版社，2004.

ZHANG Xian-da.Modern signal processing[M].Beijing：Tsinghua University Press，2004.

[8]張春森，范金健，胡平波.小波變換和Arnold變換的數(shù)字水印技術(shù)[J].西安科技大學(xué)學(xué)報，2012，32(1)：95-100.

ZHANG Chun-sen，F(xiàn)AN Jin-jian，HU Ping-bo.Digital watermark technology based on wavelet and arnold transform[J].Journal of Xi’an University of Science and Technology，2012，32(1)：95-100.

[9]唐禮忠，陳資南.礦山微震信號小波分析與研究[J].科技導(dǎo)報，2013，31(32)：29-33.

TANG Li-zhong，CHEN Zi-nan.Research and analysis on wavelet of mine microseismic signals[J].Science & Technology Review，2013，31(32)：29-33.

[10] 薛軍鵬，范海宏，徐德龍.大顆粒鼓泡流化床壓力脈動信號的小波分析[J].化工裝備技術(shù)，2006，27(2)：23-27.

XUE Jun-peng，F(xiàn)AN Hai-hong，XU De-long.Wavelet analysis of pressure fluctuation signal of larger particles bubbling fluidized bed[J].Chemical Equipment Technology，2006，27(2)：23-27.

[11] 樊計昌，劉明軍，海燕.小波(包)濾波方法的GUI及其在深地震測深數(shù)據(jù)處理中的應(yīng)用[J].科技導(dǎo)報，2009，27(2)：78-82.

FAN Ji-chang，LIU Ming-jun，HAI Yan.GUI for De-noising in Wavelet(Packet)method and its application to DSS data processing[J].Science & Technology Review，2009，27(2)：78-82.

[12] 汪峰，武風(fēng)波，黃興.Hopfield網(wǎng)絡(luò)與HVS相結(jié)合的小波域圖像水印算法[J].西安科技大學(xué)學(xué)報，2013，33(6)：731-736.

WANG Feng，WU Feng-bo，HUANG Xing.Watermarking algorithms of wavelet image for hopfield network with HVS combination[J].Journal of Xi’an University of Science and Technology，2013，33(6)：731-736.

[13] 谷立臣，閆小樂，劉澤華.小波變換的流體壓力信號自適應(yīng)濾波方法研究[J].振動·測試與診斷，2009，29(2)：205-208.

GU Li-chen，YAN Xiao-le，LIU Ze-hua.Study adaptive filter method of wavelet transform for fluid pressure signal[J].Journal of Vibration，Measurement & Diagnosis，2009，29(2)：205-208.

[14] 程凱.基于小波分解和Duffing振子的變尺度微弱信號檢測[J].計算機(jī)測量與控制，2014，22(6)：1 732-1 734.

CHENG Kai.Weak signal detection based on scale transformation of doffing oscillator and wavelet decomposition[J].Computer Measurement & Control，2014，22(6)：1 732-1 734.

[15] 范留明.地震動信號的小波分析[J].物探化探計算技術(shù)，2000，22(1)：1-4.

FAN Liu-ming.Computing techniques for geophysical and geochemical exploration[J].Computing Techniques for Geophysical and Geochemical Exploration，2000，22(1)：1-4.

[16] 王大凱，彭進(jìn)業(yè).小波分析及其在信號處理中的應(yīng)用[M].北京：電子工業(yè)出版社，2006.

WANG Da-kai，PENG Jin-ye.The wavelet analysis and its application in signal processing[M].Beijing：Electronic Industry Press，2006.

[17] 馬建國，郭遼原，任國富.套管井電纜地層測試新技術(shù)[J].測井技術(shù)，2003，27(2)：95-98.

MA Jian-guo，GUO Liao-yuan，REN Guo-fu.New methods for cased hole wireline formation test[J].Well Logging Technology，2003，27(2)：95-98.

[18] 張家富，畢全福，胡廣軍，等.利用試井分析方法判斷MDT壓力測試的準(zhǔn)確性[J].新疆石油地質(zhì)，2013，34(6)：672-675.

ZHANG Jia-fu，BI Quan-fu，HU Guang-jun，et al.Judgment of accuracy of MDT pressure test using the method of well test analysis[J].Xinjiang Petroleum Geology，2013，34(6)：672-675.