[發明專利]基于地質構造散點數據自動求取含烴面積的方法

有效

• 申請號: 202011592388.0
• 申請日: 2020-12-29
• 公開（公告）號: CN112687008B
• 公開（公告）日: 2022-12-06
• 發明（設計）人: 侯東梅;李彥來;郭敬民;權勃;張小龍;劉春艷;李博;汪巍;楊東東;孫恩慧
• 申請（專利權）人: 中國海洋石油集團有限公司;中海石油（中國）有限公司天津分公司
• 主分類號: G06T17/20
• 分類號: G06T17/20;G06T17/05;G06T7/62;G06Q50/02
• 代理公司: 天津三元專利商標代理有限責任公司 12203
• 代理人: 胡婉明

1.一種基于地質構造散點數據自動求取含烴面積的方法，其特征在于，包括如下步驟：

第一步：利用道格拉斯-普克抽稀算法對原始構造散點進行抽稀，去掉冗余的點；

第二步：依據第一步抽稀后的散點，求取該散點的外包矩形，作為構造等值線網格的邊界；

第三步：求取所有散點間距集合，從集合中獲得最小間距，取整后作為網格的縱橫間距(或者說網格步長)；

第四步：基于掃描線算法原理，求抽稀后的等值線與各條網格線的交點集合；

第五步：基于所有交點集合，依次計算交點所在單元格的頂點值；先利用加權平均法計算單元格左上頂點，再依據左上頂點所在四個單元格的邊，進一步利用線性關系計算各邊的另一個頂點，如此遞歸，直到完成所有交點所在單元格的頂點值；

第六步：將第五步獲得的所有單元格頂點作為樣本點，利用反距離插值算法對未賦值的單元格頂點賦值，得到構造等值線網格；

第七步：根據流體界面深度，在第六步中得到的構造網格中追蹤該深度的平面位置，進而得到含烴面積。

2.根據權利要求1所述的基于地質構造散點數據自動求取含烴面積的方法，其特征在于，

所述第一步是利用現有的道格拉斯-普克抽稀算法(Douglas–Peucker algorithm)對原始構造散點進行抽稀，去掉冗余的點；

所述第二步是依據抽稀后的散點，求取其外包矩形，作為構造等值線網格的邊界是依次取出抽稀后的點，比較其X，Y坐標，將獲得的最小X、最大X、最小Y、最大Y坐標作為外包矩形四條邊的坐標值；

所述第四步的基于掃描線算法原理，求抽稀后的等值線與各條網格線的交點集合是取出每條等值線的線段，求取線段兩個端點內的橫、縱網格線，利用相似三角形快速計算線段與橫、縱坐標的交點，即為網格的交點；

第六步：把所有獲得的單元格頂點作為樣本點，利用反距離插值算法對未賦值的單元格頂點賦值，得到構造等值線網格，對任意一個單元格頂點A，以三個樣本點為例，分別為樣本點B、樣本點C和樣本點D，獲得頂點A賦值的反距離插值算法為：

設：B1＝1/(A到B的距離)，C1＝1/(A到C的距離)，D1＝1/(A到D的距離)，單元格A頂點的值Z＝B的Z值*B1/(B1+C1+D1)+C的Z值*C1/(B1+C1+D1)+D的Z值*D1/(B1+C1+D1)；

第七步：根據流體界面深度，在得到的構造網格中追蹤該深度的平面位置，進而得到含烴面積是：在網格中提取所有深度大于流體界面深度的網格頂點及構造等值線，根據等值線兩側構造點的深度，將一側為大于界面深度、一側為小于界面深度的線提取出來，對其做閉合處理，得到含烴面積。