汪 偉，羅周全，秦亞光，姚 曙，顏克俊

深部開采初始地應(yīng)力場非線性反演新方法

汪 偉1，羅周全1，秦亞光1，姚 曙2，顏克俊2

(1.中南大學(xué)資源與安全工程學(xué)院，湖南長沙，410083；2.凡口鉛鋅礦，廣東韶關(guān)，512325)

針對凡口礦深部開采面臨的巖體脆?延性過渡、地壓活動頻繁等現(xiàn)狀，采用改進(jìn)型LUT套孔應(yīng)力測定系統(tǒng)開展礦山深部原巖應(yīng)力實測工作。對514組國內(nèi)地應(yīng)力實測值進(jìn)行擬合分析，綜合考慮復(fù)雜地質(zhì)體建模、反演載荷表達(dá)、載荷系數(shù)與應(yīng)力實測值之間的非線性映射，編寫非線性加載命令及載荷系數(shù)搜索程序，提出以混沌搜索理論為核心的地應(yīng)力場非線性反演新方法。將該方法應(yīng)用于凡口礦深部初始地應(yīng)力場反演。研究結(jié)果表明：應(yīng)力反演值與實測值平均擬合精度達(dá)90%；將非線性混沌方法用于地應(yīng)力場反演是可行的，能為深部安全開采提供有效支撐。

深部開采；初始地應(yīng)力場；非線性載荷；混沌搜索；反演分析

隨著社會經(jīng)濟(jì)發(fā)展對礦產(chǎn)需求量的大幅度增加，淺部資源日趨枯竭，許多礦山已進(jìn)入深部開采階段。深部開采面臨“三高一擾動”的特殊環(huán)境，工程及巖體力學(xué)響應(yīng)特征發(fā)生極大變化，地壓顯現(xiàn)更為明顯，嚴(yán)重威脅井下安全。大量的工程實踐與科學(xué)研究表明，地應(yīng)力是引起地下巖體工程變形和破壞的根本作用力[1]。因此，解決深部開采所面臨的各項技術(shù)難題，前提在于掌握地應(yīng)力場分布特征及變化規(guī)律。自1932年對美國胡夫大壩完成地應(yīng)力測量后，地應(yīng)力測量及相關(guān)研究迅速成為巖體力學(xué)的重要研究方向。目前，現(xiàn)場實測是獲取地應(yīng)力狀態(tài)最可靠的方法，然而，實測存在成本高、周期長、數(shù)據(jù)離散性大等不足，限制了其大規(guī)模開展，因此，基于有限地應(yīng)力實測資料而開展的地應(yīng)力場反演分析是當(dāng)前地應(yīng)力研究的重點。其中，以GUO等[2]提出的應(yīng)力多元回歸方法應(yīng)用最廣泛，眾多研究均以該方法為基礎(chǔ)。陳章華等[3]采用ANSYS建立甘肅省烏鞘嶺隧道數(shù)值模型，并設(shè)定線性的自重和構(gòu)造載荷，通過偏最小二乘回歸法反演初始地應(yīng)力場。趙德安等[4]在有限元數(shù)值模擬基礎(chǔ)上，編制應(yīng)力多元回歸程序MEBA，實現(xiàn)了地應(yīng)力場反演的拓展分析。金長宇等[5]通過真實還原地表剝蝕及斷層構(gòu)造運動，實現(xiàn)了對四川省白鶴灘廠區(qū)復(fù)雜地應(yīng)力場的多元回歸反演。袁海平等[6]構(gòu)建了地應(yīng)力場多元線性回歸模型，并用Fortran語言編寫了實現(xiàn)程序。謝紅強等[7]將多元回歸方法與灰色控制理論結(jié)合，實現(xiàn)了某水電站初始地應(yīng)力場的反演。景鋒等[8]針對地層變化趨勢，提出分次施加對應(yīng)載荷，提高了多元回歸方法的反演精度。應(yīng)力多元回歸方法原理簡單、適用性強，但存在一定的局限性：1)構(gòu)造應(yīng)力分布規(guī)律復(fù)雜，往往表現(xiàn)出典型的非線性特征。該方法對構(gòu)造應(yīng)力的反演是在數(shù)值模型邊界施加隨埋深線性變化的應(yīng)力或位移載荷，與其表現(xiàn)的非線性分布特征不符。2)假定自變量(載荷)與因變量(實測應(yīng)力)之間存在線性關(guān)系，在此基礎(chǔ)上求解自變量回歸系數(shù)，并未考慮載荷因素的非線性作用機理。鑒于應(yīng)力多元回歸方法的不足，一些學(xué)者開始引入遺傳算法、位移不連續(xù)法、粒子群算法等[9?11]，或構(gòu)建載荷因素與地應(yīng)力實測值之間的非線性聯(lián)系，或在數(shù)值模型邊界施加非線性變化的載荷條件，以期實現(xiàn)更準(zhǔn)確的初始地應(yīng)力場反演。從原理上而言，地應(yīng)力場非線性反演方法無疑更接近于實際情況。本文作者在借鑒已有非線性反演方法的基礎(chǔ)上，從復(fù)雜地質(zhì)體精確建模、反演載荷合理表達(dá)、載荷因素與地應(yīng)力實測值之間的非線性映射等方面綜合考慮，提出以混沌搜索理論為核心的地應(yīng)力場反演新方法，并以凡口鉛鋅礦深部地應(yīng)力場反演為例，以期為深部安全開采提供有效支撐。

