[發明專利]基于同位多相協同克里金的地震彈性參數相控建模方法

在審

• 申請號: 201710098888.0
• 申請日: 2017-02-23
• 公開（公告）號: CN106772587A
• 公開（公告）日: 2017-05-31
• 發明（設計）人: 張宏兵;韓飛龍;尚作萍;郭強;曹呈浩
• 申請（專利權）人: 河海大學
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 南京經緯專利商標代理有限公司32200
• 代理人: 徐瑩

1.基于同位多相協同克里金的地震彈性參數相控建模方法，其特征在于，包括以下步驟：

步驟S1、在目標層內采用等比例網格剖分方法對建模網格進行剖分，在相同比例厚度的位置上形成若干個平面，作為彈性參數建模的目標位置；

步驟S2、分別提取平面上已知點測井彈性參數以及分布于整個區域的地震屬性參數和沉積相信息作為建模數據，即選擇測井參數的低頻參數背景值作為作為主變量，選擇地震屬性參數作為第一協同變量，及利用各向異性擴散法對沉積相信息進行概化處理得到沉積相數字模型作為第二協同變量；

步驟S3、將步驟S2所得主變量、第一協同變量及第二協同變量作為計算參數，利用同位多相協同克里金方法進行插值計算，獲得平面上所有待估計點的參數值作為彈性參數的多源數據融合建模結果；

步驟S4、使用各向異性擴散法對步驟S3所得多源數據融合建模結果進行概化處理，得到概化處理后的地震彈性參數建模結果。

2.根據權利要求1所述的基于同位多相協同克里金的地震彈性參數相控建模方法，其特征在于，所述步驟S2中選擇測井參數的低頻參數背景值作為主變量，具體包括：

使用多項式擬合方法對提取的測井彈性參數進行平滑處理，根據不同的多項式擬合次數得到若干組平滑后的低頻參數背景值；

選擇大于預設低頻參數背景值變化范圍的一組低頻參數背景值作為主變量。

3.根據權利要求1所述的基于同位多相協同克里金的地震彈性參數相控建模方法，其特征在于，所述步驟S2中，根據主變量與地震屬性參數間的相關系數選擇地震屬性參數作為第一協同變量。

4.根據權利要求1所述的基于同位多相協同克里金的地震彈性參數相控建模方法，其特征在于，所述步驟S2中利用各向異性擴散法對沉積相信息進行概化處理，采用公式：

$Z_{i,j}^{(m+1)}={\mathrm{AZ}}_{i,j}^{\left(m\right)}+\frac{B}{8}\underset{k,l=-1}{\overset{1}{\Σ}}{C}_{i+k,j+l}\▿Z_{i+k,j+l}$

其中，Z為待處理點參數值，▽Z為各方向臨近點梯度值，A、B、C為各方向的權重系數，i、j、k為所處理點的相對位置，m為點的迭代次數。

5.根據權利要求1所述的基于同位多相協同克里金的地震彈性參數相控建模方法，其特征在于，所述步驟S3中利用同位多相協同克里金方法插值進行計算，具體包括：

步驟S31、利用同位協同多相克里金方法，建立計算平面上待估計點的參數值模型；

步驟S32、根據無偏最優條件得到克里金方程組，及將克里金方程組轉換為矩陣；

步驟S33、計算不同位置待估計點的各參數間的協方差和互協方差，以得到步驟S32所得矩陣中所有變量值；

步驟S34、通過求解步驟S32所得矩陣得到各變量的權重系數，將其代入步驟S31所建立待估計點的參數值模型計算得到待估計點的參數值。