王明柱，郭廣禮，王 磊，王 彬

（1.中國礦業(yè)大學江蘇省資源環(huán)境信息工程重點實驗室，江蘇徐州221116;

2.中國礦業(yè)大學國土環(huán)境與災害監(jiān)測國家測繪局重點實驗室，江蘇徐州 221116）

基于嶺估計的概率積分法預計參數的求取

王明柱1，2，郭廣禮1，2，王 磊1，2，王 彬1，2

（1.中國礦業(yè)大學江蘇省資源環(huán)境信息工程重點實驗室，江蘇徐州221116;

2.中國礦業(yè)大學國土環(huán)境與災害監(jiān)測國家測繪局重點實驗室，江蘇徐州 221116）

采用最小二乘擬合求取概率積分法預計參數時，經常遇到發(fā)散問題，致使求參結果失真。而求參過程中法矩陣的病態(tài)是造成結果發(fā)散的主要原因之一。將嶺估計應用于預計參數的求取中，能夠有效地克服傳統最小二乘擬合求參時法矩陣的病態(tài)。在給出基于嶺估計求參模型的同時，使用Matlab實現了參數的自動求取。最后通過模擬算例驗證了該求參模型克服法矩陣病態(tài)的有效性及求參結果的可靠性。

嶺估計;概率積分法;求參;開采沉陷;L曲線法

概率積分法是進行開采沉陷預計的主要方法，預計過程中預計參數的準確性決定了地表移動變形預計的精度。參數的影響因素主要有:覆巖的巖性和地層結構、采深、采厚等［1］。我國大部分礦區(qū)都根據所設的地表移動觀測站反演出了適合自己礦區(qū)的概率積分預計參數。

預計參數的求取通常采用最小二乘擬合的方法，該方法計算簡單，且求參精度能滿足工程需求。而最小二乘擬合求參過程中法矩陣一旦病態(tài)，就致使迭代結果發(fā)散、失真，造成求參困難。為解決由于法方程系數陣病態(tài)而導致最小二乘估計不穩(wěn)定問題，統計學家們提出了嶺估計、廣義嶺估計、主成分估計等有偏估計來改善這種情況下的最小二乘估計。嶺估計法從減小均方誤差的角度出發(fā)，最初由A.E.Hoerl于1962年提出，并于1970年由他和R.W.Kennard對該方法做了系統的發(fā)展。該方法成為目前使用最多的一種有偏估計［2］。

采用嶺估計可有效克服病態(tài)矩陣干擾，使求參結果更可靠，確保概率積分法開采沉陷預計精度。

1 嶺估計

1.1 嶺估計基本原理

設有觀測方程:

式中，L是觀測值;A是系數陣;X是待估計參數的真值;Δ是噪聲，Δ～N（0，δ20I）。則上式常規(guī)的最小二乘解為

若法矩陣N=ATPA條件數很大，即為病態(tài)矩陣時，求逆就會出現不穩(wěn)定，致使最小二乘估計的方法不可能得到未知參數的精確估值［3］。嶺估計法可解決這類問題［4］。其基本思想是利用原最小二乘估計的數學模型，在其法方程系數矩陣N的對角線上加上一個很小的正數來抑制法方程的病態(tài)，是一種在均方誤差意義下優(yōu)于最小二乘估計的有偏估計方法，是最小二乘估計的線性組合。的嶺估計解為［2］

式中，k為嶺參數。可見，嶺估計結果與嶺參數k有關，選擇不同的嶺參數，估計結果可能差異很大，所以嶺估計的關鍵是選擇合適的嶺參數。

1.2 嶺參數的確定

關于嶺參數的確定，目前已經有不少的方法，如嶺跡法、L曲線法、兩步解法及GCV法等，但是到目前為止還沒有一種公認的好方法。相比較而言，L曲線法易于確定嶺參數，是一種確定嶺參數的良好方法［5］。本文使用L曲線法來確定嶺參數，在此進行簡述，其他方法可參見文獻 ［5－7］。

根據正則化原理，（1）式的嶺估計準則為

求上式最大值就可得到最大曲率KLmax，其對應點即為所求點。然后求出該點所對應的嶺參數k，即為 （1）式所需要的k。具體求解過程可參見文獻 ［8］。

2 概率積分法預計參數的嶺估計模型

2.1 求參的嶺估計模型

概率積分法預計參數可由地表點實測移動變形資料求出。當布置的測點為正規(guī)觀測站時，可采用最小二乘曲線擬合法;當測點為非正規(guī)觀測站即一系列散點時，可采用最小二乘曲面擬合法。若最小二乘的法矩陣病態(tài)，則易造成求參結果失真、不可靠?；趲X估計的曲面擬合法不僅可以有效地抵御法矩陣病態(tài)的干擾，還可適用于非正規(guī)觀測站和殘缺觀測站資料的求參。

概率積分法的預計參數有下沉系數q、主要影響角正切tanβ、主要影響傳播角θ、拐點偏移距S1，S2，S3，S4及水平移動系數b。其任一點移動變形預計具體計算公式可參見文獻 ［1］。進行基于嶺估計的曲面擬合法求參時，首先利用實測地表下沉值來求取q，tanβ，θ，S1，S2，S3，S4這7 個參數。然后，認為這些參數已知，再利用地表點的水平移動觀測值進行迭代擬合求取b。

根據任一點移動變形預計公式，任一點下沉Wi可表示為任一點坐標 （xi，yi）和參數q，tanβ，S1，S2，S3，S4和 cotθ的函數，即

由（4），（5）式確定嶺參數k，則7個參數改正數的嶺估計解為