1 深部原巖應(yīng)力實測

凡口鉛鋅礦是我國最大的鉛鋅采選企業(yè)，為實現(xiàn)年產(chǎn)18萬t金屬的目標(biāo)，礦山現(xiàn)已進(jìn)入深部主采階段。與上部中段相比，深部面臨巖體脆—延性過渡、地壓活動頻繁等問題，安全生產(chǎn)受到較大威脅。為此，于深部中段開展原巖應(yīng)力實測工作。

1.1 測量儀器

套孔應(yīng)力解除法和水壓致裂法是應(yīng)用最廣泛的2種地應(yīng)力測量方法，但水壓致裂法本質(zhì)上是二維應(yīng)力測量方法，其精度和可靠性比應(yīng)力解除法的低[12]。因此，采用套孔應(yīng)力解除法進(jìn)行地應(yīng)力實測，測量儀器為改進(jìn)型LUT套孔應(yīng)力測定系統(tǒng)。該系統(tǒng)在應(yīng)變計定位、溫度補償和數(shù)據(jù)采集等方面較國內(nèi)設(shè)備均有較大優(yōu)勢，特別是其應(yīng)變計粘結(jié)時間僅需2～3m in，與國內(nèi)應(yīng)變計20 h以上的粘結(jié)時間相比，顯著提高了測量結(jié)果的可靠性。圖1所示為LUT應(yīng)變計探頭的原理示意圖。

圖1 LUT應(yīng)變計探頭Fig.1 LUT strain-gauge detector

如圖1所示，探頭包含3組應(yīng)變花活塞，每個活塞含有2個應(yīng)變花，應(yīng)變花由2個相互正交垂直應(yīng)變片覆蓋組成，因此，應(yīng)變計探頭上共12個應(yīng)變片。3組應(yīng)變花活塞沿探頭周向分布，空間位置分別為270°，30°和150°，其上的4個電阻應(yīng)變片與探頭軸線的夾角分別為90°，45°，0°和135°。

1.2 測點布置

凡口礦井下已開拓至?750m埋深水平，但在?650m至?750m水平之間可采礦量較少。為了解深部主采區(qū)地應(yīng)力場特征，測點應(yīng)布置在?650m水平以上。經(jīng)分析，測點布置于?550m中段和?650m中段，測點信息見表1。

1.3 測量過程及結(jié)果

采用LUT系統(tǒng)測量地應(yīng)力的步驟見文獻(xiàn)[13]。其中最需要注意的2個技術(shù)環(huán)節(jié)為：一是應(yīng)變計的快速定位，因膠結(jié)劑固化時間短(3m in內(nèi))，故操作過程需緊密銜接，一步到位；二是必須套取完整巖芯，否則測量失敗。待取出巖芯冷卻，由儀器自動記錄其解除應(yīng)變。取不同測段同方向應(yīng)變(空間上有12個不同角度)的平均值為最終解除應(yīng)變。進(jìn)行雙軸率定實驗，獲取巖芯彈性參數(shù)，見表2。將巖芯的最終解除應(yīng)變和彈性參數(shù)輸入系統(tǒng)自帶的LUT-str應(yīng)力計算程序，得到凡口礦深部原巖應(yīng)力實測結(jié)果，見表3。

表1 地應(yīng)力測點信息Table1 Information of geostressmeasuring points

