[發明專利]一種改進列聯表法的巖體均質區劃分方法

有效

• 申請號: 202110480033.0
• 申請日: 2021-04-30
• 公開（公告）號: CN113284225B
• 公開（公告）日: 2023-05-16
• 發明（設計）人: 劉鐵新;李嘉琪;姜諳男;梁正召;張正虎;詹必雄;韓鞠;唐衛平;侯拉平;張霄漢;周立飛;劉林濤
• 申請（專利權）人: 大連海事大學
• 主分類號: G06T17/00
• 分類號: G06T17/00;G06T7/187
• 代理公司: 大連東方專利代理有限責任公司 21212
• 代理人: 李洪福

1.一種改進列聯表法的巖體均質區劃分方法，其特征在于，包括：

分別記錄的兩個區域的節理的產狀、跡長數據，所述產狀包括傾向和傾角；

改進施密特圖建立新坐標系：所述施密特圖為半球面在水平面的投影，投影包括中心圓和以所述中心圓為圓心的多個圓環；多個所述圓環和所述中心圓以半徑由小至大的方式代表多個傾角帶；每個傾角帶內代表節理的傾角范圍；每個所述傾角帶內均等劃分為多個傾向塊；每個所述傾向塊代表節理的傾向范圍；在施密特圖上增加縱向垂直于所述施密特圖的跡長坐標軸，且跡長坐標軸用于顯示跡長；所述跡長坐標軸按照跡長劃分為n個跡長范圍區間；所述傾向塊、所述傾角帶和所述跡長范圍區間的交集構成胞元；

分別提取兩個區域的節理的產狀和跡長數據在所述新坐標系下投影得到兩個投影圖：相同傾角范圍、相同傾向范圍、相同跡長范圍區間的節理位于其所對應的胞元內；

以兩個所述投影圖為基礎建立列聯表：以所述胞元為列、所述投影圖數為行建立列聯表，所述胞元中的數值為所述胞元中節理的投影頻數；并根據所述列聯表得到卡方檢驗的卡方值；

確定卡方檢驗的顯著性水平和自由度；

根據卡方檢驗對兩個區域的節理進行均質區劃分。

2.根據權利要求1所述的一種改進列聯表法的巖體均質區劃分方法，其特征在于，所述卡方值χ²為：

式中：i為所述投影圖，j為胞元，c為胞元總數，f為胞元中節理的觀測頻數；e為胞元中節理的期望頻數，e＝(R_iC_j)/N，其中R為所述區域中的所有胞元中節理的總觀測頻數，C為兩個所述投影圖中同一胞元中的節理的總觀測頻數。

3.根據權利要求2所述的一種改進列聯表法的巖體均質區劃分方法，其特征在于，所述卡方檢驗的顯著性水平的計算公式為：

f(n)＝ae^-a(n+b)+d

式中，系數a＝0.8078，b＝2.4926，d＝0.0000055。

4.根據權利要求1所述的一種改進列聯表法的巖體均質區劃分方法，其特征在于，所述自由度為kn-1，其中k為所述傾向快的數量與兩個所述投影圖在同一縱列中胞元均為空的縱列數的差值。

5.根據權利要求1所述的一種改進列聯表法的巖體均質區劃分方法，其特征在于，采用卡方檢驗的p值法進行兩個區域的均質區判定。

6.根據權利要求1所述的一種改進列聯表法的巖體均質區劃分方法，其特征在于，分為三個傾角帶，內圈的傾角帶的傾角或等于0度，且小于28.1度；中間環的傾角帶的傾角大于或等于28.1度，且小于59度；外圈環的傾角帶的傾角大于或等于59度，且小于或等于90度；所述內圈的傾角帶按照順時針等分為四個傾向塊；所述中間換的傾角帶按照順時針等分為12個傾向塊；所述外圈的傾角帶按照順時針等分為18個傾向塊。