接下來基于嶺估計的迭代求參過程同上述由實測地表下沉值求參，在此不贅述。

2.2 求參結果的精度評定

所求取參數的效果如何，可通各參數的均方誤差和各點擬合中誤差來考查。

2.2.1 計算各參數的均方誤差

2.2.2 計算各點擬合中誤差

根據上述的求參嶺估計模型，使用Matlab程序實現了各參數的自動求取及參數的精度評定。

3 算例分析

通過模擬算例來檢驗該模型抵御法矩陣病態(tài)的效果及求參的可靠性。

模擬開采一個近水平煤層，采厚為3m，平均采深為300m，工作面詳情見表1。工作面上方布設2條地表移動觀測線 （圖1）:一條走向觀測線，一條與走向成45°夾角的觀測線。共設有觀測點94個，點間距約20m。確定 7個參數 q，tanβ，θ，S1，S2，S3，S4（表2），用概率積分任一點下沉預計公式計算出各測點下沉值，然后在各下沉值中加入標準差為10mm的隨機誤差，以此作為下沉觀測值進行參數反演。

表1 工作面信息

表2 預計參數真值

圖1 觀測線布設

表3 基于最小二乘估計、嶺估計的求參結果比較

最小二乘估計法第1次迭代計算時法矩陣N極度病態(tài)，其條件數數量級達1020，經第1次計算后部分參數已很大程度上偏離真值。繼續(xù)迭代計算后，各參數仍不收斂，最終計算結果發(fā)散。而嶺估計法大大改善了法矩陣的病態(tài)，第1次迭代計算時，其條件數數量級降為105。且迭代計算過程中，隨著各參數收斂于其真值，法矩陣條件數也愈小。

將基于嶺估計求得的各預計參數與表2中的真值相比，其精度足夠滿足預計工作的需要。最終各點下沉值的擬合中誤差為±10.3mm，接近于在原真值基礎上所加入的標準差為10mm的高斯白噪聲。驗證了該模型求參的可靠性。

4 結論

（1）在使用最小二乘估計求取概率積分法預計參數時，存在由于法矩陣的極度病態(tài)，致使求參結果發(fā)散的問題。

（2）嶺估計法能夠很好地改善法矩陣的病態(tài)性，較最小二乘估計其難點僅在于嶺參數的確定。所采用的L曲線法能夠方便準確地求出嶺參數，使嶺估計模型抗病態(tài)效果顯著。

（3）基于嶺估計的概率積分法求參模型，具有較好地抵御矩陣的病態(tài)性能，使求參結果更可靠。所建立的模型適用于全盆地的擬合求參，能夠處理不規(guī)則以及不完整的觀測站資料。

［1］何國清，楊 倫，凌賡娣，等.礦山開采沉陷學 ［M］.徐州:中國礦業(yè)大學，1991.

［2］崔希璋，於宗儔，陶本藻，等.廣義測量平差 （第二版）［M］.武漢:武漢大學出版社，2009.

［3］黃幼才.嶺估計及其應用［J］.武漢大學學報 （信息科學版），1987，12（4）:64－73.

［4］Hoerl A E，Kennard R W.Ridge Regression:Biased Estimation for Non－orthogonal Problems［J］ .Technometrics，1970，12（1）.

［5］汪明瑾，王靜龍.嶺回歸中確定K值的一種方法［J］.應用概率統計，2001，17（1）:7－13.

［6］黃海蘭，牛 犇.嶺參數確定的研究［J］.測繪科學，2011，36（4）:31－32.

［7］王振杰，歐吉坤.用L－曲線法確定嶺估計中的嶺參數 ［J］.武漢大學學報 （信息科學版），2004，29（3）:235－238.

［8］ Hansen P C，O＇Leary D P.The Use of The L－curve in The Regularization of Discrete Ill－posed Problems［J］ .SIAMJ.Sci.Com put，1993，14（6）:1487－1503.

［9］中華人民共和國煤炭工業(yè)部.建筑物、水體、鐵路及要井巷煤柱留設與壓煤開采規(guī)程 ［M］.北京:煤炭工業(yè)出版社，1986.

［10］張翠英.基于MapX的煤礦開采沉陷預計及地表移動觀測站數據處理系統的研發(fā)［D］.淮南:安徽理工大學，2009.

Prediction Parameters Calculation of Probability Integration Method Based on Ridge Estimate

WANG Ming-zhu1，2，GUO Guang-li1，2，WANG Lei1，2，WANG Bin1，2
（1.Key Laboratory of Resources Environment Information Engineering of Jiangsu Province，China University of Mining ＆ Technology，Xuzhou 221116，China;2.Key Laboratory of Land Environment and Disaster Monitoring of State Bureau of Surveying and Mapping，China University of Mining＆Technology，Xuzhou 221116，China）

Divergence problem will often result in parameters distortion in solve prediction parameters of probability integration method with least square method.Morbidity of normal matrix is one of main divergence reasons.Ridge estimate was applied in solving prediction parameters，which could effectively overcome distortion of normal matrix.At the same time，Matlab was used to realize parameters＇automatic solving.A simulation example was present to verify the model.

ridge estimate;probability integration method;parameter solve;mining subsidence;Lcurve method

TD325.2

A

1006-6225（2012）02-0017-03

2011－11－05

國家自然科學基金重點項目 （50834004）;國家公益性行業(yè)基金資助項目 （200811050）;中國礦業(yè)大學青年科研基金資助項目（09091252）

王明柱 （1986－），男，江蘇豐縣人，碩士研究生，主要研究方向為開采沉陷與“三下”采煤。

［責任編輯:王興庫］

基礎研究