由表3可知：測點某2個主應(yīng)力軸傾向水平或與水平面夾角小于30°，另一主應(yīng)力軸接近于垂直或與水平面夾角大于60°，表現(xiàn)出典型的水平應(yīng)力場特征。測點處最大主應(yīng)力均超過25MPa(高應(yīng)力)，最大主應(yīng)力方向均為北北東向，平均方位角為29.50°(北29.50°東)。

2 地應(yīng)力場非線性反演新方法

2.1 復(fù)雜地質(zhì)體建模方法

礦山地質(zhì)體賦存條件多變、差異明顯，普通建模方法無法滿足其復(fù)雜建模的需求。為此，提出如下Surpac-Flac3D的復(fù)雜地質(zhì)體建模方法。

1)將生探地質(zhì)剖面預(yù)處理為Surpac線串文件，按多種三角網(wǎng)連接方法將地質(zhì)線串文件連成實體模型(dtm或3dm格式)，用于表達(dá)復(fù)雜的地層、斷層及礦體模型。建立Surpac塊體模型，用實體模型數(shù)據(jù)對其進(jìn)行約束賦值，得到經(jīng)六面體規(guī)則剖分的地質(zhì)塊體模型。

2)Surpac塊體數(shù)據(jù)格式與FLAC3D并不兼容，因此，提出如下數(shù)據(jù)格式轉(zhuǎn)換技術(shù)：

①通過Surpac將約束后的塊體模型導(dǎo)出為質(zhì)心文件(.csv文件)，采用Access數(shù)據(jù)庫打開塊體質(zhì)心文件，重新命名字段，去除主鍵，生成“導(dǎo)出表”。

②在Access中(SQL視圖)創(chuàng)建“轉(zhuǎn)換表”，以存儲經(jīng)坐標(biāo)運算后的塊體單元節(jié)點信息；將“導(dǎo)出表”中的單元數(shù)據(jù)由六面體質(zhì)心點轉(zhuǎn)換為FLAC3D中的P0～P3這4點坐標(biāo)，并將轉(zhuǎn)換結(jié)果寫入“轉(zhuǎn)換表”中?！稗D(zhuǎn)換表”創(chuàng)建及寫入信息命令見圖2。

圖2 數(shù)據(jù)格式轉(zhuǎn)換命令Fig.2 Transform wordsof data format

表2 巖芯最終解除應(yīng)變及彈性參數(shù)Table2 Final release strain and elastic parameters of coring-rocks

表3 主應(yīng)力大小及方向?qū)崪y結(jié)果Tab le 3 Measured resultsof principal stress value and direction

③將寫入單元節(jié)點信息的“轉(zhuǎn)換表”導(dǎo)出為.txt格式的文本文件。

④通過FLAC3D導(dǎo)入生成的文本文件，合并單元節(jié)點，獲得復(fù)雜地質(zhì)體的數(shù)值分析模型。

2.2 反演載荷表達(dá)

為便于分析，建模選用笛卡爾坐標(biāo)系，并將地質(zhì)體所受的復(fù)雜地應(yīng)力作用分解為如圖3所示的邊界應(yīng)力形式。

邊界應(yīng)力包括：1)x方向邊界水平擠壓構(gòu)造運動應(yīng)力σx；2)y方向邊界水平擠壓構(gòu)造運動應(yīng)力σy；3)自重應(yīng)力G；4)xz平面內(nèi)豎向剪切變形構(gòu)造運動應(yīng)力τzx；5)xy平面內(nèi)剪切變形構(gòu)造運動應(yīng)力τxy；6)yz平面內(nèi)豎向剪切變形構(gòu)造運動應(yīng)力τyz。基于獨立應(yīng)力張量唯一確定主應(yīng)力的原理[14]，通過對這6組應(yīng)力分量的最優(yōu)逼近，即可實現(xiàn)對主應(yīng)力大小和方向的反演。

為確定應(yīng)力分量的函數(shù)表達(dá)通式，搜集整理國內(nèi)近幾年的地應(yīng)力實測值，共計514組。測點平均埋深480.9m，最大埋深1 220.0m。測量方法以套孔應(yīng)力解除法和水壓致裂法為主(共496組)，其他方法有聲發(fā)射法、差應(yīng)變分析法和孔徑變形法等。采用Maltab-cftool工具箱擬合實測應(yīng)力分量隨深度h的變化關(guān)系。表4所示為擬合效果較高的幾組模型信息，應(yīng)力最優(yōu)擬合曲線見圖4～6。

