[發明專利]一種近地表品質因子波形反演方法

在審

• 申請號: 202110509073.3
• 申請日: 2021-05-11
• 公開（公告）號: CN115327638A
• 公開（公告）日: 2022-11-11
• 發明（設計）人: 崔慶輝;尚新民;王修敏;趙勝天;刁瑞;金昌坤;冮明川
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/50
• 分類號: G01V1/50
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 劉亞寧

1.一種近地表品質因子波形反演方法，其特征在于，其包括以下步驟：

步驟1.采集微測井數據；

步驟2.微測井初至波拾取；

步驟3.微測井初至波時差校正；

步驟4.近地表品質因子模型組及修正值初始化；

步驟5.初至波吸收衰減補償；

步驟6.計算各初至波形之間的相似度；

步驟7.相似度或迭代次數判斷，獲得最佳近地表品質因子模型。

2.根據權利要求1所述方法，其特征在于，步驟1中，在近地表深度范圍內采集微測井數據，方法為：通過在井中不同深度布置檢波器，井口激發地震波，或者在井口布置檢波器，在井中不同深度激發地震波，進行微測井并采集數據。

3.根據權利要求1所述方法，其特征在于，步驟2中，拾取微測井各道信號初至時間及波形。

4.根據權利要求1所述方法，其特征在于，步驟4中，近地表模型組是一組近地表品質因子模型Q_i(h_j)，i＝1，2,3，……,N，j＝1,2,3,……,M；N為設置的模型數量，M為設置的近地表品質因子模型深度網格數；

對應于每一個模型Q_i(h_j)有一個用于迭代修改的值，用M_i(h_j)表示，即對應于第i個模型在當前迭代中進行更新的值，其初始值設置為M_i(h_j)＝0，i＝1,2，3,……,N，j＝1，2,3，……,M。

5.根據權利要求1所述方法，其特征在于，步驟5中，根據各近地表品質因子模型進行初至波吸收衰減補償，首先將各道信號初至波形進行傅里葉變換，然后按照以下公式進行近地表吸收衰減補償：

公式中，U(f)為接收(激發)深度為h的初至波信號在頻率域的表示，V(f)為U(f)經過近地表吸收衰減補償后的數據，f為頻率，f₀為參考頻率，n為深度h包含的模型深度網格數，Q_k為第k個模型網格的Q值，Δt_k為模型網格大小與第k個模型網格速度的商。

6.根據權利要求1所述方法，其特征在于，步驟6中，采用以下公式計算各道信號初至波形之間的相似度：

公式中，A為各道信號初至波形相似度，X_j(i)表示第j個信號的第i個采樣點；N為每道信號初至波形的采樣點個數；M為道數。

7.根據權利要求1所述方法，其特征在于，步驟7中，若相似度滿足預設的閾值或迭代次數滿足預設值，獲得最佳近地表品質因子模型；若相似度不滿足預設的閾值或迭代次數不滿足預設值，則更新近地表品質因子模型組，重復步驟5-7，直至相似度滿足預設的閾值或迭代次數滿足預設值，獲得最佳近地表品質因子模型。

8.根據權利要求7所述方法，其特征在于，按照以下公式更新近地表品質因子模型組：

公式中，迭代次數k＝1,2,3,……,I，I為最大迭代次數；為近地表品質因子模型組中第i個模型在第k次迭代完成時的修正值；為近地表品質因子模型組中第i個模型在第k次迭代完成時的值；ω、c₁、c₂均為常數；為近地表品質因子模型組中第i個模型在歷次迭代中獲得最大相似度的模型值；Q^best(h_j)為近地表品質因子模型組中所有模型在歷次迭代中獲得最大相似度的模型值。

9.根據權利要求8所述方法，其特征在于，最佳近地表品質因子模型是迭代終止時的Q^best(h_j)。