[發明專利]一種混合二維地震波走時計算方法

在審

• 申請號: 201810077621.8
• 申請日: 2018-01-23
• 公開（公告）號: CN108072897A
• 公開（公告）日: 2018-05-25
• 發明（設計）人: 孫輝;孫章慶;張志厚;張平;孟繁昌;鄒時貴;夏龍龍;鄧瑞;李夢;黎佳賓;劉丁毅;鐘雨田;張杰
• 申請（專利權）人: 西南交通大學
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 成都信博專利代理有限責任公司 51200
• 代理人: 張輝

1.一種混合二維地震波走時計算方法，其特征在于，包括以下步驟：

步驟1：讀入相關參數文件、速度模型，其中，所述參數文件包含速度模型的網格點數、網格間距和震源位置；

步驟2：從炮點沿著不同方向發射射線并計算中心射線信息；

步驟3：使用波前構建法計算射線覆蓋范圍內網格節點走時；

步驟4：對所有網格點進行分類，將它們分為接受點、窄帶點或遠離點，即如果一個網格點的地震波走時已經計算過，并且它周圍點的走時全部計算過，則這個點為接受點；如果一個網格點的地震波走時計算過，而它周圍點的走時并非全部計算過，則這個點為窄帶點；如果一個點的走時沒有被計算過，則這個點為遠離點；

步驟5：將所有窄帶點移入窄帶內，構建初始窄帶；

步驟6：擴展窄帶，直至窄帶為空；在這一過程中，網格點走時是通過迎風差分法求解二維程函方程獲得的，所述二維程函方程為：

|▽τ|＝s

其中，τ為地震波走時，s為慢度，▽為梯度符號，上述公式中梯度項的迎風差分表達式為：

| ▿ τ | = max ( D i , j - x τ , 0 ) 2 + min ( D i , j + x τ , 0 ) 2 + max ( D i , j - z τ , 0 ) 2 + min ( D i , j + z τ , 0 ) 2 1 / 2 ]]>

其中，分別為走時在網格節點(i,j)處x或z方向上的向前或向后差分表達式，具體形式分別為：

2.如權利要求1所述的一種混合二維地震波走時計算方法，其特征在于，在步驟2中，使用龍哥庫塔法求解運動學射線追蹤方程組得到中心射線信息，如下式所示：

dx i d τ = v 2 · p i dp i d τ = ∂ ∂ x i ( 1 v ) = - v - 1 ∂ v ∂ x i ]]>

其中，xi表示空間位置，pi表示慢度分量，τ表示地震波走時，v表示離散點處的速度值。