圖3 邊界應(yīng)力分解示意圖Fig.3 Decomposition of boundary stress

表4 擬合效果較好的模型信息Table4 Information of optimal fittingmodels

圖4 最大水平應(yīng)力擬合曲線Fig.4 Fitting curveof themaximum horizontalstress

由表4可知：最大水平應(yīng)力σHmax和最小水平應(yīng)力項和h項系數(shù)均為負(fù)值。一般地，地應(yīng)力隨深度增加而增大(應(yīng)力異常區(qū)除外)，即地應(yīng)力與深度h呈正相關(guān)，應(yīng)力擬合曲線(圖4、圖5)也驗證了這一規(guī)律。而對2次項擬合時，h2項和h項系數(shù)均為正值，更符合地應(yīng)力與深度正相關(guān)的規(guī)律。應(yīng)力轉(zhuǎn)換公式為[15]

式中：σN和τN分別為法向應(yīng)力和剪應(yīng)力；l和m分別為σHmax和σHmin與不同邊界面夾角的余弦。

根據(jù)式(1)所示的應(yīng)力分量與σHmax和σHmin之間為線性關(guān)系，σx和σy隨埋深h變化的關(guān)系式為

圖5 最小水平應(yīng)力擬合曲線Fig.5 Fitting curve of them inimum horizontal stress

圖6 垂直應(yīng)力擬合曲線Fig.6 Fitting curveof vertical stress

σHmin采用非線性模型擬合時，各項精度指標(biāo)(誤差平方和、均方差及校正后決定系數(shù))均明顯比線性模型的高，3次項擬合精度更高。但采用3次項擬合時，h3

式中：a1，b1，c1，a，b和c為載荷系數(shù)。當(dāng)實測主應(yīng)力軸與x或y方向一致時，可按式(2)直接施加構(gòu)造應(yīng)力載荷，此時邊界剪應(yīng)力為0MPa。而當(dāng)實測主應(yīng)力軸與x或y方向不同時，根據(jù)正應(yīng)力與剪應(yīng)力的轉(zhuǎn)換關(guān)系，施加對應(yīng)剪應(yīng)力荷載。

對于式(2)中的非線性載荷，F(xiàn)LAC3D中未提供初始的加載命令，故通過軟件內(nèi)置的FISH語言，編寫如圖7所示的非線性加載命令。

垂直應(yīng)力采用線性擬合時的精度均最高，見表4。，特別是當(dāng)埋深小于700m時，應(yīng)力實測值密集分布于擬合曲線兩側(cè)(圖6)，擬合精度很高。采用線性擬合時，h項系數(shù)為0.025 4，對應(yīng)25.4MN/m3的線性變化梯度，符合三大類巖石容重的取值范圍，表明自重應(yīng)力與垂直應(yīng)力十分接近。故定義自重修正系數(shù)ag，在FLAC3D中通過“set0 0?10*ag”命令修正默認(rèn)自重梯度，以實現(xiàn)對自重應(yīng)力的反演。

圖7 非線性反演載荷編程實現(xiàn)語言Fig.7 Program w ordsof applying nonlinear inversion load

2.3 混沌反演方法

基于復(fù)雜地質(zhì)體建模方法和反演載荷構(gòu)造表達(dá)，提出如下地應(yīng)力場非線性反演新方法。

1)構(gòu)建自重修正系數(shù)ag與非線性載荷系數(shù)(a1，b1，c1，a，b，c)的組合樣本，據(jù)此開展初始地應(yīng)力場試算。提取實測點處應(yīng)力試算值，用載荷系數(shù)組合與應(yīng)力試算值的樣本對訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)，構(gòu)建二者之間的非線性映射關(guān)系。

2)輸入任意載荷系數(shù)組合p時，定義神經(jīng)網(wǎng)絡(luò)輸出(應(yīng)力試算值))(pijσ與應(yīng)力實測值ijσ逼近的目標(biāo)函數(shù)為

式中：i為測點編號；j為應(yīng)力分量編號。

3)考慮Logistic映射[16]：

式中：μ為控制參量。當(dāng)μ=4時，方程是[0,1]區(qū)間的滿映射，且進(jìn)入完全混沌狀態(tài)。

4)任意給定7個載荷系數(shù)[0,1]區(qū)間內(nèi)相異的初值，代入式(4)，得到初始混沌變量pi，并將pi“放大”至對應(yīng)載荷系數(shù)取值范圍[di,ei]：

