[發明專利]計算頁巖各向異性參數的方法及計算機可讀存儲介質

在審

• 申請號: 201710508417.2
• 申請日: 2017-06-28
• 公開（公告）號: CN109145340A
• 公開（公告）日: 2019-01-04
• 發明（設計）人: 周楓;劉衛華
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G06F17/50
• 分類號: G06F17/50;G06T17/05
• 代理公司: 北京思創畢升專利事務所 11218
• 代理人: 孫向民;廉莉莉

1.一種計算頁巖各向異性參數的方法，其特征在于，該方法包括：

獲得頁巖巖心的三維圖像；

根據三維圖像中的層理方向，測量出取芯方向與層理的夾角參數θ；

利用所述夾角參數θ計算頁巖各向異性參數。

2.根據權利要求1所述的計算頁巖各向異性參數的方法，其中，利用所述夾角參數θ計算頁巖各向異性參數包括利用所述夾角參數θ，通過公式(1)，計算第五彈性常數：

其中，公式(1)為修正公式，C₁₃為第五彈性常數，C₄₄為第四彈性常數，C₁₁為第一彈性常數，C₃₃為第二彈性常數，ρ為所述頁巖巖芯的密度，V_pθ為所述取芯方向與層理成θ角的頁巖巖芯的縱波速度，D為求解中間量。

3.根據權利要求2所述的計算頁巖各向異性參數的方法，其中，利用所述夾角參數θ計算頁巖各向異性參數還包括如下步驟：

通過頁巖彈性常數計算公式，計算第一彈性常數、第二彈性常數、第三彈性常數和第四彈性常數；

基于所述第一彈性常數、所述第二彈性常數、所述第三彈性常數、所述第四彈性常數和所述第五彈性常數，通過湯姆森彈性常數公式，計算各向異性參數。

4.根據權利要求1所述的計算頁巖各向異性參數的方法，其中，所述三維圖像通過對所述取芯方向與層理成θ角的頁巖巖芯的二維CT切片圖進行疊加組合，獲得所述取芯方向與層理成θ角的頁巖巖芯的三維圖像。

5.根據權利要求4所述的計算頁巖各向異性參數的方法，其中，所述疊加組合通過VGStudio軟件實現。

6.根據權利要求3所述的計算頁巖各向異性參數的方法，其中，所述頁巖彈性常數計算公式包括：

其中，C₁₁為所述第一彈性常數，C₃₃為所述第二彈性常數，C₆₆為所述第三彈性常數，C₄₄為所述第四彈性常數，V_p11為所述取芯方向垂直層理的頁巖巖芯的縱波速度，V_p33為所述取芯方向平行層理的頁巖巖芯的縱波速度，V_s66為所述取芯方向垂直層理的頁巖巖芯的橫波速度，V_s44為所述取芯方向平行層理的頁巖巖芯的橫波速度，ρ為所述頁巖巖芯的密度。

7.根據權利要求3所述的計算頁巖各向異性參數的方法，其中，所述湯姆森彈性常數公式包括：

其中，ε、γ，δ為所述湯姆森各向異性參數，C₁₁為所述第一彈性常數，C₃₃為所述第二彈性常數，C₆₆為所述第三彈性常數，C₄₄為所述第四彈性常數，C₁₃為所述第五彈性常數。

8.一種計算機可讀存儲介質，其上存儲有計算機程序，其特征在于，所述程序被處理器執行時實現以下步驟：

步驟1：獲得頁巖巖心的三維圖像；

步驟2：根據三維圖像中的層理方向，測量出取芯方向與層理的夾角參數θ；

步驟3：利用所述夾角參數θ計算頁巖各向異性參數。

9.根據權利要求8所述的計算機可讀存儲介質，其中，利用所述夾角參數θ計算頁巖各向異性參數包括利用所述夾角參數θ，通過修正公式，計算第五彈性常數：

其中，C₁₃為第五彈性常數，C₄₄為第四彈性常數，C₁₁為第一彈性常數，C₃₃為第二彈性常數，ρ為所述頁巖巖芯的密度，V_pθ為所述取芯方向與層理成θ角的頁巖巖芯的縱波速度，D為求解中間量。

10.根據權利要求8所述的計算機可讀存儲介質，其中，所述三維圖像通過對所述取芯方向與層理成θ角的頁巖巖芯的二維CT切片圖進行疊加組合，獲得所述取芯方向與層理成θ角的頁巖巖芯的三維圖像。