黃 婧 劉 亮

(1.天津港東疆建設開發(fā)有限公司，天津 300463； 2.交通運輸部天津水運工程科學研究所，天津 300456；3.天津水運工程勘察設計院 天津市水運工程測繪技術企業(yè)重點實驗室，天津 300456)

0 引言

地面沉降對社會、經(jīng)濟和民生具有嚴重的不可逆轉的危害，它初期并不容易被發(fā)現(xiàn)，但是一旦發(fā)現(xiàn)則表明地面沉降給該地區(qū)的基礎設施、建筑物等帶來了極大的損壞。在早期，人們關注點主要是對地面沉降的補救和防治工作，往往忽略對地面沉降的預測研究。只有根據(jù)已有數(shù)據(jù)提前預知該地區(qū)未來的地面沉降趨勢情況，才能對該地區(qū)資源開采和基礎建設等做出輔助指導，才能更加有針對性的進行地面沉降防治工作。

1 概述

天津濱海新區(qū)是我國地面沉降問題十分嚴重的地區(qū)之一。據(jù)統(tǒng)計分析[1-3]，從20世紀50年代擁有有效記錄以及水準點監(jiān)測以來，到2008年，濱海新區(qū)大部分地區(qū)的地面沉降已經(jīng)超過1 m，不少于180 km2的區(qū)域的累計地面沉降量超過2.5 m，其最大累計地面沉降量已經(jīng)超過3 m，甚至部分區(qū)域由于地面沉降使得其地面標高已經(jīng)低于海平面。天津港位于濱海新區(qū)東部，由于大部分陸地區(qū)域是在海河入?？诘哪酁┥洗堤钤礻懚桑以诮畮啄赀M行了大量基礎設施建設，其地面沉降問題非常嚴峻。天津港地面沉降現(xiàn)狀在文獻[4]中已詳細闡述，本文在其研究基礎上進行沉降預測分析。

2 沉降預測模型

國內(nèi)外進行地面沉降的預測研究所用的模型按原理可以分為三類：數(shù)理統(tǒng)計與隨機統(tǒng)計模型、水土模型以及人工智能模型[5-7]。水土模型和人工智能模型理論復雜、不容易實現(xiàn)。數(shù)理統(tǒng)計與隨機統(tǒng)計模型僅僅需要已有的歷史數(shù)據(jù)資料，而無需其他相關的信息(如地下水流失數(shù)據(jù)、土質結構、工程施工等信息)，就能建立起預測模型，是一種較為簡單的模型方法，但是卻對數(shù)據(jù)擬合和預測有著十分有效的方法。地面沉降的數(shù)理統(tǒng)計模型包括時間序列模型、回歸分析模型、Verhulst、灰色理論模型等[8,9]。

2.1 模型選取

已有2002年—2018年天津港基礎控制網(wǎng)水準點的外業(yè)實測數(shù)據(jù)，未收集到近些年天津港詳細的地質、水文、地下水水位等資料，且由于天津港吹填造陸以及港區(qū)的快速發(fā)展和基礎建設，導致其中許多控制點監(jiān)測數(shù)據(jù)中斷，只得到其中幾年的監(jiān)測數(shù)據(jù)，造成數(shù)據(jù)的不完整問題。故選擇數(shù)理統(tǒng)計與隨機統(tǒng)計模型進行天津港未來幾年的沉降預測。張云霞[10]對隨機統(tǒng)計模型中的時間序列模型、灰色模型以及多元回歸模型進行分析比較，認為時間序列模型只考慮了地下水開采對沉降的影響，而天津港地面沉降還由吹填造陸引起，所以時間序列模型不能很好地代表天津港地面沉降，同時張云霞[10]在其天津地面沉降預測實驗中，發(fā)現(xiàn)灰色模型和多元回歸模型都能較好地預測天津地面沉降，且兩者精度相差不多。鑒于搜集的天津港地面沉降數(shù)據(jù)不夠全面，只有少數(shù)控制點具有從2002年到2018年的完整高程數(shù)據(jù)，且天津港地面沉降成因復雜，所以決定選用回歸模型中擬合效果更佳的多項式擬合來進行天津港地面沉降預測。而該多項式回歸模型也被許多學者用于地面沉降預測[11-13]。

2.2 多項式回歸模型算法介紹

回歸分析模型是一種十分重要的統(tǒng)計分析方法，在數(shù)據(jù)內(nèi)部具有某些規(guī)律的情況下，建立起一個或者多個自變量與因變量之間的關系，從而預測數(shù)據(jù)隨著時間的變化。一般根據(jù)自變量的數(shù)目將回歸分析模型分為一元回歸模型和多元回歸模型，而根據(jù)自變量的線性關系將回歸模型分為線性回歸模型和非線性回歸模型。采用的多項式回歸模型則屬于非線性回歸模型。

在研究或生活中常常會遇到這樣一類問題：假設有一堆離散點(xi,yi)(i=0,1,…,m)，需要找到一個函數(shù)，使得離散點到函數(shù)上的距離和最小。針對這類問題，數(shù)據(jù)擬合應運而生，而在所有的數(shù)據(jù)擬合算法中，多項式擬合無疑是被應用最廣泛的擬合算法。多項式擬合就是解求一個形如式(1)的函數(shù)，使得離散點到曲線上的距離和最小，如式(2)所示。

(1)

(2)

根據(jù)多元函數(shù)解算極值的原理，關于a0,a1,…,an的線性方程組，可以用矩陣表示為：

(3)

式(3)也叫作正規(guī)方程組或法方程組。

早有證明，方程組(3)的系數(shù)矩陣屬于對稱正定矩陣，故存在唯一解。從式中求得ak(k=0,1,…，n)，進而可以得到一個多項表達式：