式中：p′i為“放大”后的混沌變量；ri和si為放大系數(shù)；i=1,2,…,7。

5)將混沌變量p′i代入式(3)，計算每一步迭代的目標(biāo)函數(shù)值F(pi′)。設(shè)第1步迭代目標(biāo)函數(shù)值為，則置，繼續(xù)迭代直至保持極小，輸出對應(yīng)最優(yōu)解p′*。

i

6)通過Matlab將上述步驟編制成混沌搜索程序，根據(jù)尋優(yōu)搜索得到的載荷系數(shù)進(jìn)行數(shù)值正算，獲得反演初始地應(yīng)力場。

3 工程應(yīng)用

3.1 數(shù)值模型及材料參數(shù)

深部研究范圍包含：地層(D2db,D2da,D3ta)，鉛鋅礦體和黃鐵礦體以及主要的F3控礦斷層?？紤]邊界效應(yīng)，確定建模范圍(x，y，z)為(2 440～2 890，7 830～8 730，?750～?450)m，構(gòu)建地應(yīng)力場反演數(shù)值模型。圖8所示為控礦斷層與礦體的位置關(guān)系(深色為礦體，淺色為斷層)。斷層走向為NE15°～25°，傾向為105°～115°，傾角為65°～85°。

圖8 礦體與斷層復(fù)合模型Fig.8 Compoundmodelof orebody and fault

采集代表性巖樣，進(jìn)行室內(nèi)物理力學(xué)參數(shù)實驗，得到巖石物理力學(xué)參數(shù)。經(jīng)霍克?布朗準(zhǔn)則對其折減處理，得到深部巖體物理力學(xué)參數(shù)，見表5。

表5 深部巖體力學(xué)強度參數(shù)Table5 Mechanical parametersof deep rockmass

3.2 混沌搜索反演載荷系數(shù)

構(gòu)建載荷系數(shù)(ag，a1，b1，c1，a，b，c)的典型樣本組合，按FISH編程命令施加非線性構(gòu)造載荷，采用“set0 0?10*ag”命令設(shè)定自重梯度，開展初始地應(yīng)力場數(shù)值試算。提取載荷系數(shù)(輸入)和實測點應(yīng)力試算值(輸出)，通過二者的樣本對訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)，見圖9。建立載荷系數(shù)到應(yīng)力試算值的非線性映射，并用于代替數(shù)值試算過程。

圖9 載荷系數(shù)到應(yīng)力試算值的BP神經(jīng)網(wǎng)絡(luò)映射Fig.9 BPnetw ork mapping between load parametersand simulated stress values

至此，問題轉(zhuǎn)換為求解滿足神經(jīng)網(wǎng)絡(luò)輸出(應(yīng)力試算值)與應(yīng)力實測值最優(yōu)逼近的載荷系數(shù)條件。目標(biāo)函數(shù)m in F(p)為非凸函數(shù)，線性或隨機搜索方法在解決此類尋優(yōu)問題時，容易陷入局部極值[17]?；煦缦到y(tǒng)具有典型的遍歷性特點，可在確定邊界內(nèi)不重復(fù)自身軌跡地運動，故引入混沌載波進(jìn)行載荷系數(shù)全局尋優(yōu)。根據(jù)初始地應(yīng)力場試算結(jié)果，定義載荷系數(shù)(ag，a，b，c，a1，b1，c1)的可行解區(qū)間分別為[0.500,1.500]，[?0.001,0.001]，[?0.100,0.100]，[?10.000,10.000]，[?0.001,0.001]，[?0.100,0.100]，[?10.000,10.000]。采用編寫的Matlab搜索程序進(jìn)行載荷系數(shù)尋優(yōu)，混沌變量載波及搜索過程見圖10和圖11。圖10中μ為控制參數(shù)，x為取值區(qū)間。

由圖11可知：搜索過程歷經(jīng)400步迭代后逐漸趨于穩(wěn)定；500步后，搜索平穩(wěn)收斂，輸出尋優(yōu)載荷系數(shù)為：ag=1.012 000，a1=0.000 045，b1=0.012 100，c1=0.490 000，a=0.000 035，b=0.011 700，c=0.280 000。

3.3 反演初始地應(yīng)力場分析

將經(jīng)混沌搜索獲得的最優(yōu)載荷系數(shù)代入式(2)，得到非線性反演載荷的函數(shù)表達(dá)式為

