1.一種從陣列聲波測(cè)井資料中提取反射波的方法,其特征在于:根據(jù)陣列聲波測(cè)井資料提供的全波列數(shù)據(jù),提取出直達(dá)波信號(hào)和反射波信號(hào),結(jié)合縱波速度、橫波速度及斯通利波速度及儀器結(jié)構(gòu)參數(shù),對(duì)提取的直達(dá)波信號(hào)再次進(jìn)行反射波提取,得到平行井軸的反射體產(chǎn)生的反射波,合成后得到完整的反射波。
2.根據(jù)權(quán)利要求1所述的從陣列聲波測(cè)井資料中提取反射波的方法,其特征在于:所述方法具體包括如下步驟:
a、以陣列聲波原始聲波波形數(shù)據(jù)作為輸入,采用中值濾波進(jìn)行反射波提取,令Wtol為原始波形數(shù)據(jù),Wzg處理得到的直達(dá)波波形數(shù)據(jù),Wf1為處理得到的反射波波形數(shù)據(jù),具體計(jì)算方法如下:
Wf1(i)=Wtol(i)-Wzg(i)
上式中,Ndep為深度行數(shù),median為中值濾波器,L為濾波器取中值的跨度,Wzg(i)為第i個(gè)深度直達(dá)波形數(shù)據(jù),為第i個(gè)深度下的中值濾波器中第j個(gè)波的原始波形數(shù)據(jù),Wf1(i)為第i個(gè)深度反射波形數(shù)據(jù);
b、對(duì)陣列聲波測(cè)井?dāng)?shù)據(jù)進(jìn)行處理,得到斯通利波速度Vst,橫波速度Vs及縱波速度Vp;
c、利用步驟a中處理得到的直達(dá)波波形數(shù)據(jù)Wzg,結(jié)合步驟b得到的波速及儀器結(jié)構(gòu)參數(shù),分步濾除斯通利波Wst、直達(dá)橫波Ws、直達(dá)縱波Wp,最終得到反映平行井軸的反射體的反射波Wf2;
d、對(duì)提取的兩個(gè)反射波Wf1、Wf2進(jìn)行合成,得到完整的反射波Wf:
Wf=Wf1+Wf2。
3.根據(jù)權(quán)利要求2所述的從陣列聲波測(cè)井資料中提取反射波的方法,其特征在于:所述濾除斯通利波具體為:
ΔTstij=(dj-di)/Vst
其中,為第i個(gè)深度t時(shí)刻的斯通利波直達(dá)波,為第i個(gè)深度下第j個(gè)接收器按△Tstij移動(dòng)后的波形,△Tstij為第i個(gè)深度下第j個(gè)接收器的斯通利波時(shí)差,為第i個(gè)深度濾除斯通利波直達(dá)波后的波形,median為中值濾波器,N為接陣列聲波儀器接收器個(gè)數(shù),d為接收器距聲源的距離,△T為不同接收器接收到直達(dá)波的時(shí)差。
4.根據(jù)權(quán)利要求2所述的從陣列聲波測(cè)井資料中提取反射波的方法,其特征在于:所述濾除直達(dá)橫波具體為:
ΔTsij=(dj-di)/Vs
其中,為第i個(gè)深度t時(shí)刻的橫波直達(dá)波,為第i個(gè)深度下第j個(gè)接收器按△Tsij移動(dòng)后的波形,△Tsij為第i個(gè)深度下第j個(gè)接收器的橫波時(shí)差,為第i個(gè)深度濾橫波直達(dá)波后的波形,median為中值濾波器,N為接陣列聲波儀器接收器個(gè)數(shù),d為接收器距聲源的距離,△T為不同接收器接收到直達(dá)波的時(shí)差。
5.根據(jù)權(quán)利要求2所述的從陣列聲波測(cè)井資料中提取反射波的方法,其特征在于:所述濾除直達(dá)縱波具體為:
ΔTpij=(dj-di)/Vp
其中,為第i個(gè)深度t時(shí)刻的橫波直達(dá)波,為第i個(gè)深度下第j個(gè)接收器按△Tpij移動(dòng)后的波形,△Tpij為第i個(gè)深度下第j個(gè)接收器的縱波時(shí)差,為第i個(gè)深度濾橫波直達(dá)波后的波形,median為中值濾波器,N為接陣列聲波儀器接收器個(gè)數(shù),d為接收器距聲源的距離,△T為不同接收器接收到直達(dá)波的時(shí)差。