專利名稱：用地震表面波的波形評(píng)估土壤性質(zhì)的制作方法
技術(shù)領(lǐng)域：
本申請(qǐng)要求2008年8月11日提交的美國(guó)臨時(shí)申請(qǐng)61/087，933的權(quán)益。本發(fā)明大體涉及地球物理勘探(geophysical prospecting)領(lǐng)域，并且更特別地 涉及地震數(shù)據(jù)處理。具體地，本發(fā)明是一種用于反演地震表面波數(shù)據(jù)以獲得土壤的彈性性 質(zhì)(elastic properties)的方法，土壤的彈性性質(zhì)例如土壤剪切波(shear wave)速或剪 切模量或剪切波衰減。特別地，可以得到作為離地表的深度的函數(shù)的彈性性質(zhì)或深度分 布(cbpth profile)性質(zhì)。該方法還能夠被用來評(píng)估土壤內(nèi)異常的位置，例如洞或埋入物 (buried objects)0
背景技術(shù)：
地震表面波還被稱為地滾波(ground-roll)或瑞雷波(Rayleigh waves)或勒夫 (Love)波，地震表面波局限于接近地球表面的區(qū)域，并且因此它們的傳播依賴于近地表彈 性性質(zhì)，特別是作為離地表的深度的函數(shù)的剪切波速。當(dāng)剪切波速由土壤剪切模量給出時(shí)， 剪切波速與土壤的剛度直接相關(guān)。(剪切波速是剪切模量的平方根除以密度。)土壤的剪 切波速或剪切模量分布或其它彈性特征能夠被直接用于工程或其它目的，或能夠被直接用 來改善土壤以下或近地表區(qū)域的地球物理勘探。除了剪切模量或波速，隨距離的彈性衰減 或振幅衰減也是能夠被用于工程和地球物理勘探的有用的信息。能夠通過反演(invert)地滾波頻散曲線(相速度-頻率)來確定土壤剪切波速 或模量，以獲得地面的剪切波速分布。由于壓實(shí)效應(yīng)(compaction effects)，速度在接近 地球表面處通常較慢并且隨深度而增加。表面波的較高頻率分量在接近地表處是受限的， 并且采樣(sample)較慢的土壤層。另一方面，較低的頻率分量采樣較深、較快的層。因此， 表面波的速度隨頻率變化，即，它是頻散(dispersive)的。特別地，速度隨頻率的增加而減 少。作為頻率的函數(shù)的頻散曲線(dispersion curve)的形狀能夠與為分層的速度分布計(jì) 算的頻散曲線比較，然后分布性質(zhì)(即層厚和剪切模量)能夠被更新以更好地匹配被測(cè)量 的頻散曲線。分層產(chǎn)生共振效應(yīng)和地滾波的不同模態(tài)的捕獲(trapping)。能夠通過反演地 滾波的基階模態(tài)(fundamental order mode)禾口高階模態(tài)(higher-order modes)的頻散曲 線來獲得更高的精度。表面波被轉(zhuǎn)換成土壤剪切波速或剪切模量、剪切波衰減或其它性質(zhì)的深度分布之 后，信息能夠被直接用作設(shè)計(jì)結(jié)構(gòu)(如建筑物和橋梁)的重要的工程參數(shù)。直接使用表面波 來描述近地表特征的其它應(yīng)用包括評(píng)估地震場(chǎng)地響應(yīng)、土壤壓實(shí)控制、繪制淺表(mapping the shallow surface)、評(píng)估地下材料的強(qiáng)度、路面評(píng)價(jià)、找出埋入的文化特征或異常、評(píng) 價(jià)下水道周圍的洞以及找出到基巖的深度。近地表速度分布能夠被間接用來改善較深地下 區(qū)域的物理結(jié)構(gòu)或物理性質(zhì)的確定，以評(píng)估或提取碳?xì)浠衔?hydrocarbon)。因?yàn)榻乇?具有低的速度并且不均勻，因此這對(duì)通過近地表的地震波產(chǎn)生了較大的影響，并且能夠限 制確定深區(qū)域的結(jié)構(gòu)和性質(zhì)的能力。近地表速度分布能夠被用來給來自較深的區(qū)域的地震 反射進(jìn)行時(shí)間校正或靜態(tài)校正，或者速度信息能夠被用于成像、遷移或反演地震數(shù)據(jù)。
利用表面波、瑞雷波或勒夫波的所有現(xiàn)有技術(shù)方法的問題在于難以求解或難以區(qū) 分不同表面波模態(tài)的不同頻散曲線。一個(gè)內(nèi)在的困難產(chǎn)生于頻率和速度的不確定性，并且 當(dāng)錯(cuò)誤模態(tài)被識(shí)別或被不正確地選擇時(shí)，反演受到損害。第二個(gè)困難產(chǎn)生于缺乏對(duì)源相位 的了解和2π以上的相位變化的不確定性。此外，所有這些方法在分析中所用的接收器擴(kuò) 展距離(receiver spread distance)上對(duì)表面波性質(zhì)求平均，并且這種求平均限制了分辨 率(resolution)。最后，模態(tài)與噪聲之間的干涉以及表面波的衰減會(huì)使地震振幅變形，使識(shí) 別各個(gè)的頻散曲線變得困難。接下來，更詳細(xì)地論述確定土壤剪切模量或剪切波速的傳統(tǒng) 方法。使用表面波來表征土壤剪切模量或波速的現(xiàn)有方法涉及獲取地震數(shù)據(jù)，其后有兩 個(gè)處理階段(1)測(cè)量作為頻率的函數(shù)的頻散曲線和(2)之后反演頻散曲線以獲得作為深 度的函數(shù)的剪切模量。相似的方法可以被分別用來得到除了剪切模量之外的其它性質(zhì)，例 如衰減性質(zhì)。測(cè)量頻散曲線的方法在源的數(shù)量和獲取的接收器的數(shù)量上差別很大。最早的 方法使用了一個(gè)源和單對(duì)接收器。較新的方法使用一個(gè)源和以規(guī)則間隔隔開的多個(gè)接收器 (大約20個(gè)或更多個(gè))。多數(shù)方法使用產(chǎn)生瑞雷波的壓縮源(compression source)，但是 剪切源也能夠被用來產(chǎn)生勒夫波。使用任一種波的類型的方法是相同的。階段1 用于一對(duì)接收器的方法用瑞雷波來表征土壤剪切模量的最早應(yīng)用是“穩(wěn)態(tài)法”，其涉及使用可控震源 (seismic vibrator)以單一頻率或以緩慢變化的頻率來振動(dòng)地面(Bahjat的美國(guó)專利 3，864，667 (1975))。在每個(gè)頻率，測(cè)量?jī)蓚€(gè)檢波器(geophone)的響應(yīng)之間的相位差。根據(jù) 這些測(cè)量值推斷接收器之間的近地表的性質(zhì)。然而，這些測(cè)量值花費(fèi)長(zhǎng)時(shí)間才能得到。在 20世紀(jì)80年代，發(fā)展了 SASW(表面波波譜分析)方法(Nazarian等人，“Use of spectral analysis of surface waves method for determination of moduli and thickness of pavement system(使用表面波波譜分析方法確定路面系統(tǒng)的模量和厚度)”，Transport. Res. Record 930，38_45 (1983))。該方法通過首先計(jì)算由兩個(gè)傳感器記錄的信號(hào)之間的互 功率譜(cross power spectrum)并且之后展開相位來確定頻散曲線。SASW方法和穩(wěn)態(tài)法都遭遇到類似的問題。因?yàn)橐淮蝺H使用一對(duì)接收器，所以難以 區(qū)分不同表面波模態(tài)的影響與任何被記錄的噪聲的影響。接收器之間的距離和一對(duì)接收器 與源之間的距離被改變，以將高階模態(tài)的效應(yīng)最小化，但沒有消除。這些測(cè)量的內(nèi)在問題是 同相位的2π的倍數(shù)周圍的模糊。接收器必須足夠靠近以使得相位不會(huì)變化超過2π。能 夠區(qū)分△的相位變化與Δ+2π的變化是很重要的；后者對(duì)應(yīng)于較慢的速度。有時(shí)，少數(shù)附 加接收器被用來輔助相位展開。近來，小波變換被用作能夠改善單個(gè)地滾波模態(tài)的分離并將來自其它模態(tài)的貢獻(xiàn) 最小化的一禾中方法，如 Holschneider 等人在"Characterization of dispersive surface waves using continuous wavelet transforms (用連續(xù)小波變換描述表面波頻散特征)，，、 Geophys. J. Int. 163，463-478 (2005)中說明的。小波變換是能夠更好地對(duì)單獨(dú)模態(tài)定位的 時(shí)頻變換。該方法減小了多個(gè)模態(tài)之間的干涉，但是因?yàn)橥辔坏? π模糊，噪聲仍然是問 題。Holschneider等人建立了小波變換域內(nèi)的表面波傳播效應(yīng)的數(shù)學(xué)模型，并使用該數(shù)學(xué) 模型來首先求解一次一種模態(tài)的頻散曲線和之后的衰減曲線。階段1:多個(gè)接收器
堪薩斯州地質(zhì)調(diào)查局發(fā)展了表面波多通道分析方法(MASW)。在此方法中配 置了 20-65個(gè)或更多個(gè)接收器并記錄來自一個(gè)脈沖源或振動(dòng)源的數(shù)據(jù)。(Park等人， "Multichannel analysis of surface waves (表面波多通道分析)，，，Geophysics 64, 800-808(1999), VX R Park φ A, "Multichannel analysis of surface waves (MASff) active and passive methods (表面波多通道分析主動(dòng)禾口被動(dòng)方法)，，，The LeadingEdge 26,60-64(2007))。數(shù)據(jù)被處理為單炮集(single shot gather)，振幅被標(biāo)準(zhǔn)化或被縮放 (scaled)，并且之后被轉(zhuǎn)換到頻率波數(shù)(f_k)域或頻率慢度(frequency slowness) (f-p) 域。之后，在f_k或f-p域內(nèi)的最大振幅點(diǎn)處，一種或更多種模態(tài)的頻散曲線被選擇。一種 或更多種模態(tài)的頻散曲線被用在隨后的近地表性質(zhì)的反演中(Beaty等人，‘‘Itepeatability of multimode Rayleigh-wave dispersion studies (多模態(tài) 卷雷波頻散可重復(fù)性研究)”， Geophyics 68782-790(2003))?？梢灾貜?fù)獲取被射入新接收器擴(kuò)展(receiver spread)中 的每個(gè)新炮和被分別分析的每個(gè)炮集，以獲得每個(gè)擴(kuò)展的ID近地表速度分布。之后每個(gè)ID 分布與插值格式(interpolation scheme)合并以產(chǎn)生近地表面的2D分布。多通道方法是對(duì)SASW方法的改進(jìn)。一個(gè)優(yōu)點(diǎn)在于使用較近間隔的接收器將相 位變化的2 π模糊最小化；然而，可能根據(jù)單集評(píng)估源相位(Hermann和Ammon，"Surface Waves, Receiver Functions, and Crustal Structure :Version 3. 3 (表面波、接收器功能 以及地殼構(gòu)造版本3. 3) ”，Computer Programs in Seismology，圣路易斯大學(xué)，httpwww. eas. slu. edu/People/RBHermann/CPS330. htl. (2004))。第二個(gè)優(yōu)點(diǎn)在于變換到 f_k 或 f-p 域本身涉及軌跡的求和或堆積，這改善了頻率分辨率并減少了噪聲問題。然而，代價(jià)是橫向 分辨率的損失；可能檢測(cè)擴(kuò)展的寬度內(nèi)的速度變化。在兩個(gè)通道的情況下，選擇偏移量范 圍以強(qiáng)調(diào)或不強(qiáng)調(diào)各種模態(tài)也是很重要的(Xia等人，“Utilization of high-frequency Rayleigh waves in near-surface geophysics (在近地表地球物理學(xué)中高頻瑞雷波的利 用)”，The Leading Edge 23，753-759 (2004))。此外，接收器必須是均勻間隔的緊挨著， 以便不混淆(alias)地滾波。此外因?yàn)榈貪L波振幅被嚴(yán)重衰減，它們逐個(gè)軌跡變化，并且 總和使振幅失真。使用了歸一化(normalization)或振幅平衡，但仍然難以選擇并區(qū)分多 種干擾的地滾波模態(tài)。Lefebvre和Benhassen(美國(guó)專利申請(qǐng)公開No. 2005/0143924A1) 教導(dǎo)了用小波變換來提高區(qū)分不同模態(tài)的能力。Forbriger( "Inversion of shallow-seismic wavefields :I. Wavefield transformation ( $ M Jft ■ M 白勺 β M 1.波場(chǎng)變換)”，Geophys. J. Int. 153，719-734(2003))說明了在選擇多模態(tài)頻散曲線中的 困難和當(dāng)頻散曲線被不正確地選擇或被誤識(shí)別時(shí)后續(xù)反演的問題。Forbriger和Ryden、 Park( "Fast simulated annealing inversion of surface waves on pavement using phase-velocity spectra (用相速度譜在路面上快速模擬退火反演表面波)，，，Geophysics 71，R49-R58(2006))都避免選擇頻散曲線而是直接反演f_p變換結(jié)果。這些方法涉及f_p 域中的表面波傳播的數(shù)學(xué)模型，其涉及幾個(gè)假設(shè)和近似值(approximation)。階段1 多個(gè)接收器和多個(gè)源Ernst 等人在 “Tomography of dispersive media (頻散介質(zhì)的層析成像)”(J. Acoust. Soc. Am. 108,105-115 (2000))禾口 "Removal ofscattered guided waves from seismic data (從地震數(shù)據(jù)去除分散的導(dǎo)波)"(Geophysics 67，1240-1248 (2002))中討論 了同時(shí)用多個(gè)源位置和多個(gè)接收器來獲得橫向變化的相速度曲線。本申請(qǐng)不是關(guān)于描述近地表特征的，而是關(guān)于用于石油勘探地震應(yīng)用的分散的地滾波的緩解。他們的過程包括 操作的級(jí)聯(lián)(cascade)。它們首先基于一般行程時(shí)間使用斷層層析成像方法(tomography method)將橫向變化的相速度反演為頻率的函數(shù)。它們假設(shè)橫向變化較小并且假設(shè)它們能 夠得到時(shí)間窗口中一種模態(tài)的分離。這后一假設(shè)要求源和接收器分開得充分遠(yuǎn)從而模態(tài)能 夠及時(shí)較好地被分離，但是由于表面波模態(tài)的強(qiáng)烈衰減，很難或不可能得到近地表特征描 述所需的數(shù)據(jù)。廣義行程時(shí)間的使用涉及數(shù)據(jù)的相位的導(dǎo)數(shù)的計(jì)算，并且增加了在地滾波 相速度的確定中相位模糊的問題。階段1 衰減與頻率曲線很多現(xiàn)有技術(shù)集中在相速度頻散曲線，Xia等人(同上)討論了作為頻率的函 數(shù)的衰減曲線的使用。作為深度與剪切模量的函數(shù)的品質(zhì)因子(quality factor) (Q)也 是重要的工程量，但是衰減數(shù)據(jù)的反演具有較差的穩(wěn)定性。一般，假設(shè)衰減不依賴于頻 率(Ernst 等人，“Removal of scattered guided waves from seismic data(從地震 數(shù)據(jù)去除分散的導(dǎo)波)”，Geophysics 67，1240-1248 (2002)；和 Kulesh 等人，“Modeling of Wave Dispersion Using Continuous Wavelet Transforms II :Wavelet Based Frequency-Velocity Analysis (用連續(xù)小波變換對(duì)波的頻散建模II 基于頻率速度分析 的小波)”，Pure & Applied Geophysicsl65，255-270 (2008))。然而，對(duì)于用來描述近地表 特征的表面波，這個(gè)假設(shè)是受限的。因?yàn)樗p通常減少了土壤深度的增加，所以，作為頻率 的函數(shù)，表面波的衰減還必須以與相速度分散的相同方式減少。階段2 反演近地表速度分布的頻散曲線有許多算法可用來根據(jù)頻散曲線求解近地表速度分布，但是所有這些方法的成 功依賴于輸入頻散曲線的精度。反演是非線性模型優(yōu)化問題，其中模型是近地表速度分 布。參數(shù)包括層深和層剪切模量。算法包括線性最小二乘反演法、Levenberg Marquardt 法、擬牛頓法以及近來的模擬退火法(Beaty等人，“Simulated annealing inversion of multimode Rayleigh wave dispersion curves for geological structure (地質(zhì)構(gòu)造 的多模態(tài)瑞雷波頻散曲線的模擬退火反演)”，Geophys. J. Int. 151,622-631 (2002))?？?用軟件包括免費(fèi)軟件(Hermann 禾口 Ammon，“Surface Waves, Receiver Functions，and Crustal Structure !Version 3. 3 (表面波、接收器函數(shù)以及地殼結(jié)構(gòu):3· 3版)”，Computer Programs in Seismology,圣路易斯大學(xué)，http:www. eas. slu. edu/Peop 1 e/RBHermann/ CPS330. htl. (2004))和商業(yè)軟件(SeisOpt ReMi ，http:www. optimsoftware. com;和 KansasGeological Survey :http:/www. kgs. ku. edu/sofiware/surfseis. /index, html.)0保留對(duì)改進(jìn)的方法的需求，該改進(jìn)方法用于通過獲得到近地表面特征反演的輸入 的高分辨率、橫向變化、多模態(tài)頻散和衰減曲線，來將表面波變換成近地表面性質(zhì)的深度分 布。特別地，該方法應(yīng)該將源相位中的模糊最小化。本發(fā)明滿足了此需求。

發(fā)明內(nèi)容
在一個(gè)實(shí)施例中，本發(fā)明是一種方法，該方法用于根據(jù)地震數(shù)據(jù)軌跡(seismic data trace)來評(píng)估地球的近地表區(qū)域的剪切模量、剪切波速、剪切波衰減或另一個(gè)物理性 質(zhì)，地震數(shù)據(jù)軌跡對(duì)應(yīng)于地震勘測(cè)區(qū)域中的至少一個(gè)源位置和多個(gè)接收器位置，該方法包 括
(a)將近地表區(qū)域分成一個(gè)或更多個(gè)單元；(b)對(duì)于每個(gè)軌跡和相應(yīng)的地震源和接收器位置，計(jì)算從源位置到接收器位置穿 過每個(gè)中間單元的射線路徑距離；(c)對(duì)于每個(gè)軌跡，同時(shí)求解至少兩個(gè)地表一致性分量，每個(gè)地表一致性分量表征 穿過單元的地震表面波的傳播的濾波效應(yīng)或相應(yīng)的地震源或接收器與地球的耦合，其中所 述解通過比較預(yù)測(cè)的表面波和來自勘測(cè)的相應(yīng)數(shù)據(jù)軌跡來迭代地優(yōu)化分量，所述預(yù)測(cè)的表 面波通過采用假設(shè)的或迭代更新的分量和來自(b)的射線路徑信息來計(jì)算；(d)選擇一個(gè)或更多個(gè)地表一致性分量，并使用它們通過數(shù)值反演(numerical inversion)來計(jì)算剪切模量或近地表區(qū)域的另一個(gè)性質(zhì)。


通過參考下面詳細(xì)描述和附圖，將更好地理解本發(fā)明及其優(yōu)點(diǎn)，其中圖1A-1B圖示了對(duì)于線性應(yīng)用(IA)和面積應(yīng)用(areal application) (IB)用多 個(gè)接收器獲得地震數(shù)據(jù)；圖2A-2B圖示了對(duì)于線性應(yīng)用(2A)和面積應(yīng)用(2B)用多個(gè)源獲得地震數(shù)據(jù)；圖3A-3B圖示了對(duì)于線性應(yīng)用(3A)和面積應(yīng)用(3B)用跨過地表區(qū)域的具有不同 偏移量的多個(gè)源_接收器對(duì)獲得地震數(shù)據(jù)；圖4是物理過程的圖示，當(dāng)表面波沿地球的表面從源傳播到接收器時(shí)，物理過程 影響波形；圖5是地圖中地震勘測(cè)的一小部分的圖示，其示出了源和接收器的位置以及區(qū)域 或單元的劃分，并示出了穿過不同區(qū)域從源到接收器的示例射線路徑；圖6是本發(fā)明的方法的一個(gè)實(shí)施例中基本步驟的流程圖；圖7是本發(fā)明的一個(gè)實(shí)施例中用來執(zhí)行圖6的模型優(yōu)化部分的步驟的流程圖，其 中產(chǎn)生調(diào)整表面波傳播的地表一致性分量(傳遞函數(shù))；圖8示出了用于半空間上單個(gè)薄地面層的表面波的完整波動(dòng)方程計(jì)算機(jī)仿真的 四個(gè)源-接收器距離(偏移量)的四條軌跡；圖9A-9D示出了圖8的示例的模型優(yōu)化結(jié)果，其中波形變化已經(jīng)被分解成作為頻 率的函數(shù)的復(fù)合(振幅和相位)源和傳播分量；圖10示出了根據(jù)圖9A-9D的分量計(jì)算的預(yù)測(cè)的波形輸出；該波形可以與圖8的地 震數(shù)據(jù)輸入比較；圖11A-11B示出了用于半空間上20個(gè)薄地面層的表面波的完整波動(dòng)方程計(jì)算機(jī) 仿真的軌跡數(shù)據(jù)(IlA)和f-k譜(IlB)，其中可以看到多達(dá)6種地滾波模態(tài)的干涉；圖12示出了與圖IlA的軌跡數(shù)據(jù)有少數(shù)偏移量(軌跡到接收器的距離)的被選 擇的波形；波形沒有簡(jiǎn)單的形式，但是由多種模態(tài)的疊加組成；圖13示出了與圖12中所示的相同偏移量的預(yù)測(cè)的波形輸出；在此預(yù)測(cè)中使用了 6種模態(tài)的參數(shù)；圖14A-14F示出了來自圖12的一個(gè)數(shù)據(jù)軌跡(14A)和預(yù)測(cè)，其中預(yù)測(cè)的波形輸出 被限制為一種(14B)、兩種(14C)、三種(14D)、四種(14E)和六種(14F)模態(tài)；包含更多種 模態(tài)改善了數(shù)據(jù)軌跡(14A)的匹配；以及
圖15顯示了用本發(fā)明的方法計(jì)算的圖11A-11B的6種不同的地滾波模態(tài)的速度 頻散曲線(速度作為頻率的函數(shù))。本發(fā)明連同示例性實(shí)施方式被描述。然而，在一定程度上，下面的描述針對(duì)本發(fā)明 的具體實(shí)施例或具體使用，這意欲僅作為示例，并不被作為本發(fā)明范圍的限制。相反，其意 欲覆蓋如通過所附權(quán)利要求限定的可以被包括在本發(fā)明的范圍以內(nèi)的所有替換、修改和等 價(jià)物。
具體實(shí)施例方式本發(fā)明是從地震數(shù)據(jù)獲得作為頻率的函數(shù)的表面波的速度頻散曲線和衰減曲線 的方法。優(yōu)選地，用多個(gè)源和多個(gè)接收器來記錄數(shù)據(jù)。之后，曲線或它們的修改能夠被用來 通過已知的數(shù)值技術(shù)反演近地表性質(zhì)，例如作為深度的函數(shù)的剪切波速或剪切模量或剪切 波衰減。模型優(yōu)化和數(shù)據(jù)冗余被用來求解地表一致性模型參數(shù)，該地表一致性模型參數(shù)最 好表示當(dāng)表面波沿地表從源傳播到接收器時(shí)被記錄的地震數(shù)據(jù)波形和波形的變化。參數(shù)涉 及用于每個(gè)源位置、每個(gè)接收器位置以及通過地表的每個(gè)區(qū)域的傳播的各個(gè)分量或?yàn)V波器 傳遞函數(shù)(filter transfer function)。模型參數(shù)可以包括表面波的多種模態(tài)和地表一致 性橫向變化。因此，本發(fā)明將地震數(shù)據(jù)與兩級(jí)反演方法一起使用，地震數(shù)據(jù)優(yōu)選用多個(gè)源和多 個(gè)接收器記錄。在階段1，表面波波形的變化優(yōu)選被分解成關(guān)于每個(gè)源、每個(gè)接收器以及地 表的每個(gè)小區(qū)域的地表一致性傳遞函數(shù)。隨后，在階段2，每個(gè)區(qū)域的傳遞函數(shù)被反演以確 定作為深度的函數(shù)的土壤性質(zhì)或近地表面性質(zhì)(例如剪切波速)。該方法能夠求解具有垂 直變化性質(zhì)和橫向變化性質(zhì)的介質(zhì)的表面波的復(fù)合多模態(tài)屬性?？梢员苊鈧鹘y(tǒng)方法的求解 (resolution)中地滾波模態(tài)的錯(cuò)誤識(shí)別或假設(shè)橫向均勻的土壤性質(zhì)的錯(cuò)誤和限制。地震數(shù)據(jù)獲取首先，必須記錄或獲得地震數(shù)據(jù)。出于其它目的(例如對(duì)地表成像)被記錄的地震 數(shù)據(jù)通常能夠被使用，見下述優(yōu)選的獲取參數(shù)的論述。如果這些數(shù)據(jù)不可用，則應(yīng)該得到用 于獲得表面波性質(zhì)的專門勘測(cè)。如果僅需要沿地表線的近地表信息，則能夠用沿該線的源 和接收器記錄數(shù)據(jù)作為標(biāo)準(zhǔn)2D地震記錄。如果需要面積信息，則源和接收器應(yīng)該在區(qū)域展 開，源和接收器都跨過被調(diào)查的區(qū)域。如同任何地震獲取，源和接收器位置必須被勘測(cè)并且 必須可與被記錄的數(shù)據(jù)軌跡一起使用。能夠使用任何類型的地震源，例如甘油炸藥、振子、 氣槍、重物、錘擊、地震火炮(seismic guns)等。如在圖1A-1B中所示的線性獲取(圖1A) 和面積獲取(圖1B)，優(yōu)選地，源應(yīng)該被射入多個(gè)接收器中(例如十個(gè)或更多個(gè))。在圖IA 中，其顯示了一個(gè)源(10)和許多個(gè)接收器(11)。可以理解在源與每個(gè)接收器之間有直的 射線路徑(12)。在圖1中僅畫出了可能的射線路徑中的少數(shù)幾條。源與接收器之間的距 離被稱作偏移量。相反地，如圖2A和圖2B所示，當(dāng)源隨時(shí)間被移動(dòng)到不同的位置時(shí)，每個(gè) 接收器(20)應(yīng)該記錄來自具有不同射線路徑(22)的多個(gè)(優(yōu)選地十個(gè)或更多個(gè))炮/射 點(diǎn)(shots) (21)的數(shù)據(jù)。通過布置許多個(gè)接收器(60-1000)的接收器擴(kuò)展并且之后一個(gè)接 一個(gè)地爆破源形成整體擴(kuò)展，很容易地實(shí)現(xiàn)這兩個(gè)要求。炮/射點(diǎn)的數(shù)量不需要與接收器 的數(shù)量相同。為了減少成本，炮的數(shù)量可以比接收器的數(shù)量少很多，或反之亦然。如果可用 的儀器僅能夠記錄有限數(shù)量的通道，則記錄能夠被限制為在距離每個(gè)源的最大(最大偏移量)距離以內(nèi)的接收器。優(yōu)選的是布置源和接收器以使得多個(gè)源_接收器射線路徑(優(yōu)選地是10-20條或 更多條)通過用在特定的最小偏移量到特定的最大偏移量范圍內(nèi)的不同偏移量探測(cè)的每 個(gè)區(qū)域。在圖3A中，通過方框30示出了一個(gè)區(qū)域。附圖示出了穿過方框30從源(例如 31)到接收器(例如32)的一些射線路徑(33)。有很多具有各種偏移量的這些射線路徑 (未全顯示)。圖3B圖示了在面積勘測(cè)(areal survey)中的相同的特征。最小偏移量應(yīng) 該足夠長(zhǎng)以形成平面波傳播，平面波傳播在接近表面波波長(zhǎng)的1/2的距離處建立。這個(gè)距 離近似等于穿透深度(cbpth of penetration)。最大偏移量?jī)?yōu)選足夠短以避免表面波的過 度衰減。通常它是最小偏移量的2-3倍。能夠記錄更短的或更長(zhǎng)的偏移距離。接收器間隔 應(yīng)該優(yōu)選地不大于期望的橫向分辨率的大小，但是如果可能的話具有一半采樣可能是有用 的。源和接收器可以被布置在一個(gè)矩形格上，或可以使用非矩形格或偽隨機(jī)采樣?？赡苡?必要設(shè)計(jì)勘測(cè)(survey)以利用道路來放置源并避免結(jié)構(gòu)。還可能有用的是將源和接收器 放置在被覆蓋的期望的區(qū)域的外側(cè)以保證在邊緣處有充分的數(shù)據(jù)冗余。為了優(yōu)化勘測(cè)參數(shù)(survey parameters)，例如最小偏移量、最大偏移量以及接收 器和源采樣，有用的是具有關(guān)于面積的信息。進(jìn)行一些初步測(cè)試以確定源的頻率范圍和地 滾波的速度范圍可能是有用的。這將能夠計(jì)算表面波波長(zhǎng)。還可能有用的是在區(qū)域的中間 用好的接收器采樣將炮記錄到2-D線，以使得f_k或f-p分析能夠被用來導(dǎo)出初始參數(shù)。對(duì) 于這個(gè)特殊的線，接收器采樣(例如，間隔)應(yīng)該小于最大頻率乘以最慢的速度從而數(shù)據(jù)不 會(huì)被混淆。階段1 模型優(yōu)化除了射線路徑的不同，區(qū)分表面波與其它地震波(例如從界面反射的下行波)的 一個(gè)特征是它們的波形。表面波的波形是高振幅、低頻率和振蕩的。另一方面，震源子波 (source wavelet)是典型的高頻脈沖或短時(shí)子波。圖4例示了波形的這種變化的物理過 程，圖4示出了一塊層狀地表110。當(dāng)表面波的一種模態(tài)沿地表(射線路徑103)從位于表 面位置101處的源傳播到位于表面位置102處的接收器時(shí)，它經(jīng)歷了相當(dāng)大的大地濾波，從 源(104)的短時(shí)脈沖變化到表面波(105)的更大振蕩的低頻波形。表面波被延遲并且被衰 減。傳播(103)的距離越遠(yuǎn)，變化越大。給定源波形(104)的評(píng)估和被記錄的輸出(105)， 能夠通過計(jì)算每個(gè)軌跡的傳遞函數(shù)來量化這個(gè)變化。這可以被稱作總軌跡傳遞函數(shù)，并且 它以總大地濾波為特征，大地濾波導(dǎo)致波形從在源處產(chǎn)生的波形變化到在地震檢波器記錄 的波形。濾波器通常以它們的傳遞函數(shù)或它們的脈沖響應(yīng)為特征。傳遞函數(shù)被定義為在頻 域內(nèi)濾波器的輸出除以濾波器的輸入對(duì)于圖4中的地震采集，傳遞函數(shù)T(f)為T(f) = ^m = ^l(1)
input(f) S(f)其中D(f)是數(shù)據(jù)軌跡105的傅立葉變換，S (f)是震源子波104的傅立葉變換。同 樣地，傳遞函數(shù)能夠通過子波與數(shù)據(jù)軌跡的互相關(guān)得到，并通過震源子波的自相關(guān)劃分。W/姆)(2)
S\f)S(f)其中星號(hào)表示復(fù)共軛。傳遞函數(shù)是復(fù)合的；它必須包括振幅和相位或?qū)嵅亢吞摬恳酝暾乇碚鳛V波器效應(yīng)。計(jì)算傅立葉逆變換產(chǎn)生濾波器的脈沖響應(yīng)。傳遞函數(shù)和脈沖響 應(yīng)分別等價(jià)描述了頻域和時(shí)域中的特征。眾所周知，獨(dú)立的線性濾波器能夠被一起卷積(或在頻域內(nèi)一起相乘)以產(chǎn)生組 合的濾波器。從圖4可見，表面波的波形的整體效果由不同的物理過程組成，每個(gè)物理過程 與地表位置或區(qū)域相關(guān)聯(lián)。例如，在位置101處源能量與地滾波模態(tài)的耦合，穿過區(qū)域106 的傳播，穿過區(qū)域107的傳播，以及在位置102處接收器與地面的耦合。此外，對(duì)穿過區(qū)域 106的傳播的濾波應(yīng)該是穿過區(qū)域106的傳播距離(108)的函數(shù)，并且對(duì)穿過區(qū)域107的傳 播的濾波應(yīng)該是穿過區(qū)域107的傳播距離(109)的函數(shù)。對(duì)傳播效應(yīng)的濾波包括延遲的效 應(yīng)或作為頻率(頻散)的函數(shù)的速度和作為頻率的函數(shù)的振幅衰減。因此，對(duì)于一種表面 波模態(tài)，總的軌跡傳遞函數(shù)T (f)能夠被分解成獨(dú)立的傳遞函數(shù)或上述物理過程中的每一 個(gè)的濾波器，或T (f) = T101 (f) T106 (f, d108) T107 (f, d109) T102 (f) (3)使用方程1和方程3，能夠得到表示被記錄的地滾波的波形的模型表達(dá)式M(f)M(f) =S (f) T101 (f) T106 (f, d108) T107 (f, d109) T102 (f) (4)方程式4是單地滾波模態(tài)的波形的模型。進(jìn)一步假設(shè)數(shù)據(jù)由多種模態(tài)的線性疊加 組成。并且，首先假設(shè)各個(gè)模態(tài)不相互作用或被耦合。因此對(duì)于N種模態(tài)，復(fù)合地滾波波形 的模型能夠表達(dá)為如下形式M(/) = 5(/)|； Tmj (Z)Tl06j (/, ^108 )TW1J (/, J109 )Tmj (/)(5)
7=1其中總和覆蓋了 j = 1到N的所有不同模態(tài)。模態(tài)的數(shù)量N通常較小，大約為2-6。在方程式5中，各個(gè)物理過程被表示為頻域中的傳遞函數(shù)或?yàn)V波器。優(yōu)選地包括 所有這些濾波器，即使僅需要傳播效應(yīng)來求解近地表性質(zhì)。傳遞函數(shù)是復(fù)合的，并包括作為 頻率的函數(shù)的振幅和相位或?qū)嵅亢吞摬?。?duì)于傳播，這包括速度和衰減效應(yīng)。許多不同的過 程能夠被包含在模型表達(dá)式中。一些是已知的或能夠被測(cè)量的。例如，一個(gè)濾波器是記錄 系統(tǒng)濾波器。本領(lǐng)域內(nèi)技術(shù)人員能夠測(cè)量記錄系統(tǒng)的脈沖響應(yīng)。另一個(gè)可能是能夠從制造 商得到的傳感器的儀器響應(yīng)。第三個(gè)是振幅擴(kuò)展函數(shù)(amplitude spreading function), 振幅擴(kuò)展函數(shù)與每個(gè)軌跡的源到接收器的總距離的平方根成反比。將幾個(gè)相關(guān)的過程歸 入一項(xiàng)可能是有用的。例如，如果有接收器陣列，則陣列效應(yīng)能夠與接收器耦合項(xiàng)被歸在 一起。在兩種情況中，它們與具體的接收器地表位置有關(guān)。還可能有用的是包括被假定 為已知的一個(gè)項(xiàng)，并且之后具有第二變量校正項(xiàng)。例如，平均震源子波(average source wavelet)能夠被評(píng)估，并且之后源耦合項(xiàng)能夠包括子波的耦合和逐炮變化(shot-to-shot variation) 0對(duì)于線性勘測(cè)或面積勘測(cè)，本發(fā)明使用具有不同的射線路徑和源_接收器距離的 數(shù)據(jù)軌跡的冗余，來將地滾波波形的總變化分解成獨(dú)立的地表一致性分量或表示獨(dú)立的物 理過程的濾波器。圖5圖示了面積數(shù)據(jù)的數(shù)據(jù)冗余的這種使用。圖5示出了地圖中可能的 地震勘測(cè)的一小部分。在所顯示的部分中，有一個(gè)源，其表面位置被標(biāo)記為201。源實(shí)際上 可以位于炮眼(shot hole)中的表面以下。在圖中，有24個(gè)這種源，每一個(gè)用太陽爆炸符 號(hào)(sun burst symbol)表示。還有用梯形符號(hào)標(biāo)記的49個(gè)接收器，49個(gè)接收器中的一個(gè) 被標(biāo)記為202。假設(shè)每次源被爆破就為每個(gè)接收器做記錄，則24X49或1176組數(shù)據(jù)軌跡是可用的。對(duì)于此發(fā)明，勘測(cè)面積被分成區(qū)域或單元。在這種情況下，虛線示出矩形單元，并 且附圖的中心的那一個(gè)被標(biāo)記為203。每個(gè)源和接收器都通過射線路徑連接，射線路徑首 次逼近(on first approximation)能夠被認(rèn)為是直線。僅示出了少數(shù)示例射線路徑，例如 204,205和206，以及這三個(gè)射線路徑各自以有不同行程采樣區(qū)域203。假設(shè)表面波僅一種 模態(tài)，則在頻域中每個(gè)頻率處數(shù)據(jù)軌跡能夠被用來求解地表一致性傳遞函數(shù)的未知分量。 在這個(gè)一種模態(tài)的示例中，未知量包括24個(gè)源耦合項(xiàng)、49個(gè)接收器耦合項(xiàng)以及9個(gè)傳播項(xiàng) (因?yàn)橛?個(gè)離散單元)，總共82項(xiàng)未知量。(這圖示了圖6的步驟303中的優(yōu)化，其將在 下面討論。)由于有1176塊數(shù)據(jù)和82項(xiàng)未知量，因此問題是超定的。如果有兩種模態(tài)，則 將有164項(xiàng)未知量。本發(fā)明的關(guān)鍵部分在于涉及不同物理過程的未知量是被同時(shí)確定的。 有具有關(guān)于201處的源的信息的49個(gè)軌跡，有具有關(guān)于接收器202的信息的24個(gè)軌跡，以 及具有穿過區(qū)域203的至少460條射線路徑，其中每一條射線路徑將通過穿過區(qū)域203的 不同路徑長(zhǎng)度來被賦權(quán)重。本領(lǐng)域內(nèi)有經(jīng)驗(yàn)的人員將把這個(gè)問題識(shí)別為地表的2D直射線 層析成像問題。本發(fā)明的方法優(yōu)選地使用帶有參數(shù)的模型表達(dá)式和價(jià)值函數(shù)，其中模型表達(dá)式表 征重要的物理過程的效應(yīng)，價(jià)值函數(shù)表示數(shù)據(jù)與模型表達(dá)式之間的擬合優(yōu)度。有用的是包 括盡可能多的已知效應(yīng)，例如記錄系統(tǒng)儀器響應(yīng)。一些分量對(duì)于具體的獲取可能不重要，并 且可以被忽略。例如，如果使用了單獨(dú)的接收器并且它們通常都較好地耦合到地面，則可以 忽略耦合或用相似的標(biāo)量代替。在最小值處參數(shù)總共應(yīng)該至少為兩個(gè)，例如穿過一個(gè)區(qū)域 的傳播項(xiàng)和將能量劃分為地滾波模態(tài)的源項(xiàng)。最小平方優(yōu)化(例如)的價(jià)值函數(shù)是每個(gè)軌 跡的實(shí)際數(shù)據(jù)與模型表達(dá)式之間的差值的平方在軌跡個(gè)數(shù)范圍上的總和?？偤透采w了在求 解參數(shù)中所用的所有數(shù)據(jù)軌跡。能夠在頻域內(nèi)、在時(shí)域內(nèi)或在一些其它的傳遞域內(nèi)表達(dá)價(jià) 值函數(shù)，例如包括f-k、氡(radon)、小波、Gabor、復(fù)數(shù)道(complex trace)、Hilbert以及技 術(shù)人員了解的其它域。濾波器能夠被用到數(shù)據(jù)和模型中，例如互相關(guān)濾波。此外，能夠用干 涉法處理來自主動(dòng)源或被動(dòng)源的地震記錄。優(yōu)選地，使用頻域，并且價(jià)值函數(shù)與實(shí)部表示和 虛部表示(即數(shù)據(jù)與模型表達(dá)式的振幅和相位表示)之間的差有關(guān)。使用實(shí)部和虛部的示 例價(jià)值函數(shù)(對(duì)于步驟303)為
權(quán)利要求
1.一種將對(duì)應(yīng)于地震勘測(cè)中地下區(qū)域的至少一個(gè)源位置和多個(gè)接收器位置的地震數(shù) 據(jù)軌跡轉(zhuǎn)換成彈性剪切波傳播速度或包括剪切模量的由此可導(dǎo)出的地球其它物理性質(zhì)的 近地表地球模型的方法，所述方法包括(a)將所述地下區(qū)域的近地表部分劃分成一個(gè)或更多個(gè)單元；(b)同時(shí)求解所述地震軌跡的地表一致性分量的至少兩種類型，其中地表一致性分量 表征濾波效應(yīng)，是穿過單元的地震表面波模態(tài)的傳播即傳播濾波，或由于在源位置的表面 波能量的產(chǎn)生即源濾波，或由于在接收器位置的源波能量的接收即接收器濾波，并且其中 地表一致性分量的“類型”是以下類型之一 (i) 一種表面波模態(tài)的傳播濾波，其在被任何 所述地震軌跡穿過的所有單元被評(píng)價(jià)；(ii)源濾波，其在所述地震軌跡的所有源位置處被 評(píng)價(jià)，或(iii)接收器濾波，其在所述地震軌跡的所有接收器位置處被評(píng)價(jià)；以及(c)選擇所述地表一致性分量中的一個(gè)或更多個(gè)，并使用它們通過數(shù)值反演來計(jì)算作 為(X，1, ζ)位置的函數(shù)的剪切波速或由此可導(dǎo)出的地球的其它物理性質(zhì)的近地表地球模 型。
2.根據(jù)權(quán)利要求1所述的方法，其進(jìn)一步包括對(duì)于所述軌跡中的一些或全部以及對(duì) 應(yīng)的地震源和接收器位置，計(jì)算穿過每個(gè)中間單元從所述源位置到所述接收器位置的射線 路徑距離，并且其中通過對(duì)來自所述地震勘測(cè)和表面波預(yù)測(cè)的數(shù)據(jù)軌跡的單獨(dú)比較或迭代 比較，地表一致性分量的解優(yōu)化所述分量，其中所述表面波預(yù)測(cè)通過采用假設(shè)的或迭代更 新的分量和被計(jì)算的射線路徑距離來計(jì)算。
3.根據(jù)權(quán)利要求2所述的方法，其進(jìn)一步包括重復(fù)步驟(a)到(b)至少一次，其中(a) 中的所述單元?jiǎng)澐直桓淖?，以改善預(yù)測(cè)的波形與勘測(cè)數(shù)據(jù)波形之間的一致性。
4.根據(jù)權(quán)利要求1所述的方法，其中每個(gè)軌跡的所述分量包含源分量、接收器分量和 從所述源到所述接收器的所述射線路徑穿過的每個(gè)單元的分量。
5.根據(jù)權(quán)利要求1所述的方法，其中在所述地震數(shù)據(jù)中識(shí)別至少兩種不同的模態(tài)，并 且所述地表一致性分量包括多模態(tài)傳播的表面波的參數(shù)。
6.根據(jù)權(quán)利要求5所述的方法，其中在所述數(shù)據(jù)被轉(zhuǎn)換到頻域或頻率波數(shù)域中之后， 從所述數(shù)據(jù)中識(shí)別所述不同的模態(tài)。
7.根據(jù)權(quán)利要求1所述的方法，其中所述地表一致性分量還是方位角或傳播方向的函數(shù)。
8.根據(jù)權(quán)利要求1所述的方法，其中所述地表一致性分量包括速度頻散或衰減的參數(shù)。
9.根據(jù)權(quán)利要求1所述的方法，其中在(c)中所述數(shù)值反演通過將價(jià)值函數(shù)最小化來 被執(zhí)行，所述價(jià)值函數(shù)涉及模型表達(dá)式與所述勘測(cè)數(shù)據(jù)之間的差的一組軌跡的總和，其中 所述模型表達(dá)式是所述地表一致性分量的函數(shù)。
10.根據(jù)權(quán)利要求1所述的方法，其中模型優(yōu)化算法被用于求解所述地表一致性分量， 其中所述模型是單一地滾波模態(tài)或多模態(tài)的波形的模型。
11.根據(jù)權(quán)利要求1所述的方法，其中求解所述地表一致性分量涉及使用一個(gè)或更多 個(gè)約束、阻尼、歸一化、權(quán)重和正則化。
12.根據(jù)權(quán)利要求1所述的方法，其中(c)包括使用所述地表一致性分量來計(jì)算表面波 波形，用于與相應(yīng)的被記錄的地震數(shù)據(jù)比較。
13.根據(jù)權(quán)利要求2所述的方法，其中(c)中的所述反演是非線性優(yōu)化問題，其中所 述模型是近地表速度分布模型，并且所述模型的參數(shù)至少包括近地表層的層厚和層剪切模 量，并且所述反演從被優(yōu)化的地表一致性分量提取每種表面波模態(tài)的速度頻散曲線，并將 它們與所述速度分布模型預(yù)測(cè)的速度頻散曲線比較。
14.根據(jù)權(quán)利要求1所述的方法，其中通過加窗、通過濾波、通過互相關(guān)濾波或通過互 相關(guān)，之后通過將所述數(shù)據(jù)疊加或轉(zhuǎn)換到除了其中記錄數(shù)據(jù)的時(shí)空域之外的域中，來限定 所述地震數(shù)據(jù)。
15.根據(jù)權(quán)利要求1所述的方法，其中所使用的所述地震數(shù)據(jù)軌跡對(duì)應(yīng)于多個(gè)勘測(cè)源 位置，通過多個(gè)接收器來記錄每個(gè)源脈沖。
16.根據(jù)權(quán)利要求1所述的方法，其中所述地表一致性分量中的至少一個(gè)是頻域中的 傳遞函數(shù)或時(shí)域中的脈沖響應(yīng)。
17.根據(jù)權(quán)利要求1所述的方法，其中所述地表一致性分量是頻率的復(fù)合函數(shù)，即具有 實(shí)部和虛部的函數(shù)。
18.根據(jù)權(quán)利要求17所述的方法，其中復(fù)合地表一致性分量包括源振幅和相位、以及 被所述軌跡射線路徑穿過的每個(gè)單元和被處理的每個(gè)地滾波模態(tài)的復(fù)合慢度。
19.根據(jù)權(quán)利要求1所述的方法，其進(jìn)一步包括將來自(c)的剪切模量的計(jì)算的近地表 地球模型用于所述地下區(qū)域中或所述地下區(qū)域上的改進(jìn)工程或改進(jìn)設(shè)計(jì)。
20.根據(jù)權(quán)利要求1所述的方法，其進(jìn)一步包括將來自(c)的剪切波速的被計(jì)算的近地 表地球模型用于所述地下區(qū)域的碳?xì)浠衔镌u(píng)估或從所述地下區(qū)域提取碳?xì)浠衔铩?br> 21.一種用于從地下區(qū)域產(chǎn)生碳?xì)浠衔锏姆椒?，其包?a)獲得包括S-波地震數(shù)據(jù)的所述地下區(qū)域的地震圖像，其中所述S-波數(shù)據(jù)被處理以 便用剪切波模量值成像，所述剪切波模量值通過采用被并入此處以供參考的權(quán)利要求1所 描述的方法而得到；以及(b)至少部分在來自(a)的所述地震圖像上的基礎(chǔ)上，鉆一個(gè)井到所述地下區(qū)域中，并 且從所述井產(chǎn)生碳?xì)浠衔铩?br> 22.一種用于根據(jù)地震數(shù)據(jù)軌跡確定彈性剪切波傳播速度的近地表地球模型的方法， 所述方法包括為多個(gè)軌跡評(píng)估頻域中的至少兩個(gè)地表一致性傳遞函數(shù)或在時(shí)域中的脈沖 響應(yīng)，并且數(shù)值地反演一個(gè)或更多個(gè)所述地表一致性傳遞函數(shù)或脈沖響應(yīng)以計(jì)算作為(X， y，z)位置的函數(shù)的彈性剪切波傳播速度的近地表地球模型，當(dāng)?shù)卣鸨砻娌◤脑磦鞑サ浇邮?器時(shí)，所述至少兩個(gè)地表一致性傳遞函數(shù)或所述脈沖響應(yīng)表示地震表面波波形的變化。
全文摘要
本發(fā)明涉及一種將地震表面波變換為土壤的動(dòng)態(tài)性質(zhì)的深度分布的方法，例如土壤剪切模量或波速或波的衰減。本發(fā)明使用地震數(shù)據(jù)和兩級(jí)反演法，地震數(shù)據(jù)優(yōu)選地是用多個(gè)源和接收器記錄。首先，表面波波形的變化優(yōu)選地被分解(303)成每個(gè)源、每個(gè)接收器以及地表的每個(gè)小區(qū)域(301)的地表一致性傳遞函數(shù)。之后，每個(gè)區(qū)域的傳遞函數(shù)被反演(308)以確定作為深度的函數(shù)的土壤性質(zhì)或近地表性質(zhì)(例如剪切模量)。該方法能夠求解具有垂直變化性質(zhì)和橫向變化性質(zhì)的介質(zhì)的表面波的復(fù)合多模態(tài)特性。避免了傳統(tǒng)方法在求解中錯(cuò)誤識(shí)別地滾波模態(tài)或假設(shè)橫向均勻的土壤性質(zhì)的錯(cuò)誤和限制。
文檔編號(hào)G01V1/00GK102112894SQ200980130779
公開日2011年6月29日 申請(qǐng)日期2009年6月19日 優(yōu)先權(quán)日2008年8月11日
發(fā)明者C·E·克羅恩 申請(qǐng)人:?？松梨谏嫌窝芯抗?br>