圖10 Logistic混沌載波Fig.10 Chaotic transformation by logistic formula

圖11 反演載荷系數(shù)混沌搜索過程Fig.11 Chaotic search processof inversion load parameters

在反演模型中導(dǎo)入編寫的Fish命令流(見圖7)，施加式(6)對應(yīng)的非線性載荷，按“setgrav 0 0?10.012”命令設(shè)定自重梯度。據(jù)此開展數(shù)值正算反演，生成初始地應(yīng)力場，反演結(jié)果見表6和表7。

表6 主應(yīng)力值與方向反演結(jié)果Tab le 6 Inversion result of principal stress value and direction

表7 應(yīng)力分量反演結(jié)果Table7 Inversion resultof stress com ponents MPa

判斷地應(yīng)力場反演效果的主要依據(jù)是計算平衡后模型內(nèi)的主應(yīng)力數(shù)值和方向[18]。由表6可知：反演主應(yīng)力與實測主應(yīng)力在數(shù)值和方向上都很接近，平均擬合精度達(dá)到90%。表7中，應(yīng)力分量的反演平均誤差依次為0.89，1.82和1.08MPa，其中，σx，σy和σz的平均反演精度為97.11%，98.26%和97.59%。整體而言，反演值與實測值擬合度高，反演效果較好。但是，剪應(yīng)力τxy，τyz和τzx的反演平均誤差依次為1.29，2.23和2.57MPa，波動性較大。究其原因，一方面與礦山開采年代久遠(yuǎn)、深部地應(yīng)力場變化復(fù)雜有關(guān)；另一方面也反映了地應(yīng)力場是受多場多因素耦合作用的復(fù)雜系統(tǒng)。在當(dāng)前技術(shù)條件下，反演研究只能從主要影響因素(自重和構(gòu)造作用)入手，根據(jù)有限的“點吻合”，實現(xiàn)宏觀意義上的“場吻合”。可見，反演初始地應(yīng)力場是合理的。

圖12 ?600m中段最大主應(yīng)力分布(Z=?600m)Fig.12 Themaximum principal stress distribution on?600m level

?600m中段是礦山深部主采中段，高應(yīng)力勢必會對安全開采產(chǎn)生重大影響。?600m水平的主應(yīng)力分布特征如圖12和圖13所示。從圖12和圖13可見：主采中段地應(yīng)力場服從層狀分布，垂直應(yīng)力和水平構(gòu)造應(yīng)力均超過20MPa，且以南北向水平構(gòu)造應(yīng)力起主要控制作用；在F3控礦斷層內(nèi)，因巖體強度低，賦存應(yīng)力顯著減小，形成了沿走向長度近700m的應(yīng)力弱化帶。斷層兩側(cè)產(chǎn)生應(yīng)力集中，區(qū)域內(nèi)最大主應(yīng)力可達(dá)36MPa，最小主應(yīng)力可達(dá)14MPa，影響范圍自?600m S5#S礦房采場至N6-7#礦柱采場(走向方向Y為8 236～850 8m)。高應(yīng)力集中是誘發(fā)深部地壓災(zāi)害的重要原因，在該區(qū)域作業(yè)時應(yīng)予以重視。

圖13 ?600m中段最小主應(yīng)力分布(Z=?600m)Fig.13 Them inimum principal stress distribution on?600m level

4 結(jié)論

1)采用改進(jìn)型LUT套孔應(yīng)力測定系統(tǒng)于凡口礦?550m和?650m中段開展原巖應(yīng)力實測，獲得礦山深部實測原巖應(yīng)力分布特征為：測點處某2個主應(yīng)力軸傾向水平或與水平面夾角小于30°，另一主應(yīng)力軸接近于垂直或與水平面夾角大于60°，表現(xiàn)出典型的水平應(yīng)力場特征；測點處最大主應(yīng)力均超過25MPa，最大主應(yīng)力方向均為北北東向，平均方位角為29.50°(北29.50°東)；測點處最大水平主應(yīng)力與垂直應(yīng)力之比(側(cè)壓系數(shù))在1.16～1.49之間，與我國大陸區(qū)域側(cè)壓系數(shù)分布規(guī)律相符。

2)提出以混沌搜索理論為核心的地應(yīng)力場非線性反演方法。將該方法應(yīng)用于凡口礦深部初始地應(yīng)力場反演，反演主應(yīng)力與實測主應(yīng)力平均擬合精度達(dá)到90%，表明地應(yīng)力場非線性混沌反演方法是合理、可行的。

