[發明專利]一種地表巖體移動參數的獲取方法

在審

• 申請號: 201710236516.X
• 申請日: 2017-04-12
• 公開（公告）號: CN107014683A
• 公開（公告）日: 2017-08-04
• 發明（設計）人: 王永寶;張宗濤;劉鵬;程布康;任康;趙志華;馮上運;郭純建
• 申請（專利權）人: 山東東山王樓煤礦有限公司
• 主分類號: G01N3/08
• 分類號: G01N3/08;G06F17/50
• 代理公司: 濟南舜源專利事務所有限公司37205
• 代理人: 張亮

技術領域

本發明涉及地表巖體領域，尤其涉及一種地表巖體移動參數的獲取方法。

背景技術

計算地表移動的常用方法為概率積分法、典型曲線法、剖面函數法等。其中概率積分法使用高斯曲線作為開采的連續影響曲線，配合使用概率積分函數計算，使該方法比較嚴密符合實際，但是使用方法較為繁瑣，工作量大。而典型曲線法將同類型地質采礦條件地表下沉盆地的移動和變形分布用無因次曲線和表格表示，但是這種方法要求應用于類型相同的地質采煤條件，有局限性。剖面函數法是根據地表下沉盆地剖面形狀來選擇描述下沉盆地剖面的響應函數作為地表移動和變形的公式，這種方法建立在整理大量實測資料的基礎上，用數理統計方法確定參數。因此如何能夠既依賴于工程地質和巖石力學理論，又依托巖體工程的現場實際量測得到地表巖體移動參數是丞待解決的技術問題。

發明內容

為了克服上述現有技術中的不足，本發明提供一種地表巖體移動參數的獲取方法，方法包括：

在二維和三維地應力場分析中，回歸元素都定為三個，把地應力回歸計算值作為因變量，把數值計算求得的自重應力場和地質構造應力場相應于實測點的應力計算值作為自變量，則回歸方程的形式為：

其中，k為觀測點的序號；為第k觀測點的回歸計算值；L_i為相應于自變量的多元回歸系數；為相應應力分量計算值的單列矩陣；

其中坐標軸為：對二維地應力場分析，Y軸為計算剖面走向；對三維地應力場分析，X軸和Y軸為計算域水平邊界方向，z軸為鉛垂向上方向；

對每一個應力狀態可確定一個回歸計算值觀測值與回歸方程的偏離程度，用該觀測值與回歸計算值之差，即該觀測的殘差表示全部觀測與回歸方程的偏離程度，用全部觀測值與回歸計算值的殘差平方和表示，即：

其中，k為觀測點序號；n為觀測點的個數，n＝n₁+n₂,n₁為三維地應力測點(套鉆孔應力解除法測量)的個數，n₂為二維地應力測點(水壓致裂法測量)的個數，二維應力測點編列在三維應力測點后面；j為應力分量序號；m為應力分量的個數；為第k觀測點序號為j應力分量實測值，例如對二維地應力場分析，當k≤n₁時，m＝3，當k>n₁時，m＝1；對三維地應力場分析，當k≤n₁時,m＝6，當k>n₁時，m＝3；采用數理統計的最小二乘法原理，使殘差平方和為最小，對L_i取偏導并令其為0，得到求解回歸系數L_i的正規方程組：

解得回歸系數(L_i＝1～3)后，就可根據回歸方程和數值計算結果求得計算域內各處的應力回歸計算值；應力回歸的效果依照回歸方程和回歸元素的顯著性檢驗值F和F_i進行檢驗：

其中，R_v和R分別為回歸差平方和與殘差平方和；N為回歸差平方和的自由度；V_i為回歸方程中各自變量的貢獻，用變量的偏回歸差平方和表示：其中，C_ii為正規方程組系數矩陣逆陣的主元素。

從以上技術方案可以看出，本發明具有以下優點：

地表巖體移動參數的獲取方法適用于水電站地下廠房的建造和礦山開采中。地表巖體移動參數的獲取方法還應用于基坑工程、邊坡、大壩等各個領域。地表巖體移動參數的獲取方法的計算精度較高，對地表沉降預測有較強的實用性。

地表巖體移動參數的獲取方法是解決融合了數值方法與巖土工程，地表巖體移動參數的獲取方法既依賴于工程地質和巖石力學理論，又依托巖體工程的現場實際量測，地表巖體移動參數的獲取方法具有理論性和實踐性的實用技術。將地表巖體移動參數的獲取方法應用于地表巖移參數計算中，不但能形象的反應巖層和地表移動的力學性質，還能通過計算得到相對準確的地表移動參數。

具體實施方式