[發明專利]一種酸蝕裂縫粗糙程度定量表征方法

有效

• 申請號: 201911188876.2
• 申請日: 2019-11-28
• 公開（公告）號: CN110909477B
• 公開（公告）日: 2020-08-11
• 發明（設計）人: 陳遲;王守信;郭建春;盧聰;王琨;劉彧軒;賴杰;文帝霖;鐘燁
• 申請（專利權）人: 西南石油大學
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06T7/00;G06T7/62;G06T7/64
• 代理公司: 北京中索知識產權代理有限公司 11640
• 代理人: 房立普

1.一種酸蝕裂縫粗糙程度定量表征方法，其特征在于，包括以下步驟：

S1、采集碳酸鹽巖儲層段的井下巖心或者同層位露頭巖石，并制成多塊尺寸相同的原始方形巖板；

S2、在原始方形巖板中部沿巖板長度方向預制劃痕，再將預制劃痕的原始方形巖板劈裂為一對具有粗糙表面形態的劈裂巖樣；

S3、使用三維激光掃描儀掃描劈裂巖樣的上下兩個表面來獲取三維點云數據，再采用標準差濾波法對三維點云數據進行降噪處理；

S4、采用克里金插值法對步驟S3中獲取的三維點云數據進行插值，插值時所設置的步長分別為1mm×1mm和0.1mm×0.1mm，插值后獲得的三維點云數據分別為酸刻蝕前巖板一階點云數據、酸刻蝕前巖板二階點云數據；

S5、利用步驟S4中獲取的酸刻蝕前巖板一階點云數據、酸刻蝕前巖板二階點云數據，分別計算酸刻蝕前巖板的一階面積迂曲度和二階面積迂曲度；利用步驟S4中獲取的酸刻蝕前巖板二階點云數據，計算巖板長寬兩個方向的平均線迂曲度；

S6、對步驟S2中的具有粗糙表面形態的劈裂巖樣開展酸刻蝕物理模擬實驗，獲取酸刻蝕后巖板，并使用三維激光掃描儀掃描酸刻蝕后巖板來獲取酸刻蝕后的三維點云數據并剔除噪點數據；

S7、再采用克里金插值法對步驟S6中獲取的酸刻蝕后的三維點云數據進行插值，插值時所設置的步長分別為1mm×1mm和0.1mm×0.1mm，插值后獲得的三維點云數據分別為酸刻蝕前后板一階點云數據、酸刻蝕后巖板二階點云數據；

S8、利用步驟S7中獲取的酸刻蝕后巖板一階點云數據、酸刻蝕后巖板二階點云數據，分別計算酸刻蝕后巖板的一階面積迂曲度和二階面積迂曲度；利用步驟S7中獲取的酸刻蝕后巖板二階點云數據，計算巖板長寬兩個方向的平均線迂曲度；

S9、最后根據步驟S5和步驟S8中獲取的數據來定量表征裂縫面在酸蝕前后不同方向、不同尺度下的粗糙度程度。

2.根據權利要求1所述的一種酸蝕裂縫粗糙程度定量表征方法，其特征在于，所述步驟S3中采用標準差濾波法對三維點云數據進行降噪處理的具體過為：

步驟S31、在掃描獲取的裂縫面三維點云數據數組{X_i、Y_i、Z_i}中計算每個點{X_n、Y_n、Z_n}與之臨近的8個鄰域點的距離，8個鄰域點的x和y方向坐標分別為{X_n-1、Y_n-1}、{X_n-1、Y_n}、{X_n-1、Y_n+1}、{X_n、Y_n-1}、{X_n、Y_n+1}、{X_n+1、Y_n-1}、{X_n+1、Y_n}、{X_n+1、Y_n+1}；

步驟S32、對步驟S31中計算的距離進行統計，并計算平均距離的均值u和標準差r；

步驟S33、判斷點{X_n、Y_n、Z_n}到8個鄰域點的平均距離u與距離閾值d＝u±5r之間的關系；若大于距離d，則為噪點進行去除。