3)基于反演初始地應(yīng)力場，獲得凡口礦?600m主采中段地應(yīng)力場分布特征：垂直應(yīng)力和水平構(gòu)造應(yīng)力均超過20MPa，且以南北向水平構(gòu)造應(yīng)力起主要控制作用。F3控礦斷層內(nèi)，因巖體強度低，賦存應(yīng)力顯著減小，形成了沿走向長度近700m的應(yīng)力弱化帶。斷層兩側(cè)產(chǎn)生應(yīng)力集中，區(qū)域內(nèi)最大主應(yīng)力可達(dá)36MPa，最小主應(yīng)力可達(dá)14MPa，影響范圍自?600m S5#S礦房采場至N6-7#礦柱采場(走向方向Y范圍為8 236～8 508m)。

[1]OUYANG Zhenhua,LI Changhong,XU Wancai,et al. Measurements of in-situ stress and mining-induced stress in Beiminghe iron mine of China[J].Journal of Central South University of Technology,2009,16(1):85?90.

[2]GUO Huaizhi,MA Qichao,XUE Xicheng,et al.The analysis method of the initial stress field for rock mass[J].Chinese Journalof Geotechnical Engineering,1983,5(3):64?75.

[3]陳章華,陳磊,紀(jì)洪廣.基于偏最小二乘法的地應(yīng)力場擬合[J].北京科技大學(xué)學(xué)報,2013,35(1):1?7.

CHEN Zhanghua,CHEN Lei,JIHongguang.Geo-stress field fitting based on the partial least square method[J].Journal of University of Science and Technology Beijing,2013,35(1):1?7.

[4]趙德安,李國良,陳志敏,等.烏鞘嶺隧道三維地應(yīng)力場多元有限元回歸拓展分析[J].巖石力學(xué)與工程學(xué)報,2009,28(S1): 2687?2694.

ZHAO Dean,LI Guoliang,CHEN Zhimin,et al.Threedimensional FE regression analysis of multivariate geostress field of Wushaoling tunnel[J].Chinese Journal of Rock Mechanicsand Engineering,2009,28(S1):2687?2694.

[5]金長宇,馮夏庭,張春生.白鶴灘水電站初始地應(yīng)力場研究分析[J].巖土力學(xué),2010,31(3):845?850,855.

JIN Changyu,FENG Xiating,ZHANG Chunsheng.Research on initial stress field of Baihetan hydropower station[J].Rock and SoilM echanics,2010,31(3):845?850,855.

[6]袁海平,趙科,王繼倫,等.多元線性回歸地應(yīng)力場反演方法研究[J].有色金屬科學(xué)與工程,2011,2(1):59?62.

YUAN Haiping,ZHAO Ke,WANG Jilun,et al.On themethod of multivariate linear regression for geo-stress field[J]. NonferrousM etals Science and Engineering,2011,2(1):59?62.

[7]謝紅強,何江達(dá),肖明礫.大型水電站廠區(qū)三維地應(yīng)力場回歸反演分析[J].巖土力學(xué),2009,30(8):2471?2476.

XIE Hongqiang,HE Jiangda,XIAO M ingli.Regression analysis of 3D initial geostress in region of underground powerhouse for large hydropower station[J].Rock and Soil Mechanics,2009, 30(8):2471?2476.

[8]景鋒,尹健民,陳昊,等.考慮復(fù)雜地質(zhì)成因的工程區(qū)地應(yīng)力場反演[J].長江科學(xué)院院報,2009,26(9):51?54.

JING Feng,YIN Jianmin,CHEN Hao,etal.Regression analysis of initial geostress field considering complex geologic origin[J]. Journal of Yangtze River Scientific Research Institute,2009, 26(9):51?54.

[9]GE Hongwei,LIANG Yanchun,LIU Wei.Applications of artificial neural networks and genetic algorithms to rock mechanics[J].Rock Mechanics and Engineering,2004,23(9): 1542?1550.

[10]LI Ke,WANG Yingyi,HUANG Xingchun.DDM regression analysis of the in-situ stress field in a non-linear fault zone[J]. International Journal of M inerals,M etallurgy and M aterials, 2012,19(7):567?573.