這樣就可以得到一個離散點集的多項式擬合曲線。

多項式擬合應用廣泛，在信號處理，遙感數(shù)據(jù)處理、數(shù)學建模等各個方面都有應用，而在地面沉降預測方面，也有學者對其進行了研究。仲先強在對常澄高速公路的地面沉降預測研究中使用了多項式模型[13]，同時利用實測數(shù)據(jù)進行驗證，證明了多項式擬合方法對地面沉降檢測有著相當好的效果。除此之外，文獻[11][12][14]同樣在具體工程上證明了多項式擬合對地面沉降的預測有很好的效果，能夠較好地擬合地面沉降數(shù)據(jù)，并做出預測。

雖然多項式擬預測模型在地面沉降預測方面有較好的效果，但是也存在一些不足：多項式的階數(shù)較低時不能很好地擬合數(shù)據(jù)，但是當階數(shù)過高的時候會產(chǎn)生震蕩，預測精度會變差；多項式預測模型的各項的物理意義并不明確，不能用來表示造成地面沉降的各因素的影響程度。

3 沉降預測及分析

文章數(shù)據(jù)來源與王海剛[15]對天津市的研究一樣，只有水準點實測高程數(shù)據(jù)，所以同王海剛一樣選擇數(shù)理統(tǒng)計與隨機統(tǒng)計模型進行天津港地面沉降預測，但是由于水準點數(shù)據(jù)的殘缺，不能很好地用于多元線性回歸模型和灰色模型，所以選擇適用度更好的多項式模型來進行預測。同時為了保證預測精度，選用有連續(xù)超過10年以上高程數(shù)據(jù)的水準控制點(共5個)進行沉降預測，如表1所示，表中“—”表示沒有數(shù)據(jù)，圖1是這些水準控制點在天津港的位置。

表1 用于預測的水準點高程數(shù)據(jù) m

在地面沉降預測研究中，將表1中5個水準控制點的數(shù)據(jù)建立多項式模型，并以此預測2019年、2020年和2021年的高程值和沉降量。由于數(shù)據(jù)量有限，為了使?jié)M足多項式模型的解算條件，選擇二次和三次多項式模型作為其預測模型進行分析比較，選取精度最高的多項式模型進行預測，同時由于地面沉降的局部性特征，分別對5個監(jiān)測點建立多項式模型。

3.1 多項式模型次數(shù)選擇

為了模型預測結果顯示清晰，選取表1中G415，G408和G029的數(shù)據(jù)分別進行二次和三次多項式預測，結果如圖2，圖3所示。

對比圖2，圖3，可以發(fā)現(xiàn)：二次多項式模型的預測趨勢更符合實際情況，應選擇二次多項式模型作為天津港地面沉降預測的模型。

3.2 沉降預測及其精度分析

在選擇二次多項式模型作為天津港的地面沉降模型后，將2002年—2018年的所有數(shù)據(jù)作為輸入數(shù)據(jù)進行天津港的地面沉降預測，表2是其地面沉降預測結果。

從表2可以看出：G415和GM19更靠近塘沽老城區(qū)，預測沉降量最小，較其他點位更穩(wěn)定，G408和G029分別位于天津港北疆港區(qū)和南疆港區(qū)老港區(qū)，預測沉降量基本一致，G003位于新吹填造陸的天津港東疆港區(qū)，預測的沉降量最大，該區(qū)域最不穩(wěn)定，符合實際情況；以G415等5個點為代表的沉降趨勢反映，在天津港現(xiàn)有地下水位、工程建設等條件不變的情況下，天津港及周邊的地面依然保持下沉的趨勢，且較為嚴峻，每年的沉降量大概在10 mm～50 mm，區(qū)域整體年平均沉降速率在29 mm/年左右(取5個點2020年和2021年沉降量的平均值)。

表2 控制點高程預測結果

4 結論與建議

選擇多項式回歸模型對天津港地面沉降進行了預測，利用歷史沉降數(shù)據(jù)對多項式模型的階數(shù)進行了研究，發(fā)現(xiàn)對于現(xiàn)有的天津港歷史沉降數(shù)據(jù)，相對于三次及以上的多項式預測模型，二次多項式回歸模型更能精確模擬天津港的地面沉降規(guī)律；利用二次多項式模型對天津港監(jiān)測周期超過10年的5個水準點進行地面沉降預測，最終結果顯示：以G415等5個點為代表的沉降趨勢反映，在現(xiàn)有的地理環(huán)境、人為因素等不變的情況下，天津港的地面沉降仍將繼續(xù)，且較為嚴峻，地面會以29 mm/年左右的速度持續(xù)下沉。

基于對天津港地面沉降的預測分析，有以下幾點建議：

1)建立天津港地面沉降監(jiān)測網(wǎng)，加強地面沉降監(jiān)測?，F(xiàn)有天津港的監(jiān)測數(shù)據(jù)量非常有限，主要來源于天津港的基礎控制網(wǎng)測量數(shù)據(jù)。

2)禁止開采地下水資源。天津港靠近內(nèi)陸地區(qū)，包括濱海新區(qū)等區(qū)域，地面沉降的主要原因是地下水資源的過度開采，需要控制地下水的開采，開發(fā)其他地面水資源作為企業(yè)和居民用水，同時還應該對本地區(qū)進行地下水回灌，讓該地區(qū)的地下水資源不再下降。

3)提高地基處理技術水平和驗收標準。天津港地區(qū)由于吹填造陸形成的陸域的沉降量較大，針對該區(qū)域，應該提高地基處理技術水平和驗收標準，使該地區(qū)地基更加穩(wěn)固。