[發明專利]一種結合層位信息的微測井方位加權插值建模方法

在審

• 申請號: 201910661853.2
• 申請日: 2019-07-22
• 公開（公告）號: CN110618460A
• 公開（公告）日: 2019-12-27
• 發明（設計）人: 金昌昆;王延光;尚新民;王修銀;王修敏;徐強;孫興剛
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/40
• 分類號: G01V1/40
• 代理公司: 11350 北京科億知識產權代理事務所(普通合伙)
• 代理人: 湯東鳳

1.一種結合層位信息的微測井方位加權插值建模方法，其特征在于，包括以下步驟：

步驟一，確定微測井的主要信息，微測井的主要信息包括微測井速度曲線、解釋層位、工區范圍、反演所需的微測井數、反演的網格、約束權重，然后輸入微測井的主要信息；

步驟二，統計各微測井的距離分布，確定參考距離；

步驟三，基于微測井的參考距離和層位，插值出層位面；

步驟四，統計各層位面的平均值，并將微測井的速度曲線變換為對應層位平均值的速度曲線；

步驟五，基于微測井的參考距離和水平位置，計算徑向權重和方位權重；

步驟六，以微測井的方位系數作為未知量，構建反演方程組，以變換后的微測井速度作為右端項，并添加反演約束項；

步驟七，求解反演方程組，獲得方位系數。

步驟八，逐點進行方位插值，并將各點的插值結果進行層位變換，獲得表層速度模型。

2.根據權利要求1所述的結合層位信息的微測井方位加權插值建模方法，其特征在于：所示步驟二中，對于各微測井，統計出其它微測井與該微測井的距離，基于輸入的反演所需微測井數，確定與之相對應的距離，并作為該微測井的參考距離。

3.根據權利要求2所述的結合層位信息的微測井方位加權插值建模方法，其特征在于，在步驟三中，插值層位面是通過徑向基函數方法計算得到，具體計算方法如下：某層位面在n口微測井處測得的層位深度為s_i,(i＝1,…,n)，以作為目標函數，式中λ_j為第j口微測井的權重系數，s_j為第j口微測井的層位深度，||x_i-x_j||₂表示第i口微測井與第j口微測井的水平距離，||.||₂為二范數運算，x_i、x_j分別為第i、j口微測井的位置；R為徑向權重函數，其形式為：式中，為第i口微測井和第j口微測井的距離與參考距離的比值，r_ref為參考距離；結合最小二乘思想，求解得到各口微測井的權重，并進行全區的層面插值，插值公式如下：其中，s(x)為插值得到的層位面，x為插值點的水平位置。

4.根據權利要求1所述的結合層位信息的微測井方位加權插值建模方法，其特征在于，在步驟四中，以各層位面的平均值為基準，對微測井速度所對應的深度進行線性變換，使對應同一層位的深度統一為該層位平均深度，層位間的深度由上下層位的平均深度線性計算得到，具體計算公式為：其中，H為變換后的深度，h為原始深度，h_i、h_i+1為h對應的上下層位面深度，為相應的層位面的平均深度。

5.根據權利要求1所述的結合層位信息的微測井方位加權插值建模方法，其特征在于，所述步驟五中，每口微測井的徑向權重的計算公式為：式中，R(r)為每口微測井的徑向權重函數，為第i口微測井和第j口微測井的距離與參考距離的比值，x_i、x_j分別為第i、j口微測井的位置，r_ref為參考距離。

6.根據權利要求1所述的結合層位信息的微測井方位加權插值建模方法，其特征在于，在步驟五中，基于微測井的空間位置，對于各微測井，計算該微測井與其它微測井之間的方位角，并以正北、東北、正東、東南、正南、西南、正西、西北這八個主方位系數加權的形式表示相應方位角的權重，方位系數加權所用函數為二階三角B樣條函數，具體形式如下：其中，t∈[0,1]，為由方位角計算得到的參數，計算公式為：其中θ為方位角，{...}為提取t的小數部分的數學運算符。