[11]賈善坡,伍國軍,陳衛(wèi)忠.基于粒子群算法與混合罰函數(shù)法的有限元優(yōu)化反演模型及應(yīng)用[J].巖土力學(xué),2011,32(S2): 598?603.

JIA Shanpo,WU Guojun,CHENWeizhong.Application of finite element inverse model based on improved particle swarm optimization and mixed penalty function[J].Rock and Soil M echanics,2011,32(S2):598?603.

[12]蔡美峰,郭奇峰,李遠(yuǎn),等.平煤十礦地應(yīng)力測量及其應(yīng)用[J].北京科技大學(xué)學(xué)報,2013,35(11):1399?1406.

CAI Meifeng,GUO Qifeng,LI Yuan,et al.In situ stress measurement and its application in the 10th M ine of Pingdingshan coalgroup[J].Journalof University of Science and Technology Beijing,2013,35(11):1399?1406.

[13]陳楓,饒秋華,徐紀(jì)成,等.應(yīng)變解除法原理及其在大紅山鐵礦地應(yīng)力測量中的應(yīng)用[J].中南大學(xué)學(xué)報(自然科學(xué)版),2007, 38(3):545?550.

CHEN Feng,RAO Qiuhua,XU Jicheng,etal.Principle of strain release method and its application to the in-situ stress measurement in Dahongshan Iron M ine[J].Journal of Central South University(Science and Technology),2007,38(3): 545?550.

[14]YUE Xiaolei,LI Yong,WANG Haipeng,et al.Inversion of initial in situ stress field by themethod ofmultivariate analysis and engineering application[J].App lied Mechanics and M aterials, 2011(44/45/46/47):1203?1206.

[15]袁風(fēng)波.巖體地應(yīng)力場的一種非線性反演新方法研究[D].武漢:中國科學(xué)院武漢巖土力學(xué)研究所,2007:50?52.

YUAN Fengbo.Study on a new nonlinear back analysismethod of geostress field[D].Wuhan:Geotechnical Mechanics Institute of Chinese Academy of Sciences,2007:50?52.

[16]CHEN Guanrong,MAO Yaobin,CHUIC K.A symmetric image encryption scheme based on 3D chaotic cat maps[J].Chaos, Solitons&Fractals,2004,21(3):749?761.

[17]LIXiang,LIU Guangying,QIJianxun.Fuzzy neural and chaotic searching hybrid algorithm and its application in electric customers’credit risk evaluation[J].Journal of Central South University,2007,14(1):140?143.

[18]YIDa,XU Mingyi,CHEN Shenghong,et al.Applications of artificial neural network to back analysis of initial stress field of rock masses[J].Rock and Soil M echanics,2004,25(6): 943?945.

(編輯 陳燦華)

A new nonlinear inversionmethod of geostress field in deepm ining

WANGWei1,LUO Zhouquan1,QINYaguang1,YAOShu2,YANKejun2

(1.SchoolofResourcesand Safety Engineering,Central South University,Changsha410083,China; 2.Fankou Lead?Zinc M ine,Shaoguan 512325,China)

Aim ing at rock brittle-ductile transition and frequent underground pressure activity in deep Fankou mine, geostressmeasurementw as carried out by using developed LUT over-coring system in deep level.Considering com plex geological modeling,inversion load expression and nonlinear mapping between load factors and measured stress,the nonlinear loading words and load parameters searching program were compiled,and thus a new nonlinear inversion method of geostress field characterized by chaotic search theorywas put forward.Thenew inversionmethod wasapplied in the Fankou m ine.The results show that the inversed stress fits w ell w ith measured stress,and the average fitting accuracy exceeds 90%.The nonlinear chaotic method is reasonable and workable for geostress field inversion,w hich can lay an effective foundation for safe production in deepm ining level.

deepm ining;in-situ geostress field;nonlinear load;chaotic search;inversion analysis

TD311

A

1672?7207(2017)03?0804?09

10.11817/j.issn.1672-7207.2017.03.031

2016?04?10；

2016?06?22

國家自然科學(xué)基金資助項目(51274250)；中南大學(xué)研究生自主探索創(chuàng)新項目(2016zzts091)(Project(51274250)supported by the National Natural Science Foundation of China;Project(2016zzts091)supported by Independent Exploration and Innovation Foundation of CSU)

羅周全，教授，博士研究生導(dǎo)師，從事金屬礦床開采及災(zāi)害監(jiān)控理論與技術(shù)研究；E-mail:lzq501505@163.com