[發明專利]一種確定地質變量分位值的方法和系統

有效

• 申請號: 201710107109.9
• 申請日: 2017-02-27
• 公開（公告）號: CN108510134B
• 公開（公告）日: 2022-02-01
• 發明（設計）人: 金之鈞;盛秀杰;徐忠美;吳博
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G06Q10/06
• 分類號: G06Q10/06;G06Q50/02
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建

1.一種確定地質變量分位值的方法，其特征在于，所述方法包括以下步驟：

步驟S10：地質變量確定及數據準備，根據具體任務確定待分析的地質變量，針對確定的待分析的地質變量進行數據統計；

步驟S20：建立地質變量變差模型，根據步驟S10得到的地質變量的統計數據，建立變差模型；

步驟S30：變差函數分析，基于步驟S20中的變差模型，交互調整變差關鍵參數，進行變差分析，得到地質變量的等值線分布圖；

步驟S40：工區網格化及編碼，將步驟S30中得到的地質變量的等值線分布圖進行網格化，并對網格進行編碼，得到規則的網格化單元，網格單元包括被等值線覆蓋的網格單元和未被等值線覆蓋的網格單元；

步驟S41：確定感興趣的研究范圍，包括全區范圍、單一局部范圍和多個局部范圍三種情況；

步驟S42：針對感興趣區域，選擇網格化方式，包括行列數設置、行距及列距設置兩種模式，前者是設置固定的行數及列數來劃分網格，后者是設置不同的行距和列距來劃分網格；

步驟S43：基于前述步驟S41和步驟S42，進行網格劃分，記錄網格中心空間位置并按行序對網格進行編碼；

步驟S50：對步驟S40中得到的網格化單元進行地質變量的均值和方差的計算；

步驟S60：隨機模擬抽取網格化單元的變量值及分位值，并進行排序；

其中，步驟S50中，對于所述被等值線覆蓋的網格單元，根據覆蓋網格單元的等值線的參數值給網格單元進行賦值，基于網格的樣本點統計得到均值和方差；對于所述未被等值線覆蓋的網格單元，步驟S50包括以下子步驟：

步驟S51：基于網格劃分的結果，結合已知網格間空間位置關系，確定網格點模擬插值路徑；

步驟S52：采用克里金插值算法計算網格點均值和方差。

2.根據權利要求1所述的確定地質變量分位值的方法，其特征在于，步驟S10中所述的數據統計包括變量參數值和空間坐標。

3.根據權利要求1所述的確定地質變量分位值的方法，其特征在于，步驟S20中所述的變差模型為球狀變差模型、高斯變差模型、線性變差模型或指數變差模型。

4.根據權利要求1所述的確定地質變量分位值的方法，其特征在于，步驟S30中所述的關鍵參數包括塊金值、基臺值和/或變程值。

5.根據權利要求4所述的確定地質變量分位值的方法，其特征在于，步驟S30中所述地質變量的等值線分布圖是通過設置迭代次數、網格參數和分位值來進行校正的。

6.一種確定地質變量分位值的系統，其特征在于，包括：

地質變量確定單元，根據具體任務確定待分析的地質變量，針對確定的待分析的地質變量進行數據統計；

變差模型單元，根據地質變量確定單元得到的地質變量的統計數據，建立變差模型；

變差函數分析單元，基于變差模型單元得到的變差模型，交互調整變差關鍵參數，進行變差分析，得到地質變量的等值線分布圖；

等值線分布圖單元，用于工區網格化及編碼，將變差函數分析單元中得到的地質變量的等值線分布圖進行網格化，并對網格進行編碼，得到規則的網格化單元；

均值和方差的計算單元，對等值線分布圖單元中得到的網格化單元進行地質變量的均值和方差的計算，網格單元包括被等值線覆蓋的網格單元和未被等值線覆蓋的網格單元；以及

變量值及分位值確定單元，用于隨機模擬抽取網格化單元的變量值及分位值，并進行排序；

其中，等值線分布圖單元包括以下模塊：

研究范圍確定模塊，用于確定感興趣的研究范圍，包括全區范圍、單一局部范圍和多個局部范圍三種情況；

網格化模塊，用于針對感興趣區域，選擇網格化方式，包括行列數設置、行距及列距設置兩種模式，前者是設置固定的行數及列數來劃分網格，后者是設置不同的行距和列距來劃分網格；以及

網格單元編碼模塊，基于研究范圍確定模塊和網格化模塊，進行網格劃分，記錄網格中心空間位置并按行序對網格進行編碼；

均值和方差的計算單元中，對于所述被等值線覆蓋的網格單元，根據覆蓋網格單元的等值線的參數值給網格單元進行賦值，基于網格的樣本點統計得到均值和方差；對于所述未被等值線覆蓋的網格單元，均值和方差的計算單元包括以下模塊：

網格點模擬插值路徑模塊，用于基于網格劃分的結果，結合已知網格間空間位置關系，確定網格點模擬插值路徑；以及

均值和方差計算模塊，采用克里金插值算法計算網格點均值和方差。