[發明專利]一種多點地質統計學沉積相建模方法

有效

• 申請號: 201610265355.2
• 申請日: 2016-04-26
• 公開（公告）號: CN107316341B
• 公開（公告）日: 2020-06-23
• 發明（設計）人: 廉培慶;李劍;譚學群;牛宇虹;張俊法;徐樟有;張文彪;高慧梅;杜秀娟
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G06T17/05
• 分類號: G06T17/05
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 張文娟;朱繪

1.一種多點地質統計學沉積相建模方法，包括以下步驟：

S100，以地質概念模型為指導對研究區進行模式擬合，建立原型地質模型；

S200，根據原型地質模型建立符合沉積模式的訓練圖像并創建搜索樹，對原型地質模型進行變差函數分析得到各沉積相主方向的變程，以獲取幾何因子；

S300，判斷區域空間中是否存在條件數據：

若存在條件數據，則建立網格，將條件數據分配到最近的網格節點，引入幾何因子，在每個訪問節點的周圍構建幾何因子空間，通過指示克里金方法確定幾何因子空間內待估點處的沉積相規模概率作為幾何因子概率P(A|C)，P(A|C)表示在幾何因子空間內以事件C作為條件的A發生的規模條件概率；

S400，指定能夠訪問所有網格節點的隨機訪問路徑；

S500，對于每個訪問節點，判斷是否為條件數據：

若否，執行步驟S510；

若是，確定該節點的沉積相類型，執行步驟S520；

S510，判斷該節點處條件數據事件在搜索樹中是否有記錄：

若有記錄，執行步驟S520；

若無記錄，執行步驟S540；

S520，根據以該節點為中心的條件數據事件在訓練圖像中出現的重復次數得到沉積相估計概率作為該點的局部條件概率P(A|B)，P(A|B)表示數據事件B作為條件在搜索樹中得到的A發生的類型條件概率；

S530，通過校正比例恒等式將局部條件概率P(A|B)與幾何因子概率P(A|C)相結合，得到條件累積概率P(A|B,C)，由此獲得條件概率分布函數，執行步驟S550；

S540，令該節點的條件概率為幾何因子概率P(A|C)，由此獲得條件概率分布函數；

S550，在條件概率分布函數中提取一個模擬值，將其當作實際的觀測數據加入到原來的條件數據集中，作為隨后其他節點的模擬條件；

S600，在訪問完所有節點并確定了所有節點的模擬結果的基礎上，根據沉積相帶的空間結構和幾何形態，得到所需的沉積相模型。

2.根據權利要求1所述的多點地質統計學沉積相建模方法，其特征在于，步驟S300包括以下小步驟：

S310，在每個訪問節點的周圍構建幾何因子空間；

S320，在幾何因子空間內對沉積相類型進行指示變換；

S330，計算每種指示變量的指示變差函數；

S340，利用指示克里金方法計算幾何因子空間內待估點的沉積相估計概率；

S350，確定沉積相的任意次序，將各沉積相主變程進行歸一化處理，獲得待估點沉積相展布規模的概率作為幾何因子概率P(A|C)。

3.根據權利要求1所述的多點地質統計學沉積相建模方法，其特征在于，所述步驟S300中，構建的幾何因子空間滿足以下條件：

①幾何因子空間不小于先驗地質模型相帶連續最小閾，不超過先驗地質模型相帶連續最大閾；

②幾何因子空間內盡可能包含所有已知相帶的條件數據；

③待估點位于幾何因子空間中心。

4.根據權利要求1所述的多點地質統計學沉積相建模方法，其特征在于，所述步驟S520包括以下小步驟：

S521,根據搜索樹獲得條件數據事件和訪問節點為指定狀態在訓練圖像中同時出現的重復次數；

S522,判斷條件數據事件在訓練圖像中的重復次數是否大于給定的重復次數閾值：

若否，去掉距離中心位置最遠的條件數據，返回步驟S521；

若是，計算該節點的條件概率作為局部條件概率P(A|B)。

5.根據權利要求1所述的多點地質統計學沉積相建模方法，其特征在于，所述步驟S530中，校正比例恒等式為：

P(A|B,C)＝1/(1+((1-P(A|C))/P(A|C))×((1-P(A|B))/P(A|B))/(1-P(A))/P(A)))

式中，P(A)為事件A發生的邊緣概率。

6.根據權利要求1所述的多點地質統計學沉積相建模方法，其特征在于，所述步驟S550中，采用隨機方法對條件概率分布函數曲線進行取樣，取樣結果作為模擬值，并將其當作實際觀測數據加入到條件數據集中。

7.根據權利要求6所述的多點地質統計學沉積相建模方法，其特征在于，所述模擬值同時包含沉積相展布規模和相帶空間結構的信息。