層析成像技術論文
層析成象是在物體外部發(fā)射物理信號,接收穿過物體且攜帶物體內部信息,利用計算機圖象重建方法,重現(xiàn)物體內部一維或三維清晰圖象。 下面是小編精心推薦的一些層析成像技術論文,希望你能有所感觸!
層析成像技術論文篇一
前言
層析成象是在物體外部發(fā)射物理信號,接收穿過物體且攜帶物體內部信息,利用計算機圖象重建方法,重現(xiàn)物體內部一維或三維清晰圖象。 層析成象技術最大的特點是在不損壞物體的條件下,探知物體內部結構的幾何形態(tài)與物理參數(shù)(如密度等)的分布。層析成象與空間技術、遺傳工程、新粒子發(fā)現(xiàn)等同列為70年代國際上重大科技進展。層析成像應用非常廣泛,如醫(yī)學層析的核磁共振成像技術、工業(yè)方面的無損探傷、 在軍事工業(yè)中,層析成象用于對炮彈、火炮等做質量檢查、在石油開發(fā)中被用于巖心分析和油管損傷檢測等,層析成象是在物體外部發(fā)射物理信號,接收穿過物體且攜帶物體內部信息,利用計算機圖象重建方法,重現(xiàn)物體內部一維或三維清晰圖象。 聲波層析成像技術
聲波層析成像方法所研究的主要內容,一個是正演問題,即射線的追蹤問題,是根據(jù)已知速度模型求波的初至時間的問題;另一個問題就是反演問題,即根據(jù)波的初至時間反求介質內部速度或者慢度分布的問題。層析成像效果的好壞與解正演問題的正演算法和解反演問題的反演算法都有直接的關系。 論文詳細研究聲波層析成像的射線追蹤算法,重點探討了基于Dijkstra算法的Moser曲射線追蹤算法,并用均勻介質模型、空洞模型、低速斜斷層等模型使用Moser曲射線追蹤時的計算精度與計算效率,發(fā)現(xiàn)了內插節(jié)點是影響Moser曲射線追蹤效果的主要因素,得到了內插節(jié)點數(shù)為5~7之間,計算速度較快,計算精度較高。模型試算的結果表明,正演采用內插10個節(jié)點,
反演過程中采用內插5個節(jié)點,效果最佳。 在層析成像正演算法的基礎上,詳細研究了誤差反投影算法(BPT)、代數(shù)重建法(ART)、聯(lián)合迭代法(SIRT);研究了非線性問題線性化迭代的最速下降法、共軛梯度法(CG);重點推導和建立了層析成像的高斯—牛頓反演法(GN);詳細研究了非線性最優(yōu)化的蒙特卡洛法(MC)、模擬退火法(SA)、遺傳算法(GA);研究了將非線性全局最優(yōu)化和線性局部最優(yōu)化方法相結合的混合優(yōu)化方法,探討了基于高斯牛頓和模擬退火相結合(GN-SA)混合優(yōu)化算法。在此基礎上,以速度差為10%的低速斜斷層模型為例,詳細探討了線性化算法SIRT、GN;非線性最優(yōu)化算法SA、GA以及混合優(yōu)化算法GN-SA五種算法對該模型的計算結果,并探討了直射線和Moser曲射線追蹤的反演效果。數(shù)值試驗表明,基于Moser曲射線追蹤的高斯—牛頓反演法的層析成像效果最佳,計算效率最高。 采用基于Moser曲射線追蹤的高斯—牛頓法,對速度差為25%的等軸狀空洞構造、速度差為33%的不連通空洞模型、速度差為33%的高速巖脈進行了反演試算,對于這些理論模型,高斯—牛頓法均取得了較好的成像效果。為進一步驗證各種層析成像法,在實驗室制作了水泥臺和石膏板實物模型,并分別在水泥臺中央制作一個方形空洞,在石膏板中央制作一個倒“L”形空洞。對這兩個實物模型進行了實測,對測量的數(shù)據(jù),用高斯—牛頓法進行層析成像反演,均取得了較好的成像效果。 通過本文的研究和數(shù)值試驗,得到了以下結論:(1)基于直射線追蹤方法,適用較為簡單的地質體,亦或是測量精度要求不高的問題。由于直射線追蹤方法在成像過程中,只需要追蹤一次就可以
求得距離矩陣,這樣它的成像速度比較快,而基于Moser曲射線追蹤的SIRT成像反演法,在迭代過程中,需要不斷地進行距離矩陣的更新,計算速度相當慢。因此,當實際地質情況比較簡單時候,可以考慮先采用直射線方法進行成像,然后采用曲射線追蹤進行構造精細解釋。(2)基于Moser曲射線的射線追蹤方法,追蹤效果與內插節(jié)點的數(shù)目有較大的關系。(3)在Moser曲射線追蹤基礎上,結合SIRT層析成像反演方法,對正演模擬的旅行時間進行了成像反演。當速度差異小于15%時,基于曲射線的SIRT層析成像反演結果與直射線情況下差異不大;當速度差高達33%時,基于直射線的SIRT層析成像方法對高速異常區(qū)的成像效果仍然比較好。但當速度差異大于67%時,基于直射線的SIRT成像效果比較差,但是基于Moser曲射線的追蹤方法,仍然可以給出比較好的成像效果。(4)基于Moser曲射線的高斯—牛頓反演層析成像方法,進行了反演。(5)通過方法的比較,高斯牛頓法一般只需要迭代2~3次,就可以得到比較好的成像效果。而一般的SIRT成像方法,需要迭代10次左右才能得到比較好的成像效果。(6)理論模型的數(shù)值試驗表明,盡管非線性最優(yōu)化方法在理論上可以收斂到全局最優(yōu)解,但是在實踐過程中,非線性最優(yōu)化算法SA、GA以及混合最優(yōu)化方法,目前仍然存在搜索次數(shù)太大,搜索時間過長等弊病而無法實用。 這些理論模型和實測資料的反演結果,為建筑物構件等的無損檢測提供了理論依據(jù)。
一.井間地震層析成像技術
井間地震層析成像技術是利用地震波在不同方向投射的波場信息, 對地下介質內部精細結構(速度、衰減系數(shù)、反射系數(shù)等的分布)進行成像, 以其分辨率高、解析成果直觀等特點, 廣泛應用于工業(yè)及 民用建筑、公路、鐵路、環(huán)境等方面工程地質勘察中。井間層析成像可分為基于射線理論的走時層析成像和波動理論的繞射層析成像兩類。井間地震波場信息豐富復雜, 波場識別和分離比較困難, 而直達波至相對簡單, 故工程勘察中常常采用基于射線理論的直達波至走時層析成像。井間地震層析成像的核心問題是: 至波走時和線路徑計算, 即正演問題; 過不同的重建算法進行成像、解釋即反演問題。筆者正演采用最短路徑法射線追蹤, 反演采用基于正交分解最小二乘法(簡稱LSQR 算法)的反演算法。
1、 最短路徑法射線追蹤
它是基于Fermat最小旅行時原理和網絡理論中的最短路徑算法來實現(xiàn)。把地下介質離散成若干小單元體, 并在各單元邊界上設置一些節(jié)點, 地下速度模型就可表示成由這些節(jié)點以及它們之間的連線 所形成的網絡。網絡中的每一個節(jié)點只能與彼此相鄰的節(jié)點連接。相鄰節(jié)點之間的連接權等于地震波沿其連線傳播的旅行時。一條路徑是由相互連接的節(jié)點序列組成的, 沿著該路徑的旅行時為該路徑上所有連接權之和。從一個節(jié)點到另一個可能有無數(shù)條路徑, 按照Fermat原理, 把旅行時最小的路徑近似為地震波傳播通過的射線。網絡中, 速度場分布在離散的節(jié)點上。相鄰節(jié)點之間的旅行時為他們之間歐氏距離與其平均慢度之積。將波陣面看成是由有限個離散點次級源組成, 對于某個次級源(即某個網格節(jié)點) , 選取與其所有相鄰的點(鄰域點)組成計算網格點; 由一個源點出發(fā), 計算出從源點到計算網格點的透射走時、射線路徑和射線長度; 然后把除震源之外的所有網格點相繼當作次級源, 選取該源點相應的計算網格點,計算出從次級源點到計算網格點的透射走時、射線路徑和射線長度; 將每次計算出來的走時加上從震源到次級源的走時, 作為震源點到該網格節(jié)點的走時, 記錄下相應射線路徑位置及射線長度。
2、層析成像的反演LSQR 算法
LSQR 方法是Pa ige 和Sanders在1982 年提出的, 它是利用Lanczos迭代法求解最小二乘問題的一種方法。LSQR 方法具有計算量小的優(yōu)點, 并且能很容易地利用矩陣的稀疏性簡化計算, 因而適合 求解大型稀疏問題。
LSQR 是目前層析成像中常用的方法, 在迭代過程中,它只涉及到非零元素, 占有用存儲空間少,運算速度快, 運算穩(wěn)定, 迭代次數(shù)少。
3、結論
(1)井間地震可以把震源和檢波器的排列直接布置在鉆孔中, 這可以使接收到的地震信號能保留更高頻率的有效成分, 為提高分辨率打下物理基礎。但是, 也正是由于震源和檢波器位置分布相對固定,且數(shù)目有限, 故獲得的不同角度的數(shù)據(jù)量有限, 不能夠像醫(yī)學CT 那樣獲得全方位的數(shù)據(jù), 這使得層析成像的解不唯一, 故進行圖像解釋時, 必須通過鉆孔資料加以約束。
(2)對規(guī)模較小的破碎帶、節(jié)理裂隙發(fā)育等結構缺陷的低速異常體的探測, 由于對地震波的走時影響較小, 不足以改變地震波的射線路徑, 則層析成像的分辨率依然達不到, 不能夠將這些低速異常體分辨出來。
(3)層析成像質量不僅與異常體大小有關, 還與孔間距及孔深與孔間距之比有關, 外業(yè)數(shù)據(jù)采集質量、拾取初至波走時誤差、反演算法亦對反演圖像解釋產生直接影響。
(4)為提高探測精度, 可嘗試利用VSP和井間CT數(shù)據(jù)采集方式的相似性, 設置科學的野外觀測系統(tǒng), 同時開展工程VSP 與地震CT 聯(lián)合探測方法研究, 進行聯(lián)合反演, 這也將是我們今后努力的方向。
結束語
層析成像方法在地球物理探測方面還有許多廣泛的應用,這里我們就不一一討論。如何利用層析成像技術為地球物理探測技術的發(fā)展
做出貢獻,是每個物探人應該努力的方向。
點擊下頁還有更多>>>層析成像技術論文