福建省地震地形變學(xué)科發(fā)展研究報告

福建省地震學(xué)會

?

福建省地震地形變學(xué)科發(fā)展研究報告

福建省地震學(xué)會

該文詳細介紹了近40年來，福建省地震地形變學(xué)科發(fā)展的概況。福建省地形變觀測網(wǎng)是由全省GPS連續(xù)觀測網(wǎng)、重力、水準(zhǔn)、地磁流動觀測網(wǎng)以及重力、傾斜與應(yīng)變觀測臺站組成。經(jīng)過“十五”、“十一五”、“十二五”的發(fā)展，福建省地殼形變臺網(wǎng)建設(shè)得到了進一步加強，從而推動了我省地殼形變的地震預(yù)測與研究工作的快速發(fā)展，在推進全國地殼形變學(xué)科發(fā)展中也起到了積極的作用。

地震 地形變 發(fā)展研究 福建省

0 引言

我國地殼形變的地震預(yù)測與研究工作已開展40多年，并形成了具有中國特色的相當(dāng)規(guī)模的地殼形變預(yù)報地震觀測和研究體系，在推進地殼形變學(xué)科發(fā)展中起到了積極的作用。

地殼形變是地殼介質(zhì)在內(nèi)生的構(gòu)造應(yīng)力和外生的天體引力以及地表荷載力作用下發(fā)生的一種表現(xiàn)形式。與地震有關(guān)的地殼形變則是在孕震過程中，隨著構(gòu)造應(yīng)力的不斷積累，直到巖石發(fā)生破裂前后的一種最直接的伴隨現(xiàn)象，這種構(gòu)造應(yīng)力積累直到巖石快速破裂，一般需要數(shù)十年甚至更長的時間，因此與其伴生的地殼形變呈現(xiàn)出長、中、短、臨的時空變化圖像。

對上述具有不同周期的地殼形變，要采用具有不同頻率響應(yīng)的觀測儀器進行觀測。對于由地震激發(fā)的高頻地殼運動，可采用地震儀進行觀測，這種儀器的最佳頻率范圍通常從十分之一秒到幾分鐘，而監(jiān)測非常緩慢的地殼形變則可采用常規(guī)的大地測量手段，至于監(jiān)測介于這兩者之間的頻率范圍的地殼運動，可通過在特別感興趣的地區(qū)內(nèi)增加大地測量的觀測次數(shù)，布設(shè)全球性的長周期和超長周期地震儀臺網(wǎng)以及發(fā)展地球物理和大地測量新技術(shù)來達到。在發(fā)展地球物理和大地測量儀器方面，近20年來，我國地殼形變觀測技術(shù)有了長足的發(fā)展，特別是隨著數(shù)字化觀測儀器（重力、傾斜和應(yīng)變儀）和空間大地測量技術(shù)（如GPS）的發(fā)展，使我們可以對高至地震頻率、低至零頻這樣一個非常寬闊的頻率范圍內(nèi)的地殼運動進行監(jiān)測，使我國的地殼形變觀測研究進入了一個新的發(fā)展階段。這些新進展主要表現(xiàn)在兩方面：一是在空間尺度上把地殼形變觀測的各個手段聯(lián)成一個整體的觀測系統(tǒng)，從而有可能獲得真實的地殼運動圖像；二是幾何觀測與物理（重力）觀測相結(jié)合，使地殼形變研究工作由幾何學(xué)向運動學(xué)、動力學(xué)方向發(fā)展，使人們有可能定量化地研究地殼形變動力學(xué)過程，從而可以對地殼形變進行一定的預(yù)測，推動地震預(yù)報及地震災(zāi)害的預(yù)測研究工作。觀測系統(tǒng)的長足進展日臻完善，為采用地形變預(yù)測預(yù)報地震提供了堅實而可靠的基礎(chǔ)。

地殼形變觀測對地震理論認識的最大貢獻是在1906年美國舊金山地震前后進行的三角測量，從而導(dǎo)致了美國大地測量學(xué)家里德系統(tǒng)地闡述彈性回跳理論?，F(xiàn)在關(guān)于震源的很多知識是來自震前和震后大地測量資料的比較。比較的價值取決于震前與震后測量的精度。同時，大地測量精度的改進使對較小地震的震源機制研究成為可能，并有可能揭示地震的前兆。

實驗和理論研究表明，在應(yīng)變積累過程中由于地殼介質(zhì)的不均勻性，會在破裂（地震）之前出現(xiàn)明顯的在時間上分階段的不同周期的地殼形變，這就是目前采用大面積形變測量、跨斷層形變測量和固體潮汐形變測量方法進行地震預(yù)報的理論基礎(chǔ)。

1 地殼形變預(yù)報地震的科學(xué)思路

地震大都發(fā)生在活動板塊的邊界和活動斷裂帶上，在這些條帶上發(fā)生突發(fā)性位移（地震）之前，地面位移、應(yīng)力場會相繼發(fā)生長、中、短、臨變化，而這些變化又是與孕震過程的階段性有關(guān)，且前兆反映一般以群體形式顯示。總之，根據(jù)地形變在時空上的這些變化特征，并結(jié)合當(dāng)?shù)氐牡刭|(zhì)構(gòu)造背景，可以對大震前地殼應(yīng)力—應(yīng)變狀態(tài)進行跟蹤，以便圈定地殼形變較顯著的地區(qū)作為今后要發(fā)生地震的可能地點。因此，研究判定與地震相關(guān)的地形變前兆，可采用以下物理判據(jù)與方法：

1.1 大面積地形變與地質(zhì)構(gòu)造的一致性

短臨前兆出現(xiàn)之前的大面積地形變分布主要反映了地殼大范圍運動的概貌。正常的大面積地殼運動具有速度緩慢，變化形態(tài)與地殼淺部構(gòu)造分布一致，與大區(qū)域應(yīng)力場在方向和強度上相適應(yīng)的特點。因此構(gòu)造形變都是沿活動斷裂帶分布的，當(dāng)局部地區(qū)的構(gòu)造單元運動速率明顯改變，在其邊緣出現(xiàn)垂直形變高梯度帶和重力梯度帶，當(dāng)它們的分布與地質(zhì)構(gòu)造一致時，則反映了這些構(gòu)造單元相對運動，從而反映現(xiàn)今構(gòu)造應(yīng)力場的作用方式。構(gòu)造單元的邊緣往往容易出現(xiàn)應(yīng)變、垂直形變和重力等值線形態(tài)的轉(zhuǎn)折和密集，形成高梯度帶。因此，短臨前兆出現(xiàn)之前的大面積地殼形變及其形成的梯度帶是判定地形變短臨前兆的物理判據(jù)之一。

1.2 地形變異常的階段性

地震是地殼以突然破裂的形式表現(xiàn)出來的一種地殼變形。在主破裂（地震）前后，地殼演變過程具有長期緩慢變化、快速變化、震前突變和震后調(diào)整四個階段特征，這是與地殼的非均勻性分不開的。因此，在短臨前兆出現(xiàn)之前是否會出現(xiàn)長期和中期地殼形變是判定短臨前兆的物理判據(jù)之一。因為根據(jù)地殼介質(zhì)的非均勻性，在震前，地殼介質(zhì)會經(jīng)歷非線性—線性—非線性變形三個階段。就破裂過程而言，會經(jīng)歷張性微破裂—剪切破裂—主破裂三個基本過程。這兩者從不同角度反映了受破裂過程控制的地殼形變大致分為長期、中期和短臨三個階段。

1.3 不同地區(qū)地形變異常的差異性和復(fù)雜性

在大震發(fā)生前，不同地區(qū)的地形變異常有很大的差異，也表現(xiàn)出復(fù)雜性。這與不同地區(qū)大地構(gòu)造單元的分布及其特點，地殼介質(zhì)的性質(zhì)，主斷裂帶的分布、產(chǎn)狀與活動水平有關(guān)。一般來說，對正斷層而言，垂直形變較顯著；對走滑斷層而言，水平位移較顯著；在兩條或多條斷層交匯或轉(zhuǎn)折處，垂直形變和水平位移都較顯著。大震之前地殼形變的這種差異性和復(fù)雜性正是不同地區(qū)存在不相同構(gòu)造背景的一種表現(xiàn)形式。

1.4 地形變短臨前兆特征

震前地形變短臨前兆種類上呈現(xiàn)出群體性，空間上表現(xiàn)出集中性，時間上出現(xiàn)有序性。地形變前兆群體性是指在震前會陸續(xù)出現(xiàn)種種類型的地形變，不過各次地震前有所不同，因為從整體大地測量的角度看來，任何一種形變都可表示為地點、時間和重力位的函數(shù)，這種地形變前兆的群體性體現(xiàn)了各種地形變量、重力與應(yīng)力之間的物理上的相關(guān)性。在震前，不但地形變前兆表現(xiàn)出群體性，而且其它手段（如地下水位、水化學(xué)萬分、地電、地磁）亦會出現(xiàn)異常，這正是各種物理量和水化學(xué)組分對震源過程的不同響應(yīng)。地形變前兆在時間上的有序性是指在應(yīng)變積累的過程中，由于地殼介質(zhì)的不均勻性，會在破裂（地震）之前出現(xiàn)明顯的不同階段，會出現(xiàn)周期不同的地形變，越接近破裂，響應(yīng)就越明顯。地形變前兆在空間上的集中性是指孕震過程發(fā)展到某一特定階段后，各種地形變異常會逐漸轉(zhuǎn)移到未來震中附近的一個較小范圍內(nèi)。

2 地殼形變監(jiān)測地震的基本方法

地殼形變可理解為可變形地塊的起始狀態(tài)（未變形）和終止?fàn)顟B(tài)（已變形）之間的地塊位形的變化，這種位形變化包括地塊質(zhì)點的位置變化，質(zhì)點之間的距離及方向變化，后者即為塊體的變形。地殼形變是由于地殼質(zhì)點的位移所引起的，位移為質(zhì)點的位置向量在兩個不同時刻之差，單位時間的位移即為該質(zhì)點的地殼運動速率。根據(jù)彈性力學(xué)理論，任何變形都是由平移、旋轉(zhuǎn)和應(yīng)變?nèi)糠纸M成。地殼形變從空間分布來分，可分為全球性的、區(qū)域性的和局部性的；從時間尺度來分，可分為長期性的、周期性和間歇性的；從對作用力的響應(yīng)來分，可分為彈性形變和粘滯性形變；從地形變的方向來分，可分為垂直形變和水平形變，對這兩個分量進行矢量合成，最終可確定地殼形變的方向和大??；從形變的成因來分，又可分為構(gòu)造形變、地球潮汐形變和地殼載荷形變和突發(fā)性形變，前者與板塊運動和斷層活動有關(guān)，中者與天體起潮力有關(guān)，次者與海潮和水庫蓄水有關(guān)，后者與火山和地震有關(guān)。

不管是垂直形變觀測還是水平形變觀測，地面上的一切觀測都是在地球重力場中進行的。因此在進行高精度的地形變觀測時，配合進行高精度的重力測量，不但可以對測量結(jié)果進行固體潮改正，而且可以把幾何與物理測量結(jié)合起來對地殼形變的動力學(xué)過程進行研究。

震前地形變類型決定了監(jiān)測震前地形變應(yīng)采用的觀測手段。要監(jiān)測地面的垂直位移，需采用水準(zhǔn)儀；對水平位移，可采用基線、激光測距和GPS等手段；對于旋轉(zhuǎn)，可采用經(jīng)緯儀進行方向觀測以確定方向改變量，或者通過應(yīng)變觀測經(jīng)解算確定；要確定一個區(qū)域內(nèi)的應(yīng)變場狀態(tài)，可根據(jù)該區(qū)域內(nèi)諸點的位移改變量通過解算而獲得，對于某一點的應(yīng)變狀態(tài)，可采用應(yīng)變儀對三個以上方向進行觀測得到；對于震前因地殼介質(zhì)密度變化以及點位高程變化所引起的重力變化，可采用重力儀進行觀測。有了好的觀測手段，又有考慮縝密周到的多層次的觀測臺網(wǎng)，才能根據(jù)所監(jiān)測到的地形變資料作長、中、短、臨地震預(yù)測預(yù)報。目前我國是通過以下方法與途徑來進行預(yù)測的：

①監(jiān)測大面積重力場、位移場和應(yīng)變場的中長期趨勢變化。主要觀測手段為精密水準(zhǔn)測量、水平形變測量、流動高精度重力測量和空間大地測量?，F(xiàn)有的水準(zhǔn)網(wǎng)大部分是在原國家一、二等水準(zhǔn)網(wǎng)的基礎(chǔ)上改造、擴建而成，按一等水準(zhǔn)要求進行復(fù)測，復(fù)測和聯(lián)測精度為10-6，達到國際同類觀測水平。水平形變監(jiān)測網(wǎng)是由GPS站所構(gòu)成，分為連續(xù)觀測和定期復(fù)測兩種，GPS觀測的精度可達到毫米級。

②監(jiān)測斷層中、短期運動學(xué)和動力學(xué)特征。斷層形變觀測方法有短水準(zhǔn)測量和連續(xù)自記動態(tài)測量兩種方法，精度可達10-6~10-7。

③監(jiān)測固體潮短臨前兆變化特征。固體潮是指固體地球?qū)μ祗w起潮力的3種不同響應(yīng)方式，分別為重力、傾斜和應(yīng)變固體潮。天體起潮力對地震發(fā)生可起到一種觸發(fā)作用，在地震發(fā)生之前，這3種固體潮會發(fā)生畸變，潮汐因子會隨之變化。為了監(jiān)測這些變化，我國已建立了固體潮監(jiān)測臺網(wǎng)，該臺網(wǎng)是由重力、傾斜和應(yīng)變固體潮觀測臺站組成，重力儀的精度可達1微伽（10-8m/秒2），傾斜儀的精度可達到1ms (0.001角秒)，應(yīng)變儀精度可達10-8~10-10。

在我國，從20世紀(jì)90年代起，已初步建立了一個較完整的地殼形變監(jiān)測系統(tǒng)，該系統(tǒng)是由中國大陸地殼形變監(jiān)測、區(qū)域地殼形變監(jiān)測、斷層形變監(jiān)測和固體潮監(jiān)測網(wǎng)等4個子系統(tǒng)組成。中國大陸形變網(wǎng)主要采用VLBi、SLR、GPS、精密水準(zhǔn)測量、高精度重力測量以監(jiān)測中國大陸地殼運動全貌、次板塊相對運動，其作用是提供基礎(chǔ)背景資料，參與全球地球動力學(xué)研究。區(qū)域形變網(wǎng)主要采用GPS、精密水準(zhǔn)和流動重力測量，對該區(qū)域的地震危險區(qū)地殼形變圖像及其進程進行監(jiān)測，其作用是提供長、中期地震前兆信息。斷層形變監(jiān)測網(wǎng)主要采用短水準(zhǔn)、測邊小網(wǎng)、重力測量、連續(xù)自記動態(tài)觀測，以監(jiān)測主要斷層的活動狀態(tài)，提供中、短、臨地震信息。固體潮形變監(jiān)測網(wǎng)所采用的主要觀測方法為地傾斜、地應(yīng)變與重力，其監(jiān)測的內(nèi)容為固體潮波形畸變，潮汐因子的變化以及地殼介質(zhì)的物性變化，以期提供震前中、短、臨信息。

3 國內(nèi)外地震地形變學(xué)科發(fā)展現(xiàn)狀

3.1 國外地震地形變學(xué)科發(fā)展現(xiàn)狀

大震發(fā)生后，人們?？煽吹降缆峰e開，鐵軌扭曲，地面開裂、沉陷或隆起，冒沙噴水等一派觸目驚心的地殼變動景象。人們不禁要問，震后地面會發(fā)生如此有目共睹的地殼形變，地震前地殼亦會發(fā)生形變嗎？回答是肯定的。不過，一般說來，這種形變不是肉眼所能察覺出來的，它是一種緩慢的形變過程，要用精密的觀測手段才能檢測出來。也許美國人里德是世界上最早企圖用測量地殼運動的方法研究震前地殼形變理論的一位學(xué)者，他根據(jù)1851年~1855年間的三角測量結(jié)果和震后立即進行的測量結(jié)果分析了1906年美國舊金山8.3級大地震，根據(jù)這些分析結(jié)果提出了彈性回跳理論。這個理論認為，由于圣安德列斯斷層的相對運動導(dǎo)致了地殼中彈性應(yīng)變能的積累，當(dāng)應(yīng)變達到臨界點時，斷層回跳產(chǎn)生地震波。根據(jù)他的計算，震前50年間斷層相對位移大約3.2m，震時最大滑動量達6.5m，據(jù)此計算舊金山地震的復(fù)發(fā)周期為100年。上世紀(jì)80年代對圣安德列斯斷層的開挖調(diào)查結(jié)果，認為該區(qū)大地震的平均復(fù)發(fā)時間為150年或200年。里德所采用的長期地震預(yù)報方法基本上與現(xiàn)在一樣，亦是以預(yù)測地殼應(yīng)變程度來進行預(yù)報。

從物理學(xué)觀點看，斷層彈性回跳理論對發(fā)生地震的本質(zhì)給出了較為完善的解釋，從而被各國地震學(xué)家所接受。后來，由簡單的彈性回跳模型發(fā)展到斷層位錯理論，繼又轉(zhuǎn)向?qū)φ鹪次锢磉^程的研究，使我們對孕震過程中產(chǎn)生地形變的機制及其演化過程了解得更為清楚。既然地殼形變是由斷層的位錯所引起的，當(dāng)然就可以根據(jù)斷層的幾何參數(shù)（走向、傾向、傾角、長度、寬度、走向滑動錯距、傾向滑動錯距、斷層上下邊界距地殼的距離、斷層埋深）來計算斷層錯動在地表面的位移，反之，亦可根據(jù)大地測量資料反演斷層的震源參數(shù)。

在日本，今村明恒在第二次世界大戰(zhàn)前就曾試圖用測量的方法進行地震預(yù)報，他認為地震是由地塊傾斜引起的。應(yīng)他的要求，日本陸地測量局在第二次世界大戰(zhàn)前對很多水準(zhǔn)路線進行了復(fù)測，其中最有意義的是1944年東南海地震臨震前在靜風(fēng)岡縣御前崎附近進行的測量，測量結(jié)果發(fā)現(xiàn)震前有異常的地面傾斜。

1980年，美國地殼運動測量專業(yè)委員會在其所提出的題名為《地殼形變監(jiān)測戰(zhàn)略》的一份研究報告（下簡稱報告）中指出，從地震資料得到關(guān)于斷層活動細節(jié)的信息是非常困難的。在地震信號的記錄中，大部分能量局限于0.1～10秒的范圍內(nèi)，盡管少量長周期地震儀可測到1000秒的地震動，較長周期的斷層運動是不能用地震方法監(jiān)測的，但實際上它可改變應(yīng)力分布。非彈性效應(yīng)和余震發(fā)生率有很大關(guān)系。因此，這非常類似于緩慢的非彈性過程，此過程可改變應(yīng)力分布而不釋放地震能量，而大地測量方法的主要貢獻在于此種方法能監(jiān)測到地震頻帶寬度以外的準(zhǔn)靜態(tài)位移。如果經(jīng)典的大地測量技術(shù)能達到所要求的精度，則它們完全適用于此目的。美國地殼運動測量專業(yè)委員會的這份報告對全球范圍測量、活動板塊邊界附近測量、板內(nèi)形變和其它新構(gòu)造運動的測量，測量構(gòu)造形變的儀器設(shè)備以及有關(guān)監(jiān)測構(gòu)造形變的大地測量戰(zhàn)略作了詳細深刻的闡述。特別是對地面觀測技術(shù)中的應(yīng)變儀、傾斜儀、重力儀、激光測距裝置、氣泡水準(zhǔn)測量的原理、精度及其布設(shè)方案作了較詳細的介紹，在空間觀測技術(shù)方面，對人造衛(wèi)星激光測距（SLR）、甚長基線無線電干涉儀（VLBi）、全球定位系統(tǒng)（GPS）以及測量位置的其他空間技術(shù)的原理、精度、運用也作了全面的介紹。在這里，要特別指出，隨著GPS精度及采樣率的提高，今天GPS成了用來監(jiān)測地殼地形變的一種多時空尺度、高分辨率的實時動態(tài)觀測系統(tǒng)，且可全天候、不要求通視、費用低廉的高技術(shù)手段。上述的空間測量技術(shù)與地面觀測技術(shù)已成為上世紀(jì)80年代以來美國監(jiān)測地震帶內(nèi)地殼各種運動的主要手段，其主要的任務(wù)分別為在地震帶內(nèi)和地震帶周圍，監(jiān)測其長期應(yīng)變積累圖像及速率，監(jiān)測可能的中期應(yīng)變的變化，監(jiān)測大震前可能出現(xiàn)的短期前兆。

美國從上世紀(jì)80年代開始，除了研制地面和空間測量技術(shù)以外，還對如何制定一個切實可行的觀測方案進行了規(guī)劃，例如在研究美國西部地殼運動時，建議對所采取的三邊和其它大地測量方案作如下實質(zhì)性擴充：增大觀測頻度、大量加密觀測點、臺網(wǎng)的布設(shè)要擴大到離主要斷層更遠的地方、盡量把主要斷裂帶都包括進去；對一個大的地震帶，該報告認為地殼形變測量的主要目標(biāo)應(yīng)進一步提高精度，擴大監(jiān)測范圍，它應(yīng)包括盡可能以最高的精度對擴充的地形變網(wǎng)進行半年或一年一次的復(fù)測工作，對沿主要斷層布設(shè)點距為10公里或20公里的兩點之間的距離進行定期的復(fù)測工作（約2周復(fù)測一次），在小范圍內(nèi)，采用定點監(jiān)視儀器進行幾乎連續(xù)的觀測，以便檢測大震之前的任何短期前兆；要用空間技術(shù)把形變監(jiān)測網(wǎng)擴大到美國之外的西半球和西太平洋的一些主要地震帶上去，并和這些國家進行協(xié)作；應(yīng)該發(fā)展海底測量所必須的觀測儀器；應(yīng)該改進覘標(biāo)安裝和評價其穩(wěn)定性工作，這些建議被許多國家所借鑒。

在如何利用地形變觀測資料進行研究方面，美國學(xué)者如普雷斯提出了利用位錯理論根據(jù)震前震后大地測量資料反演斷層震源參數(shù)，保加利亞學(xué)者維尼迪柯夫提出了利用固體潮觀測資料計算潮汐因子與相位滯后的數(shù)字濾波方法，潮汐因子反映震前地殼物性的變化可作為短臨前兆，如何剔除固體潮觀測資料中有問題的觀測數(shù)據(jù)，日本學(xué)者Nakaii提出了Nakaii方法，日后被普遍采用。

在地震成因研究方面，從地殼形變的角度來看，除早期根據(jù)大地形變測量資料提出彈性回跳模型外，美國學(xué)者于上世紀(jì)80年代根據(jù)巖石實驗和實測資料進一步發(fā)展提出了包括蠕變—粘滑現(xiàn)象的斷層模型，繼后美國學(xué)者肖爾茨又提出了膨脹擴容模式，無論哪一種模式均表明在地震的孕育、發(fā)生和發(fā)展過程中，由于地殼構(gòu)造的復(fù)雜性和介質(zhì)的不均勻性，會在某些特殊部位，例如活動斷裂帶出現(xiàn)劇烈的變化。

前蘇聯(lián)的梅謝里柯夫、美國的肖爾茨、日本的藤井陽一郎等學(xué)者根據(jù)大量震例總結(jié)出形變前兆的階段性特征，即穩(wěn)定變化段α，加速變化段β，失穩(wěn)段β2及震前突變段γ。圖1為巖石應(yīng)力—應(yīng)變本構(gòu)關(guān)系示意圖，圖2是藤井陽一郎根據(jù)30個震例歸納的應(yīng)變積累階段特征圖，圖2中的α段對應(yīng)于圖1中的AB段，AB段表示地殼出現(xiàn)線彈性形變；β階段對應(yīng)于BC段，表示地殼的非彈性變形明顯，非彈性體積應(yīng)變使巖石出現(xiàn)膨脹現(xiàn)象，隨著應(yīng)力增加，應(yīng)變速率增大亦加快，C點表示巖石在一定條件下所能隨的最大載荷；γ階段可對應(yīng)于CD段，此段與地震發(fā)生有著非常密切的關(guān)系，越過C點以后，地殼承載能力下降，應(yīng)變加速，仿佛巖石隨應(yīng)力增加而變軟；δ階段對應(yīng)于DE段，表示地殼破裂已完成了，斷裂面已經(jīng)形成。

圖1 巖石應(yīng)力-應(yīng)變本構(gòu)關(guān)系示意圖

圖2 震前地形變前兆的階段性特征

隨著地形變觀測資料的增多，發(fā)現(xiàn)在同一個地點的震前地形變圖像呈現(xiàn)復(fù)雜性，沒有確定的前兆模式可言，為了解釋前兆的復(fù)雜性，西方學(xué)者分別從耗散結(jié)構(gòu)理論、分形理論以及非線性理論角度對這種復(fù)雜性進行了研究，發(fā)現(xiàn)當(dāng)外部傳遞給變形體系統(tǒng)的能量轉(zhuǎn)變?yōu)閹r石內(nèi)部應(yīng)變能，當(dāng)應(yīng)變能積累到一定程度時，巖石開始軟化，將巖石內(nèi)的應(yīng)變能轉(zhuǎn)變?yōu)闊崮芎驼饡r的彈性能，即將利用率較高的能量變?yōu)槔寐瘦^低的能量，這個過程就是巖石能量的耗散過程，在此過程中，地殼形變圖像呈現(xiàn)無序和有序同在，自組織和混沌并存局面，使對地形變分析推進到一個新局面。

鑒于地面觀測系統(tǒng)（例如傾斜儀、應(yīng)變儀）易受到外界干擾，許多學(xué)者轉(zhuǎn)而致力于深井地球特理綜合觀測，尤其是在科技發(fā)達的國家得到了迅速的發(fā)展，建立了大量深井觀測站或長期觀測實驗點。在“地震之國”日本，深井觀測站的深度達1000～3800m，其中取得成功經(jīng)驗的深井觀測研究項目有伊豆半島地震群、神岡礦山活斷層以及南非金礦的觀測與研究等。在2003年，東京大學(xué)、東濃地震研究所、名古屋大學(xué)等在日本本洲中部岐阜縣屏風(fēng)山活斷層附近1020m的深井中，設(shè)置了世界上最先進的長期綜合地球物理深井觀測儀器，并進行了初期應(yīng)力測量，結(jié)果表明，與地面觀測相比較，深井觀測能提高觀測精度1~2個數(shù)量級，如地應(yīng)變能觀測到10-10的極微弱應(yīng)變變化。美國從1992年底在著名的圣安德列斯斷層開展了深部長期觀測研究項目，又稱“地球透鏡計劃”，其中一個重要的組成部分為“板塊邊界觀測（PBO）計劃”，使用了大量的鉆孔應(yīng)變儀，目前約有80個Gladwin 三分量鉆孔應(yīng)變觀測點。德國亦已實施德國大陸超深鉆井計劃，簡稱KTP項目，該項目是德國地球科學(xué)領(lǐng)域規(guī)模最大，耗資最多的一個科學(xué)項目。

3.2 國內(nèi)地震地形變學(xué)科發(fā)展現(xiàn)狀

3.2.1臺網(wǎng)建設(shè)與觀測技術(shù)

我國地殼形變監(jiān)測系統(tǒng)按現(xiàn)有的觀測（臺）網(wǎng)綜合分類，它是由固體潮汐形變監(jiān)測、斷層形變監(jiān)測、區(qū)域地殼形變監(jiān)測和中國大陸地殼形變監(jiān)測四部分組成。區(qū)域地殼形變測量和中國大陸地殼形變測量，習(xí)慣上稱為大面積形變測量，主要觀測手段為精密水準(zhǔn)測量、水平形變測量、流動高精度重力測量和空間大地測量。精密水準(zhǔn)測量采用德國蔡司廠的Ni002精密水準(zhǔn)儀施測，高精度重力測量采用拉柯斯特—隆貝格重力儀、德國的GS重力儀和國產(chǎn)的DZW型重力儀，觀測精度達1微伽，空間大地測量采用GPS。上世紀(jì)90年代起普遍推行了數(shù)字化觀測，水準(zhǔn)測量采用數(shù)字水準(zhǔn)儀，數(shù)字水準(zhǔn)儀有德國蔡司廠的DiNi系列、瑞士徠卡廠的NA3003和日本拓普康生產(chǎn)的DL101等。我國在“九五”期間，選定了DiNi系列水準(zhǔn)儀作為地震系統(tǒng)斷層形變定點臺站水準(zhǔn)測量的監(jiān)測儀器，DiNi系列中的DiNi11數(shù)字水準(zhǔn)儀每公里往返測觀測中誤差為±0.30mm，測高最小讀數(shù)單位為0.01mm，測距最小讀數(shù)單位為1mm，單個測量時間為4s，而空間大地測量采用的GPS接收機則為美國AHTECH CGRS GPS接收機，該款GPS接收機為美國AHTECH公司出品的連續(xù)大地參考站觀測系統(tǒng)，是高精度科學(xué)研究、陸地測量、地理信息系統(tǒng)和工程應(yīng)用永久性GPS基準(zhǔn)站廣泛采用的設(shè)備之一，該接收機系統(tǒng)誤差（GPS衛(wèi)星間），偽距（全波段）<1.00cm，載波相位（全波段）<0.01cm。在國家計委、科技部支撐的國家重大科學(xué)工程“中國地殼運動觀測網(wǎng)絡(luò)”項目中，除上述的AshtechZ-X11外，還引進了如Trimble4000SSI、Rouge8100等GPS接收機。

中國地震局牽頭的“中國地殼運動工程網(wǎng)絡(luò)”由VIBI.SLR、連續(xù)與超導(dǎo)重力站與GPS二期基準(zhǔn)站、一期區(qū)域站、二期區(qū)域站與二期基本站組成，見圖3。

圖3 中國地殼運動工程網(wǎng)絡(luò)臺站分布圖

固體潮汐形變監(jiān)測系統(tǒng)是由地傾斜、地應(yīng)變和重力固體潮等3個子系統(tǒng)構(gòu)成。這樣一個系統(tǒng)經(jīng)歷了3個階段的發(fā)展，在第一個階段，以裝備低靈敏度、目視觀測儀器（第一代）為主，第二個階段以裝備高靈敏度（10-8~10-9）、能清晰自動連續(xù)記錄潮汐形變的儀器（第二代）為主，第三階段則是以模擬記錄為特征的第二代儀器向數(shù)字化儀器轉(zhuǎn)換為主，在實現(xiàn)數(shù)字化觀測后，固體潮汐監(jiān)測系統(tǒng)就可通過不同時間的采樣，在一個非常寬闊的頻帶范圍內(nèi)對地殼動力學(xué)現(xiàn)象進行觀測，以監(jiān)測地殼運動的方向、速度及其巖石特性的變化，從而填補了測震學(xué)和大地測量學(xué)之間頻率范圍的空白，有利于監(jiān)測短期前兆向臨震前兆的過渡和發(fā)展，并有利于監(jiān)測震后的地殼調(diào)整過程（蠕變或松弛過程）。

地傾斜固體潮監(jiān)測子系統(tǒng)目前所采用的數(shù)字觀測技術(shù)有DSQ型水管傾斜儀和VS垂直擺傾斜儀，前者的靈敏度為1ms（0.001"），日漂移量0.005",標(biāo)定精度優(yōu)于1%，全自動標(biāo)定；后者的靈敏度為0.1ms，日漂移小于0.005"。

2010年全國傾斜觀測臺網(wǎng)正式運行并納入學(xué)科管理的臺站共161個。運行儀器包括水管傾斜儀、垂直擺傾斜儀、水平擺傾斜儀、鉆孔井下傾斜儀等類型，共計有252套儀器，主測項513個分量，觀測方式有數(shù)字化觀測與模擬觀測兩種。2010年度傾斜儀分布圖見圖4。

地應(yīng)變固體潮監(jiān)測子系統(tǒng)目前所用的數(shù)字化儀器有SS-Y型伸縮儀、TJ型體積式鉆孔應(yīng)變儀、YRY-4型4分量鉆孔應(yīng)變儀和RZW型4分量鉆孔應(yīng)變儀。SS-Y型伸縮儀的觀測精度優(yōu)于10-9，日漂移<1×10-8；TJ型體積應(yīng)變儀的精度≤1×10-9，且自身穩(wěn)定性≈1×10-7/a；YRY-4型鉆孔應(yīng)變儀精度為10-9；RZB型精度為10-8。

2010年度應(yīng)變觀測臺網(wǎng)正式運行并納入學(xué)科管理的洞體應(yīng)變臺站96個、觀測儀器96套。在運行儀器包括伸縮儀、體應(yīng)變儀、分量應(yīng)變儀等類型，共計有227套儀器，主測項496個分量，觀測方式分?jǐn)?shù)字化觀測與模擬觀測兩種。2009年度應(yīng)變儀分布圖見圖5。

圖4 全國地傾斜臺站分布圖

圖5 全國洞體應(yīng)變臺站分布圖

重力固體潮監(jiān)測子系統(tǒng)所采用的數(shù)字化觀測儀器有DZW型重力儀和GS型重力儀。DZW型重力儀的讀數(shù)靈敏度為0.1Mgal，GS數(shù)字化重力儀是在對從德國ASKANIA公司引進的GS重力儀改造而成，改造后的精度為1×10-8ms-2。我國現(xiàn)有重力臺網(wǎng)包括39臺（套）重力儀，分布在14個省份，3個直轄市、5個自治區(qū)內(nèi)。重力臺站分布見圖6。

圖6 全國重力臺站分布圖

上述SS-Y型伸縮儀、DSQ型水管傾斜儀、VS垂直擺傾斜儀以及DZW型重力儀分別為中國地震局地震研究所呂寵吾、聶磊、胡國慶研究員研制，對GS重力儀改造也是該所進行的，YRY-4型鉆孔應(yīng)變儀由河南省鶴壁市地震局池順良教授研制，DZW型鉆孔應(yīng)變儀由中國地震局應(yīng)力研究所歐陽祖熙教授研制，TJ型體應(yīng)變儀由應(yīng)力所蘇愷之研究員研制。上述地形變儀器的性能與精度均已達到國際同類儀器的先進水平。有些儀器例如中國地震局地震研究所蔡惟鑫研究員等研制的石英伸縮儀曾于上世紀(jì)80年代~90年代遠赴歐洲比利時、西班牙、盧森堡進行地球動力學(xué)觀測，獲得好評，并獲得盧森堡大公國勛章。

3.2.2天體起潮力位展開及固體潮理論研究

起潮力為天體（月亮）對地球表面某一點的萬有引力與月亮對地心的萬有引力的矢量差，在地心起潮力為零。月亮在該點的起潮力位等于月亮在該點的引力位與該點繞月、地公共質(zhì)心作圓周運動的向心力位之差。1922年，杜森（Doodson）將起潮力位按照幅角數(shù)對潮汐波進行了自然分類，把起潮力位表示成振幅為常數(shù)的386個分諧波之和，每個諧波的振幅都大于0.00001D，其中D稱為杜森常數(shù)，1967年前D=26206cm2/s2，1967年后D=26277 m2/s2。起潮力位是推導(dǎo)重力、傾斜與應(yīng)變固體潮理論值的基礎(chǔ)，起潮力位展開的諧波個數(shù)越多，則固體潮理論值的精度越高，對于高精度固體潮觀測，對其觀測資料作調(diào)和分析，則要求具有更高精度的起潮力位展開式，即要求組成該展開式的諧波數(shù)要遠遠大于杜森展開式的386個諧波。在1922年那個年代，沒有電子計算機，杜森僅靠手算完成了如此復(fù)雜費時的諧波展開，確非易事。到了上世紀(jì)80~90年代，重力固體潮觀測的精度已達到微伽級，傾斜觀測的墨水記錄已達到1ms/mm，應(yīng)變觀測已達到10-9~10-10，此時就要求有諧波展開個數(shù)大于杜森展開個數(shù)更多的起潮力位展開式。

中國地震局分析預(yù)報中心郗欽文研究員從上世紀(jì)80年代開始，利用早已發(fā)展的計算機演繹方法，在國家自然科學(xué)基金（批準(zhǔn)號：49574217）的資助下，自行發(fā)展計算機公式演繹推導(dǎo)程序，成功地解決了起潮位調(diào)和展開，最后得到一個包括日、月、行星以及地球形狀影響的行星地球潮汐模型，其中日月潮波3070項。郗欽文教授得到的起潮位展開的新成果已被國際會議推薦用于高精度潮汐數(shù)據(jù)分析，國際上稱為“Xi(郗)引潮位展開”。經(jīng)國外多位學(xué)者檢驗、比較，作為完全調(diào)和展開，無論在時間域還是頻率域上都具有最高的精度，已被德、比、法、美、加等國學(xué)者使用。國際固體潮中心梅爾基奧爾（Melchior）教授曾經(jīng)應(yīng)用該研究成果處理他們的超導(dǎo)重力儀資料，發(fā)現(xiàn)了地核中的慣性重力波。朱盛源和格羅頓（Zhu & Groteu）利用“Xi(郗)引潮位展開”系數(shù)導(dǎo)出了新的天文章動系數(shù)。

郗欽文教授的行星地球潮汐模型——“Xi(郗)引潮位展開”這項研究成果獲得1999年國家自然科學(xué)二等獎。

為了能從清晰記錄到固體潮的觀測中提取震前潮汐波畸變信息，必須推導(dǎo)出重力、傾斜和應(yīng)變固體潮理論值。為了進行維尼柯夫調(diào)和分析，還必須給出二階位與三階位的日波、半日波與1/3日波的大地系數(shù)及其相位改正。1982年，郗欽文在《地球物理學(xué)報》第25卷增刊上發(fā)表了《固體潮理論值計算》，1986年，劉序儼等在《地球物理學(xué)報》第29卷第5期上發(fā)表了《應(yīng)變固體潮理論值計算及調(diào)和分析》，給出了大地系數(shù)與相位滯后，并給出了湖北省黃石地震臺伸縮儀觀測資料的調(diào)和分析結(jié)果。當(dāng)時比利時皇家天文臺國際固體潮中心正與中國地震局地震研究所進行協(xié)作，劉序儼應(yīng)比利時皇家天文臺臺長P.梅爾基奧爾教授邀請，赴天文臺進行學(xué)術(shù)訪問。1987年，張先華與張勇分別在《地殼形變與地震》刊物上發(fā)表了《計算應(yīng)變固體潮理論值的旋轉(zhuǎn)橢球模型》和《地殼潮汐應(yīng)變理論值的計算》。

重力、傾斜與應(yīng)變固體潮是由月亮、太陽起潮力位所引起的三種同源異象現(xiàn)象，它們之間的內(nèi)在聯(lián)系可用來表征地球彈性力學(xué)性質(zhì)的三個勒夫數(shù)（h、l、k）。劉序儼1994年對固體潮觀測資料的維尼迪柯夫調(diào)和分析結(jié)果的地球物理含義進行了分析，指出調(diào)和分析結(jié)果為一復(fù)數(shù)R(w)ejΨ(w),該復(fù)數(shù)等于固體地球的傳遞函數(shù)與地球模型的傳遞函數(shù)之比，該比值的實數(shù)R(w)即為潮汐因子，Ψ(w)為潮汐相位滯后。若地球模型所采用的勒夫數(shù)越接近于固體地球的實際勒夫數(shù)，則R(w)越接近于1，Ψ(w)越接近于零。從這一點出發(fā)，可以說潮汐因子反映了地球模型對固體地球的擬合程度，間接反映了地球的物性變化，并推導(dǎo)出了根據(jù)應(yīng)變固體潮觀測值按最小二乘原理反演勒夫數(shù)h和l的方法。

3.2.3中國地殼形變圖像及其分析方法

上世紀(jì)80年代末，中國地震局開始在全國布設(shè)GPS網(wǎng)，對全國地殼形變進行監(jiān)測。到上世紀(jì)末，利用已積累的GPS觀測資料對中國大陸及鄰區(qū)現(xiàn)今地殼運動和形變進行了分析。1999年，吳云等人在《地震學(xué)報》發(fā)表了《用GPS觀測結(jié)果對中國大陸及鄰區(qū)現(xiàn)今地殼運動和形變的初步探討》一文，該文利用1994年到1996年的全國GPS網(wǎng)21個測點的復(fù)測數(shù)據(jù)，推算了測點的運動速度，采用轉(zhuǎn)動函數(shù)冪級數(shù)展開法，計算了水平速度和應(yīng)變率場，初步探討了中國現(xiàn)今地殼運動和形變的態(tài)勢與特征。結(jié)果表明：印度、太平洋和菲律賓海板塊聯(lián)合控制著中國大陸的現(xiàn)今地殼運動與形變。在三大板塊中，印度板塊起著主導(dǎo)作用，南北地震帶在中國大陸現(xiàn)今地殼運動和形變過程中起著重要的主導(dǎo)作用。到2000年，經(jīng)1992年、1994年及1996年3次GPS會戰(zhàn)，我國用以測量與監(jiān)測中國地殼運動的GPS網(wǎng)業(yè)已建好，其測量地殼運動精度已優(yōu)于5mm/a。這些GPS復(fù)測結(jié)果是聯(lián)合重力測量、水準(zhǔn)測量、地震和地形變觀測手段取得的，足夠建立關(guān)于大尺度變形的定量化模型。以中國科學(xué)院上海天文臺朱文耀為第一作者，包括中國地震局、中國科學(xué)院、國家測繪局、總參等部委等15人在2000年英文版《中國科學(xué)（D）》發(fā)表了“Crustal motion of Chinese mainland monitored by GPS ”一文，分別給出了GPS與NNR-NUVELIA的中國GPS站的水平速度，兩者的水平速度大多數(shù)測點上都大致相同。2001年，王琪等10人在《中國科學(xué)（D輯）》上發(fā)表《中國大陸現(xiàn)今地殼運動和構(gòu)造》一文，該文認為，全球定位系統(tǒng)揭示的中國大陸現(xiàn)今運動場清晰地表明了以活動地塊為單元的分塊運動特征，不同的活動地塊具有不同的運動和變形方式，中國大陸以活動地塊為單元的現(xiàn)今構(gòu)造變形可能與大陸巖石圈的結(jié)構(gòu)和性質(zhì)有關(guān)。2003年，王琪在《地震學(xué)報》（Vol.16,No.5）發(fā)表了《中國大陸現(xiàn)今地殼運動》一文，該文著重闡述了對東亞和中亞的構(gòu)造形變的定量化研究在全球板塊運動和巖石圈動力學(xué)研究中的重要意義。2003年，江在森等在《地球物理學(xué)報》（Vol.46,No.3）發(fā)表《GPS初動結(jié)果揭示的中國大陸水平應(yīng)變場與構(gòu)造變形》一文，該文分析了中國大陸水平運動、應(yīng)變場空間分布特征及其與強震的關(guān)系。2007年，李延興在《地球物理學(xué)報》（Vol.50,No.2）上發(fā)表了《由空間大地測量得到的太平洋板塊現(xiàn)今構(gòu)造運動與板內(nèi)形變場》一文，在該文中，根據(jù)太平洋板塊的空間大地測量觀測結(jié)果，建立了太平洋板塊的彈性運動模型。結(jié)果表明，該模型與板塊實際運動狀態(tài)的符合程度明顯地優(yōu)于剛體運動模型。在太平洋板塊的西北部，主壓應(yīng)變軸為NW—SW方向，主壓應(yīng)變率大于主張應(yīng)變率；在板塊的東南部，主壓應(yīng)變軸為NE—SW方向，主張應(yīng)變率大于主壓應(yīng)變率；板塊的東南邊界為擴張邊界，主張應(yīng)變軸基本上與洋中脊的擴張方向一致；板塊的西北邊界為俯沖邊界，主壓應(yīng)變軸基本上與太平洋板塊的俯沖方向一致。2008年，李延興等人在《地球物理學(xué)報》上發(fā)表了《汶川Ms8.0地震孕育發(fā)生的機制與動力學(xué)》一文。在該文中，作者根據(jù)GPS的多期觀測成果，分析了震前震中及周圍地區(qū)速度場和應(yīng)變場的變化，對汶川大震孕育發(fā)生的機制和動力學(xué)問題進行了探討。根據(jù)中國大陸GPS網(wǎng)2004年和2007年GPS觀測資料處理，獲得巴顏喀拉塊體和四川盆地的應(yīng)變場，發(fā)現(xiàn)震前震中區(qū)的主壓與主張應(yīng)變率分別為-30.840×10-9/a與13.95610-9/a，主壓應(yīng)變軸方向為N105.4°E，與震源機制解得到的主壓應(yīng)力軸的方向N103°E一致。由GPS數(shù)據(jù)得到的發(fā)震斷層龍門山斷裂的斷層滑動方向和走向與地表破裂調(diào)查的震源機制解的結(jié)果一致。

對地應(yīng)變分析方法，武艷強等于2009年在《地球物理學(xué)報》上發(fā)表了《利用最小二乘配置在球面上整體解算GPS應(yīng)變場的方法及應(yīng)用》一文，由于應(yīng)變計算方法的不同，致使不同的處理者即使利用相同的數(shù)據(jù)也很難得到完全相同的結(jié)果，究竟如何評價方法是否有效顯得尤為重要。作者在對最小二乘配置平面解算方法進行深入研究的基礎(chǔ)上，對最小二乘配置在球面上求解應(yīng)變場的方法進行討論，并且在對模擬數(shù)據(jù)進行計算的基礎(chǔ)上討論了該方法的有效性，最后對實際數(shù)據(jù)進行了處理。

3.2.4非線性理論在地殼形變異常分析中的應(yīng)用

根據(jù)巖石破裂實驗，地殼線性變形、非線性變形直至變形局部化現(xiàn)象是與巖石破裂過程的階段性相對應(yīng)的。地形變變化出現(xiàn)的這種階段性是巖石破裂過程在地殼的一種力學(xué)過程的刻劃。當(dāng)然，地殼變形發(fā)生的地點和形式隨每次地震而有所不同，這正是孕震系統(tǒng)具有非線性行為的一種客觀體現(xiàn)，并由斷層模擬研究所證實，用滑塊模型對斷層位移進行模擬研究后發(fā)現(xiàn)，隨著滑塊模型參數(shù)的大小或取值不同，斷層位移會出現(xiàn)復(fù)雜的圖像。劉式達、劉式適（1994）在所著的《孤波與湍波》一書中認為，雖然滑塊的微分方程是確定性的，但隨著參數(shù)的變化，微分方程會出現(xiàn)多解。換句話說，在孕震過程中，由于巖石參數(shù)或初始條件稍微發(fā)生變化，都有可能使斷層位移出現(xiàn)不同的變化。雖然微分方程本身是確定性的，但由于該微分方程含有非線性項，隨著巖石介質(zhì)發(fā)生變化，由該方程表征的斷層位移在其演化過程中會發(fā)生分岔，經(jīng)自組織奔向最后的歸宿吸引子，并不存在一種固定不變的歸宿（中心點）。這正是耗散系統(tǒng)有別于保守系統(tǒng)的一種顯著特點。而孕震過程就是這樣一個具有耗散結(jié)構(gòu)的系統(tǒng)。周碩愚在其所著的《系統(tǒng)科學(xué)導(dǎo)引》一書中對自組織、耗散結(jié)構(gòu)等非線性理論作了較詳細的闡述，在周碩愚此后的一系列論文中，采用非線性理論對強震前后地殼形變場動力學(xué)圖像及其參量特征進行了研究，認為在斷層即將破裂奔向一個新平穩(wěn)態(tài)之前，地殼形變圖像會出現(xiàn)減熵、降維或有序演進趨向，同時由于系統(tǒng)的整體行為，源和場之間的相干與反饋，使近場與遠場地形變圖像在時空中發(fā)現(xiàn)一種有序結(jié)構(gòu)和協(xié)同演進方式，為判別大震前長、中、短、臨前兆提供依據(jù)。

4 福建省地震地形變學(xué)科發(fā)展現(xiàn)狀

4.1 福建省地形變觀測臺網(wǎng)

福建省地形變觀測臺網(wǎng)由固定臺站和流動觀測組成，前者在地震臺站進行連續(xù)觀測，后者不為連續(xù)觀測，而是定期觀測。福建省地形變觀測網(wǎng)是由全省GPS連續(xù)觀測網(wǎng)、重力、水準(zhǔn)、地磁流動觀測網(wǎng)以及重力、傾斜與應(yīng)變觀測臺站組成。

網(wǎng)點的布設(shè)是根據(jù)福建省主要構(gòu)造斷裂而定的，以便這些網(wǎng)點可以監(jiān)測北東向分布的長樂—詔安斷裂帶、政和—海豐斷裂帶與邵武—河源斷裂帶以及北西向分布的閩江斷裂帶、沙縣—南日島斷裂帶、永安—晉江斷裂帶、九龍江下游斷裂帶、上杭—云霄斷裂帶的構(gòu)造活動，在這些斷裂以及北東與北西斷裂交匯部位往往是地震發(fā)生的地點。

4.1.1 GPS網(wǎng)

從上世紀(jì)70年代初，在福建沿海進行了大面積水準(zhǔn)測量、跨斷裂測量，激光測距和重力測量，取得了大量的地殼形變測量數(shù)據(jù)，對研究局部地區(qū)的地殼形變發(fā)揮了重要作用。但由于上述測量手段是常規(guī)的，監(jiān)測范圍分散、邊長短、觀測要求通視、耗時費用高，不可能對福建的大尺度地殼水平運動同時進行高精度的觀測。尤其要想跨過海峽對臺灣島相對于大陸的相對運動進行觀測研究，常規(guī)大地測量方法更是無能為力、只能望洋興嘆。通過多年的論證，認識到采用GPS是克服上述常規(guī)測量方法缺陷最佳的技術(shù)手段。從1993年下半年起，福建省地震局成立了GPS課題組進行調(diào)研論證。于1995年1月向福建省計委有關(guān)部門呈交《臺灣海峽GPS聯(lián)測與臺灣海峽地球動力學(xué)特征研究》項目。同年7月，省計委組織有關(guān)專家對項目進行了開題論證。專家論證后，一致同意該項目的開題報告，建議省計委給予資助，并建議將課題名稱改為《利用GPS技術(shù)開展福建沿海地球動學(xué)特征應(yīng)用研究》。翌年2月，該項目正式立項，簽訂了專題合同書。該項目由劉序儼、林繼華等負責(zé)實施。1995年進行了GPS網(wǎng)優(yōu)化設(shè)計、選點埋石建立了沿海GPS監(jiān)測網(wǎng)。同年12月進行了第一期GPS聯(lián)測，參加聯(lián)測測點的福州、水口、古田、南平、湄州、泉州、德化、廈門、永安、東山、尤溪；1997年12月進行了第二期GPS聯(lián)測；1999年12月進行了第3期聯(lián)測。

臺灣海峽兩岸的GPS網(wǎng)在福建一側(cè)有12個GPS測點，與臺灣地球科學(xué)研究所商定以臺北、臺中、臺南、恒春、臺東、花蓮作為聯(lián)測點，以研究臺灣相對于大陸的碰撞過程，雖然在GPS聯(lián)測期間，兩岸皆進行同步觀測，但由于眾所周知的原因，雙方均不可能擁有兩岸的聯(lián)測數(shù)據(jù)，因此無法取得臺灣海峽地殼水平運動的信息。

在“十五”、“十一五”、“十二五”期間，分別在“福建省地殼形變觀測臺網(wǎng)建設(shè)”、“防震減災(zāi)二期工程建設(shè)GPS連續(xù)觀測臺網(wǎng)建設(shè)”以及“福建省地震監(jiān)測預(yù)警與社會服務(wù)系統(tǒng)之地殼形變觀測臺網(wǎng)建設(shè)”帶動下，福建省地殼形變臺網(wǎng)建設(shè)得到了進一步加強。到目前為止，GPS連續(xù)觀測臺站已達到11個。圖7為福建省GPS連續(xù)觀測臺站分布圖。到“十二五”結(jié)束，福建省連續(xù)GPS觀測站將共有61個。

圖7 福建省GPS連續(xù)觀測臺站分布圖

4.1.2流動GPS觀測網(wǎng)

為了獲取福建省現(xiàn)今水平形變場時空變化特征，福建省地震局廈門地震勘探研究中心建立了一個基本覆蓋全省的GPS基本網(wǎng)，該網(wǎng)共有16個GPS基本觀測點，見圖8。

圖8 福建省流動GPS觀測網(wǎng)

4.1.3流動重力觀測網(wǎng)

為了觀測福建省全省的重力背景場動態(tài)變化及加密跟蹤監(jiān)測重點區(qū)域斷裂構(gòu)造活動區(qū)的重力場的短期變化時空演化特征，我省布設(shè)了覆蓋全省的流動重力監(jiān)測網(wǎng)（見圖9）。該網(wǎng)共有流動重力測點200個，從2008年起每年復(fù)測2期，采用CG5重力儀觀測。

圖9 福建省流動重力觀測網(wǎng)

4.1.4流動地磁觀測網(wǎng)

為監(jiān)測福建省地磁背景場動態(tài)變化特征，我省布設(shè)了基本覆蓋全省的流動地磁監(jiān)測網(wǎng)。該網(wǎng)共布設(shè)有95對流動地磁測點，分為3個磁測區(qū)：東南部沿海區(qū)、中西部測區(qū)和東北部測區(qū)。流動地磁每年復(fù)測2～3次，采用G856磁力儀。流動地磁觀測網(wǎng)布設(shè)見圖10。

圖10 福建省流動地磁觀測網(wǎng)

4.1.5大面積水準(zhǔn)觀測網(wǎng)

為了獲取福建省垂直形變場時空演化特征，2002年布設(shè)了全省垂直形變網(wǎng)，水準(zhǔn)路線布設(shè)在北緯23.6°~26.2°、東經(jīng)117.0°~119.4°范圍，覆蓋福州、莆田、石獅、廈門、漳州、龍巖等城市，構(gòu)成3個水準(zhǔn)閉合環(huán)，共13條水準(zhǔn)路線（其中支線1條），一等水準(zhǔn)1529.7公里，385個測段，水準(zhǔn)點454個，見圖11。

圖11 福建省垂直形變觀測網(wǎng)

廈門地區(qū)進行加密觀測，共13條測線，分島內(nèi)8條和島外5條兩個部分，共245公里，見圖12。

圖12 廈門垂直形變觀測網(wǎng)

4.1.6流動跨斷層短水準(zhǔn)觀測

為了監(jiān)測福建主要斷裂垂直運動，在全省主要斷裂重點部位布設(shè)了17處跨斷層短水準(zhǔn)觀測場地（見圖13）。觀測周期每年3期。采用德國蔡司廠的Ni004精密水準(zhǔn)儀施測。

圖13 福建省跨斷層短水準(zhǔn)觀測場地分布圖

4.1.7跨斷層綜合觀測

為了綜合監(jiān)視長詔斷裂帶北段與中段的斷裂運動特征，分別在福清、泉州各布設(shè)了一個跨斷層綜合場地，集GPS、水準(zhǔn)、重力等觀測手段進行綜合觀測。福清綜合場地觀測路線總長度為69.5 公里，共20個水準(zhǔn)點，10個GPS點、10個重力點。該場地跨過長詔帶的東張—詔安斷裂、長樂—東山斷裂和平潭—東山澳角斷裂。泉州場地觀測路線總長度為69.8公里，共20個水準(zhǔn)點、10個重力點。跨過長詔帶的東張—詔安斷裂、長樂—東山斷裂。從2008年起兩個場地每年復(fù)測1期一等水準(zhǔn)、流動GPS、流動重力。一等水準(zhǔn)采用Ni002精密水準(zhǔn)儀施測，流動GPS采用雙頻大地型GPS接收觀測，流動重力采用CG5流動重力儀觀測。

4.1.8數(shù)字化形變觀測臺站

福建省數(shù)字化形變臺站分布見圖14，各臺站觀測項目和儀器見表1，可以看出，福建省數(shù)字化形變觀測項目比較齊全，但臺站分布不均，主要分布在福建中部和沿海，而北部僅有南平臺。

圖14 福建省數(shù)字化形變觀測臺站分布圖

表1 福建省形變臺站基礎(chǔ)信息表

4.2 福建省地震地形變科學(xué)研究

4.2.1福建沿海地球動力學(xué)特征探討

根據(jù)福建省上述地殼形變觀測分析結(jié)果，發(fā)現(xiàn)福建沿海的地球動力學(xué)有以下特征：

1）福建地區(qū)隨同中國大陸較均勻地向東南方向運動，方向為NW（NWW）~SE（SEE），年速度約為數(shù)毫米，由GPS測定的現(xiàn)時福建沿海的應(yīng)變場與震源機制解給出的區(qū)域應(yīng)力場以及地質(zhì)學(xué)推測的構(gòu)造應(yīng)場三者基本相符，方向皆為NW（NWW）~SE（SEE）。福建沿海與中國大陸地殼運動的這種相似性，說明了它們有同一個力源，即印度板塊、歐亞板塊、太平洋板塊和菲律賓板塊之間的相互作用為這種大尺度的地殼運動提供了力源。

2）由GPS和常規(guī)大地測量給出的福建沿海的地殼形變圖像是南強北弱，東強西弱。由GPS給出的福建沿海的面應(yīng)變剪應(yīng)變都是由北西方向向南西方向呈壓縮狀態(tài)，表明福建沿海的閩南地區(qū)地形變程度較強，這與福建沿海的斷裂活動有關(guān)。因為福建沿海有北東走向的長樂—詔安斷裂帶，該斷裂帶活動比西部的政和—海豐、邵武—河源斷裂帶較活躍，而在閩南沿北西走向的永定—東山、九龍江和永安—晉江斷裂又比閩北的閩江斷裂帶活躍，因此福建沿海的構(gòu)造形變呈現(xiàn)出南強北弱、東強西弱的特征。

3）福建沿海的地殼形變南強北弱、東強西弱的特征與該地區(qū)的地震活動特征吻合，說明斷裂活動是造成構(gòu)造形變的主要原因。因此，采用GPS技術(shù)可以迅速而精確地監(jiān)測到地震孕育過程中地殼變形演化過程，從而可以為地震的中長期趨勢提供科學(xué)依據(jù)。

4）根據(jù)重力和人工地震資料分析表明，福建沿海地區(qū)的地殼結(jié)構(gòu)具有由陸殼向洋殼過渡的特征。

5）造成福建沿海地殼運動和臺灣海峽、臺灣及東部海域大震的動力來源于我國西部受到來自烏茲別克斯坦、帕米爾以及青藏高原向東的推動力以及東亞大陸與臺灣島東部同時受到太平洋板塊和菲律賓海板塊向西的推動力聯(lián)合作用的結(jié)果。

福建沿海地球動力學(xué)特征這一科研初步成果是1996年省計委列為全省重點項目《利用GPS技術(shù)開展福建沿海地球動力學(xué)特征應(yīng)用研究》中的主要成果，該項目獲2001年福建省科學(xué)技術(shù)進步三等獎。

圍繞福建沿海地球動力學(xué)特征這一項目，福建省地形變科技人員共撰寫了20篇論文，有的在《地殼形變與地震》上發(fā)表，有的在有關(guān)國際學(xué)術(shù)會議上宣讀，有的在海峽兩岸測繪學(xué)術(shù)研究會上宣讀。例如劉序儼、林繼華等撰寫的《福建沿海地球動力學(xué)特征淺析》、梁全強等《利用GPS技術(shù)研究福建省孕震構(gòu)造網(wǎng)絡(luò)活動特征》、林繼華等《福建沿海平面變化及其影響》、劉序儼等《福建省現(xiàn)代地殼垂直運動與斷裂活動》等論文發(fā)表在1999年《地殼形變與地震》增刊Ⅱ，陳園田等《臺灣海峽兩岸新生代地殼運動的對比》發(fā)表在《臺灣海峽》1999年第16卷第3期上，林繼華等《利用GPS技術(shù)開展福建沿海與臺灣海峽地球動力學(xué)特征研究》一文在由臺灣成功大學(xué)主辦的第二屆測繪學(xué)術(shù)研究會暨第十七屆測量學(xué)術(shù)及應(yīng)用研討會（1998年）上宣讀，劉序儼、郭逢英、林繼華撰寫的《The Approach of Geodynamic Features and Research Methods in Fujian Coast of South China and Taiwan Strait》一文在1997年在中國武漢市召開的“東亞與東南亞現(xiàn)今地殼運動和減災(zāi)國際學(xué)術(shù)會議”上宣讀。

4.2.2地形變基礎(chǔ)理論的譯著

自福建省地震局成立以來，除福建省建立起一個能覆蓋全省地殼形變的臺網(wǎng)以外，從事地形變的科技人員在做好日常的數(shù)據(jù)處理工作以外，還追蹤了國內(nèi)外地形變動態(tài)，及時掌握前沿知識，努力做到博采眾長、為我所用、有所創(chuàng)新。上世紀(jì)80年代初起，美國開始應(yīng)用空間技術(shù)與地面技術(shù)致力于構(gòu)造形變之研究，以探明震前的長、中、短、臨地形變異常。1983年劉序儼等翻譯了美國地殼運動測量委員會撰寫的一份有關(guān)美國在上世紀(jì)80年代10年間的現(xiàn)代地殼運動尤其是地震活動帶的監(jiān)測工作提出的戰(zhàn)略和規(guī)劃，書的譯名為《構(gòu)造形變》，該書于1983年由地震出版社出版。1993年，劉序儼與徐菊生等人翻譯了德國重力學(xué)家Wolfgaug Torge 所著的《Gravimetry》一書，書的譯名為《重力測量學(xué)》，該書由地震出版社出版。該書包括了重力測量學(xué)的基礎(chǔ)概念、基礎(chǔ)理論、應(yīng)用研究、重力觀測方法和資料處理及其發(fā)展簡史。該書還特別安排了專門的章節(jié)對現(xiàn)今重力測量學(xué)研究領(lǐng)域的重點課題“重力隨時間變化的測量與解釋”，進行了精辟地闡述。1997年，劉序儼等編寫了《地殼形變分析預(yù)報方法》一書，該書系統(tǒng)總結(jié)和全面闡述了有關(guān)地形變預(yù)報地震的科學(xué)原理、預(yù)報思路，各類地形變數(shù)據(jù)處理及其預(yù)報方法，匯集了國家地震局“七五”和“八五”期間地形變攻關(guān)成果中的精華，以及當(dāng)時國際上在這一研究領(lǐng)域的前沿性最新研究成果，并集中反映了該書作者的一些科學(xué)見解。該書作為當(dāng)時國家地震局地震預(yù)報系列教材之一，于1998年由地震出版社出版。1991年，劉序儼參加了國家地震局科技監(jiān)測司委托中國地震學(xué)會地殼形變測量專業(yè)委員會會同科技監(jiān)測司聯(lián)合編著的《地殼形變監(jiān)測整體規(guī)劃》的有關(guān)編寫工作，該《規(guī)劃》重點提出了在這一領(lǐng)域中當(dāng)前國際上一些突出的前沿課題。

4.2.3地形變基礎(chǔ)理論研究

福建省地震局從事地震地形變工作的科技人員對地形變基礎(chǔ)理論作了較深入的研究，從1989年至今發(fā)表了50余篇論文，其中主要有：劉序儼等在2002年第4期《大地測量與地球動力學(xué)》上發(fā)表的《平動及其在起潮力中的作用》一文，重新對天體起潮力作了新定義。經(jīng)過嚴(yán)密的數(shù)理推導(dǎo)，給出了新定義：“地球上某一點的天體起潮力為天體對該點的萬有引力與天體對地心的萬有引力之矢量差”。該定義比原定義“地球上某一點的天體起潮力為天體對該點的萬有引力與該點隨同地心繞地球與天體的公共質(zhì)心作圓周運動的慣性離心力的矢量和”顯得更簡潔，并摒棄了慣性離心力之提法，更揭示了起潮力之本質(zhì)。在固體潮整體研究方面，劉序儼等發(fā)表在《地殼形變與地震》（1997年）上的《固體潮整體綜合分析方法研究》一文給出了如何利用重力、傾斜與應(yīng)變固體潮觀測數(shù)據(jù)進行綜合分析的方法。在2003年第2期《大地測量與地球動力學(xué)》上發(fā)表的《論潮汐波與地震波》一文對潮汐波與地震波的異同及兩者的特征進行了較全面深入的闡述。劉序儼等在2005年第1期《大地測量與地球動力學(xué)》上發(fā)表了《地面傾斜的形變特征及計算方法》一文，該文從理論上論證地面傾斜是一種偽矢量，給出了根據(jù)任何2個方向的地傾斜分量計算地面傾斜的普適表達式，并給出了描述地傾斜面的數(shù)學(xué)表達式，同時指出地傾斜面僅是一個相對平面，并不是唯一的。在該刊2007年第4期上發(fā)表了《潮汐傾斜的幾何物理含義及其剪應(yīng)變屬性》一文，該文給出了潮汐傾斜的幾何物理含義，并對為什么剪應(yīng)變可作為其屬性進行了闡述。在《地震學(xué)報》2007年第3期的中英文刊物上發(fā)表了《旋轉(zhuǎn)橢球面上的應(yīng)變與轉(zhuǎn)動張量表達》一文，該文推導(dǎo)出了在大地坐標(biāo)系中的應(yīng)變與轉(zhuǎn)動張量表達式，可供GPS觀測數(shù)據(jù)在全球大地參考系（WGS84）中作應(yīng)變與轉(zhuǎn)動張量分析之用。在目前的彈性力學(xué)有關(guān)的文獻中僅給出了在直角坐標(biāo)系、柱坐標(biāo)系、球坐標(biāo)系下的應(yīng)變與轉(zhuǎn)動張量表達式，缺失橢球坐標(biāo)下的表達式，該文給出的橢球面上的表達式則可填補這方面的缺失。翌年2月在《大地測量與地球動力學(xué)》上發(fā)表了《正交曲線坐標(biāo)系的應(yīng)變張量轉(zhuǎn)換》一文，給出了在直角坐標(biāo)系、柱坐標(biāo)系、球坐標(biāo)系、旋轉(zhuǎn)橢球坐標(biāo)系以及旋轉(zhuǎn)拋物線等坐標(biāo)系相互之間的應(yīng)變張量轉(zhuǎn)換表達式。同年4月，又在同一刊物上發(fā)表了《正交曲線坐標(biāo)系應(yīng)變張量的普適表達》一文，經(jīng)嚴(yán)格推導(dǎo)，給出了在任何一種正交曲線坐標(biāo)系下的應(yīng)變與轉(zhuǎn)動張量普適表達式，這是在考慮到任何一種正交曲線坐標(biāo)系下的一點的單位標(biāo)架為正交標(biāo)架的前提下，以已知的直角坐標(biāo)系的應(yīng)變與轉(zhuǎn)動張量矩陣為平臺的基礎(chǔ)上推導(dǎo)出來的，推導(dǎo)簡潔明了。2010年、2011年，劉序儼等先后在《大地測量與地球動力學(xué)》上發(fā)表了《地形變旋轉(zhuǎn)張量探討》與《地形變應(yīng)變張量矩陣的不變量分析》兩篇論文，分別對旋轉(zhuǎn)與應(yīng)變張量矩陣的幾何物理特征及其所蘊藏的幾何物理不變量進行了較深入的闡述。李華等在2008年第2期《大地測量與地球動力學(xué)》上發(fā)表了《漳州臺地傾斜矢量圖的潮汐導(dǎo)納分析》一文，該文利用地傾斜觀測與理論矢量圖，對兩者在同一時刻的矢量圖的矢經(jīng)與方位角進行了數(shù)值對比計算，從而取得了潮汐因子與相位滯后值，發(fā)現(xiàn)它們與觀測值的維尼迪柯夫調(diào)和分析結(jié)果較一致。季穎峰等在2010年第2期《華南地震》上發(fā)表了《大地地坐標(biāo)系應(yīng)變張量表達及其與地心直角坐標(biāo)系的相互轉(zhuǎn)換》一文，給出了WGS84坐標(biāo)系與地心直角坐標(biāo)系的具體應(yīng)變張量矩陣轉(zhuǎn)換式。在2011年的《大地測量與地球動力學(xué)》增刊上發(fā)表了黃聲明的《坐標(biāo)系的物理基礎(chǔ)及其空間結(jié)構(gòu)》一文，對坐標(biāo)系的物理基礎(chǔ)及其空間結(jié)構(gòu)作了深刻的闡述。

眾所周知，地形變測量可以給出地殼是如何變形的，但想要知道地殼變形的力源是什么，還要從地球物理與地質(zhì)學(xué)角度進行探討，福建省地震局李祖寧高級工程師在這方面進行了探討。在2002年第45卷第4期的《地球物理學(xué)報》上發(fā)表了《剝蝕及地幔作用下青藏高原隆升過程的數(shù)值模擬》一文，在該文中，作者修改了England和Mckenzie 的黏性薄層流變模型中控制大陸形變的連續(xù)性方程，將剝蝕作用對高原隆升演化的影響直接引入該方程，并考慮下伏地幔小尺度對流對增厚巖石層的搬離作用對高原隆升演化后期的影響，用有限差分法直接模擬青藏高原隆升過程。數(shù)值模擬結(jié)果所顯示的高原隆升演化過程與實際觀測資料吻合較好，揭示了高原隆升過程的非平穩(wěn)和多階段的特性，同時還表明上地幔小尺度對流對巖石層底部的搬離作用可能是最近8Ma 以來高原快速隆升的主導(dǎo)機制。李祖寧等在同年《地震學(xué)報》上發(fā)表了《多種驅(qū)動力作用下東亞大陸形變及應(yīng)力場演化》一文，該文認為控制東亞大陸西部形變和應(yīng)力場基本格局的主要驅(qū)動力來源于印度板塊對歐亞板塊的碰撞、擠壓，而對東部地區(qū)還應(yīng)當(dāng)考慮其與太平洋板塊和菲律賓海板塊的相互作用。2005年，李祖寧在《地震地質(zhì)》上發(fā)表了題名為《利用數(shù)字地震臺網(wǎng)資料聯(lián)合反演福建地區(qū)Q值、場地響應(yīng)和震源參數(shù)》一文，該文利用福建省數(shù)字地震臺網(wǎng)1998年建網(wǎng)以來記錄到的2.5級以上地震的波形資料，采用Atkinson 等提出的方法計算了福建地區(qū)地震波非彈性衰減系數(shù)、幾何擴散系數(shù)和臺站場地響應(yīng)，對震源譜的低頻水平和拐角頻率進行了聯(lián)合反演，并取得了有關(guān)地震的地震矩、應(yīng)力降和破裂半徑，并取得了福建地區(qū)地殼平均Q值與頻率的關(guān)系式，而場地響應(yīng)表現(xiàn)為與頻率相關(guān)并基本上在1附近波動，這與福建數(shù)字地震臺站基本上建立在基巖之上的實際情況相符合，震源參數(shù)的計算結(jié)果表明，地震矩與震級以及震源譜拐角頻率的相關(guān)性較好。

4.2.4地形變數(shù)據(jù)處理及提取地震前兆信息研究

在取得了大量地形變觀測資料后，如何進行數(shù)據(jù)處理、排除干擾、去偽存真，從中提取地殼構(gòu)造形變信息，是擺在科技人員面前的重大課題，為此，福建省科技人員作了不懈的努力和探討。

在用地殼形變觀測資料對福建地殼形變的時空變化特征研究方面，郭逢英等撰寫了《GPS測量揭示的福建沿海地區(qū)的現(xiàn)今地形變》（發(fā)表在2003年第3期《大地測量與地球動力學(xué)》），丁學(xué)仁等撰寫的《福建地區(qū)地塊運動的時空變化特征研究》（發(fā)表在2006年第2期《大地測量與地球動力學(xué)》），梁全強等撰寫的《臺灣海峽西岸重力場變化梯度場圖像分析》（發(fā)表在2006年第1期《大地測量與地球動力學(xué)》），陳光等撰寫的《福建地區(qū)地殼形變特征研究》（發(fā)表在2009年第3期《高原地震》），鐘繼茂撰寫的《福建近岸海域濱海斷裂的應(yīng)力分析》（發(fā)表在2009年第5期《地震地磁觀測與研究》），陳光等的《福建地殼運動與地震關(guān)系初探》（發(fā)表在2010年第1期《華北地震科學(xué)》），上述這些論文，從不同時期的地形變觀測資料對福建地區(qū)的現(xiàn)今形變進行了分析，除不同年份的地殼形變速率稍有變化外，總的福建內(nèi)陸運動狀態(tài)仍為南強北弱，東強西弱，因為地殼運動狀態(tài)基本上為活動斷裂帶所制約。

除對福建內(nèi)陸地殼形變的時空演變狀態(tài)進行分析以外，利用地形變資料對福建水口的應(yīng)力分布與地震活動性進行了探討，鐘繼茂發(fā)表在2009年第3期《華南地震》的《福建水口庫區(qū)次級斷裂的應(yīng)力分析和活動特征研究》一文，除選用庫區(qū)15個P波初動清晰可靠的地震震源機制解推算了庫區(qū)構(gòu)造應(yīng)力場的主軸方向外，還根據(jù)庫區(qū)布設(shè)的6個跨斷層短水準(zhǔn)場地觀測資料，通過線性擬合發(fā)現(xiàn)西瓜洲斷裂為正斷活動，其余為逆斷活動；陳光等撰寫的《水口水庫區(qū)地震活動特性分析》（發(fā)表在2010年《大地測量與地球動力學(xué)》增刊），該文根據(jù)庫區(qū)的短水準(zhǔn)地形變資料與震源機制解結(jié)果進行了分析，分析結(jié)果表明，水口庫區(qū)的地震活動前期主要為水的誘發(fā)作用，后期與構(gòu)造應(yīng)力場作用加強有關(guān)。

在對地形變觀測資料提取前兆信息方面，也作了很多探討。例如在利用小波分析方法對地形變觀測資料中所隱藏的豐富信息進行多層次分析方面，劉序儼、鐘繼茂、陳光分別采用不同的地形變資料進行了分析，劉序儼等撰寫的《定點形變數(shù)字觀測中高頻信號特征及提取方法》（發(fā)表在2006年第6期《大地測量與地球動力學(xué)》）一文根據(jù)廈門臺傾斜與伸縮儀的分鐘值進行從周期為2分鐘到超低頻這樣一個寬廣的頻域內(nèi)對信息進行了多層分析，分析結(jié)果表明，在大震前會有高頻信息出現(xiàn)；鐘繼茂撰寫的《基于小波分析的閩臺地區(qū)地震形變異常研究》（發(fā)表在2011年第1期《華北地震科學(xué)》）一文，采用小波分析，發(fā)現(xiàn)在強震前2~3年出現(xiàn)周期2年左右的異常信號，在中強地震前2~5個月出現(xiàn)周期1個月左右的異常信號，在震前幾天到十幾天則出現(xiàn)周期半小時到1小時的異常信號；陳光等在利用福建GPS網(wǎng)2004年3月~2008年10月的連續(xù)觀測資料所得到的位移時間序列后，利用小波分析方法提取時間序列的時頻特征，從而得到了非線性變化信息。

除了采用小波分析方法提取前兆信息外，還探討了利用地震震源解、潮汐應(yīng)變場前兆圖像法、主應(yīng)變提取前兆方法。例如鐘繼茂的《由多個地震震源機制解求川滇地區(qū)平均應(yīng)力方向》論文（發(fā)表在《地震學(xué)報》2006年第4期），劉序儼等撰寫的《天馬和謝坊形變試驗場地殼形變特性分析》（發(fā)表在1993年第1期《地殼形變與地震》上），以及發(fā)表在1994年《地球物理學(xué)報》增刊Ⅱ的《應(yīng)變固體潮主應(yīng)變的計算》，鐘繼茂的《區(qū)域地震狀態(tài)的損傷張量分析》（《中國地震》，2005年第21卷第3期）。

長期在臺站從事的科技人員，除搞好儀器維護以外，還開展了科研工作，發(fā)表了不少論文。例如廈門臺的熊先保高級工程師先后在《大地測量與地球動力學(xué)》刊物上發(fā)表了多篇論文，例如《廈門臺地傾斜觀測資料中海潮干擾的綜合特征》、《海潮對廈門臺地傾斜干擾的可能機制》、《廈門臺數(shù)字化潮汐資料背景變化及其主要干擾因素》，以及分別發(fā)表在《華南地震》、《地震研究》、《地震地磁觀測與研究》上的《廈門臺地形變觀測資料中的“海潮”效應(yīng)》、《從三次地震來討論廈門臺地形變儀器的監(jiān)測能力》、《固體潮觀測數(shù)據(jù)預(yù)處理中的一些常用方法及計算機實踐》。廈門臺的黃曉華高級工程師對臺站的地傾斜與應(yīng)變觀測資料的異常進行了深入的分析，先后在《大地測量與地球動力學(xué)》與《華南地震》刊物上發(fā)表了4篇論文，分別是《廈門臺地震前兆異常分析》、《廈門地震臺數(shù)字化水管傾斜儀觀測資料分析》、《數(shù)字化形變觀測典型曲線畸變分類與說明》、《臺灣海峽地震定點形變的異常分析》。廈門臺的楊婕工程師從2008年至今分別在《防災(zāi)科技學(xué)院學(xué)報》、《國際地震動態(tài)》、《地震地磁觀測與研究》、《華南地震》、《華北地震科學(xué)》、《大地測量與地球動力學(xué)》等刊物上發(fā)表了10篇論文，論文集中討論了地形變儀器性能、故障分析與處理、強降雨干擾以及靜電干擾影響及預(yù)防、同震響應(yīng)與異常圖像分析處理等問題，對如何保證地形變儀器正常運轉(zhuǎn)具有參考價質(zhì)。劉其壽在龍巖地震臺從事地傾斜管理工作，他除了維護儀器正常運轉(zhuǎn)以外，同時還從事觀測數(shù)據(jù)的處理與軟件編程工作，先后在多家刊物上發(fā)表了許多論文，例如在《華南地震》上發(fā)表了《測震臺址勘選數(shù)據(jù)處理輔助程序設(shè)計與應(yīng)用》（2010年第2期）、在《大地測量與地球動力學(xué)》（2010年第6期）上發(fā)表《龍巖地震臺中強地震前地傾斜異常分析》等論文。

李祖寧等利用1995年~2003年福建GPS測站與臺北和桃園核心基準(zhǔn)站資料，計算了福建GPS各測站和桃園之間的邊長及其相對變化，計算了臺北和桃園測站的主應(yīng)變和方位角。計算結(jié)果表明，在1995年~1997年和2001年~2003年間臺海兩岸的地殼相對運動都處于拉張狀態(tài)，但兩期拉張的主應(yīng)變大小和方向有較大變化，這種變化可能反映了1999年臺灣集集7.6級地震前后應(yīng)變積累與釋放過程，此研究成果發(fā)表在2007年第5期《大地測量與地球動力學(xué)》刊物上，其論文名為《臺灣海峽兩岸地殼相對運動分析》。該論文曾獲得福建省科協(xié)優(yōu)秀論文三等獎。

4.2.5綜合研究

除地形變觀測以外，與地形變有關(guān)的觀測系統(tǒng)還有井水位與測震儀。從本質(zhì)上講，地形變與地震皆為地殼質(zhì)點的一種振動。地震發(fā)生后，井水位能記錄到所謂的水震波，體應(yīng)變儀與伸縮儀也能記錄到由地震所激勵的體應(yīng)變，同樣，地震儀能記錄到由震源發(fā)出的地震波傳播到觀測地點的瑞雷波。體應(yīng)變儀、伸縮儀、地震儀與井水位記錄到的同震響應(yīng)之間到底有什么關(guān)系？又如何進行綜合分析呢？對此，劉序儼等進行了研究。劉序儼、鄭小菁等認證了承壓井水位記錄到的所謂“水震波”實際上為體應(yīng)變波。該種波可由井水位觀測值乘以其格值而得到。井水位的格值可由井水位觀測資料的維尼迪柯夫調(diào)和分析結(jié)果的振幅比（潮汐因子）取得。他們推導(dǎo)了承壓井水位觀測系統(tǒng)對體應(yīng)變的響應(yīng)機制，給出了井水位與體應(yīng)變關(guān)系的表達式，該研究成果發(fā)表在2009年第12期的《地球物理學(xué)報》上，論文名為《承壓井水位觀測系統(tǒng)對體應(yīng)變響應(yīng)機制分析》。由地震儀記錄到的瑞雷波也可計算出體應(yīng)變，因為體應(yīng)變等于瑞雷波的位移與哈密爾頓微分矢量算子的矢量積。劉序儼等根據(jù)福建5口井水位的觀測數(shù)據(jù)與5臺地震儀的瑞雷波記錄計算了福建地塊對汶川大震的體應(yīng)變響應(yīng)，并進行了對比分析，其研究成果發(fā)表在2009年《大地測量與地球動力學(xué)》，論文名為《福建井水位對汶川大震體應(yīng)變響應(yīng)》。由伸縮儀觀測也可求出體應(yīng)變，因為相互正交方向上伸縮儀的觀測值之和的2/3即為體應(yīng)變。在劉序儼指導(dǎo)下，陳瑩根據(jù)井水位、體應(yīng)變、伸縮儀與地震波觀測資料對福建地塊對汶川大震的體應(yīng)變響應(yīng)進行了計算分析，撰寫了兩篇論文，分別為《福建省井水位、體應(yīng)變與伸縮儀對汶川大震的響應(yīng)》（2010年《大地測量與地球動力學(xué)》）、《汶川大震對福建四種觀測系統(tǒng)的體應(yīng)變響應(yīng)分析》（發(fā)表在2009年《華南地震》）。

此外，劉序儼還對地形變、數(shù)字觀測原理以及觀測系統(tǒng)的逼真度問題進行了探討，先后撰寫了兩篇論文，分別為《地形變數(shù)字觀測基礎(chǔ)》（發(fā)表在2002年第2期《大地測量與地球動力學(xué)》）、《觀測系統(tǒng)的逼真度特性分析》（《大地測量與地球動力學(xué)》2010年增刊）。并先后參與了中國地震局2001年頒布的《地殼形變數(shù)字觀測技術(shù)規(guī)范》、國標(biāo)GBT19531.3-2004地震臺站觀測環(huán)境技術(shù)要求（地殼形變）、中國地震局科技司編寫的《地震地形變觀測技術(shù)》與2003年編寫的《地殼形變數(shù)字觀測技術(shù)》兩書的審核與有關(guān)內(nèi)容的編寫工作。

5 福建省地震地形變學(xué)科發(fā)展存在的問題

5.1 觀測環(huán)境保護

隨著國民經(jīng)濟的高速發(fā)展，各地修建高速公路、高鐵、住宅以及水利工程，都有可能對GPS、水準(zhǔn)、重力測點以及臺站的觀測儀器產(chǎn)生干擾甚至損壞，保護好地形變觀測環(huán)境是一個值得重視的問題。

5.2 外界干擾

地形變觀測儀器精度高，對外界各種干擾（例如人工爆破、氣溫、氣壓、降水等）極其靈敏，如何消除或降低干擾是提取與識別前兆信息的一個非常重要的環(huán)節(jié)。

5.3 儀器老化

臺站和用以流動觀測的儀器，有的老化現(xiàn)象嚴(yán)重，對取得連續(xù)可靠的觀測數(shù)據(jù)造成嚴(yán)重影響。

5.4 人員培訓(xùn)

由于上世紀(jì)70～80年代開始從事地形變工作的老科技人員陸續(xù)退離工作崗位，新上崗的年輕科技人員業(yè)務(wù)、技能、對《觀測規(guī)范》的了解熟悉不夠，對各種地形變儀器的原理、維修、管理以及數(shù)據(jù)處理較生疏，對這些年輕科技人員加強培訓(xùn)工作應(yīng)擺在議事日程上。

5.5 地震預(yù)報要知難而進

雖然目前還未找到根治癌癥的良方，但醫(yī)生仍需懸壺濟世、救治病人，地震宛如是地球的“癌癥”，雖然目前地震預(yù)報遠未過關(guān)，地形變用以捕捉地震前兆的理論還處于探索階段，地震地形變學(xué)科任重道遠，但地震預(yù)報仍要知難而進。

6 福建省地震地形變學(xué)科未來發(fā)展重點

針對上述存在問題，前四個是一些較具體的技術(shù)性問題。通過具體措施并與臺站所在地有關(guān)部門協(xié)作可較好解決環(huán)境保護和外界干擾問題；對于儀器老化問題，省局有關(guān)部門正著手對儀器維修、標(biāo)定與更新作出規(guī)劃；對人員培訓(xùn)，省局人事處、監(jiān)測處與科技處正聯(lián)手制訂人員的培訓(xùn)計劃，對新招聘的科技人員，對博士與碩士優(yōu)先錄用且要進行嚴(yán)格的面試，采取這一系列的措施，有望提高在職科技人員的業(yè)務(wù)素質(zhì)與水平。在解決處理好這四個問題的基礎(chǔ)上，重點是如何解決好地震地形變學(xué)科知難而進的問題。地震預(yù)報是個世界性難題，不是短期內(nèi)能解決的，必須要有長期的戰(zhàn)略性眼光。根據(jù)目前福建省地震局地震地形變學(xué)科的現(xiàn)狀，在今后幾年內(nèi)應(yīng)加強以下幾方面的分析研究工作：

6.1 加強地形變預(yù)報地震機理的研究

6.1.1地形變異常產(chǎn)生的機理。目前已有不少地形變異常機理，實質(zhì)上大都屬于某種科學(xué)假設(shè)。為此，仍然需要深入研究地形變的前兆異常機理問題，提出能夠解釋實際觀測到的異常變化的新理論，力求理論與實際相結(jié)合。

6.1.2地形變時空變化與地殼構(gòu)造活動的關(guān)系。地殼現(xiàn)今構(gòu)造形變具有多種形式，一方面要研究在震源過程中地形變的演變形式，另一方面要研究在地震孕育、發(fā)展以至發(fā)生的各階段地形變的時空變化特征及其指標(biāo)。

6.1.3地殼應(yīng)力應(yīng)變關(guān)系。地殼形變是個力學(xué)問題，應(yīng)力與應(yīng)變是個因果關(guān)系，但又不是一個線性關(guān)系，地形變是在孕震過程中應(yīng)力積累直至斷層發(fā)生破裂的一種表現(xiàn)形式，因此要研究孕震過程中地殼應(yīng)力的演變圖像，以及伴隨應(yīng)力所產(chǎn)生的應(yīng)變圖像，著重要探討如何根據(jù)地形變數(shù)據(jù)反演地殼應(yīng)力。

6.2 開展多測項和多學(xué)科的綜合研究。地形變的各種測項以及其他學(xué)科的多種測項的前兆異常都是對震源應(yīng)力場、區(qū)域應(yīng)力場的不同響應(yīng)，他們都是同源異?，F(xiàn)象，在成因上都有內(nèi)在聯(lián)系，研究不同學(xué)科的各測項的內(nèi)在成因關(guān)系，找出表征這種內(nèi)在成因關(guān)系的定量指標(biāo)有利于提高預(yù)報的可信度，有利于識別異常，因此除了要開展地形變不同測項的綜合研究外，還要開展地形變與地下流體、地震波、地應(yīng)力、重力等前兆異常成因關(guān)系的研究。

6.3 開展福建沿海及海峽地殼運動模型與變形分析。在多期GPS觀測成果的基礎(chǔ)上，采用連續(xù)體變形理論建立經(jīng)緯格網(wǎng)的福建沿海與臺灣海峽的地殼運動模型，以便根據(jù)有限個GPS測點的地殼運動速率來計算未設(shè)測點的其他區(qū)域的地殼運動速率，給出經(jīng)緯格網(wǎng)的線應(yīng)變、最大剪應(yīng)變、面應(yīng)變與轉(zhuǎn)動速率的等值線圖以及主應(yīng)變等值線圖。

6.4 福建沿海地殼垂直形變及其特征分析。地殼垂直運動差異性顯著地區(qū)往往為構(gòu)造塊體相互作用的交接帶，找出這些垂直運動差異性顯著地帶，對判斷地形變前兆具有參考價值。

6.5 福建沿海地球動力學(xué)特征探討。根據(jù)地殼形變觀測數(shù)據(jù)探討產(chǎn)生形變的力源，以達到了解構(gòu)造應(yīng)力分布及積累情況，促進對地震中、長、短、臨變化的認識。

6.6 臺灣海峽兩岸地殼相對運動速率及方向研究。根據(jù)《國務(wù)院關(guān)于支持福建省加快建設(shè)海峽西岸經(jīng)濟區(qū)的若干意見》精神，福建省地震局將加大和臺灣地震研究部門的合作力度，有望在不久的將來共同就兩岸GPS聯(lián)測進行深度合作和交流，屆時將根據(jù)兩岸的GPS同步觀測資料對兩岸的相對地殼運動速率及其方向進行研究?？梢灶A(yù)料的是，這塊由于眾所周知原因尚未開墾的地學(xué)處女地，屆時將成為國內(nèi)外進行地球動力學(xué)研究的熱點地區(qū)，不但可以推進和加深兩岸GPS聯(lián)測與合作研究，而且有利于進一步提升兩岸地震地形變科技水平。

[1] [美]地殼運動測量委員會. 劉序儼,鄒其嘉, 譯.構(gòu)造形變測量[M]. 北京: 地震出版社,1984.

[2] 國家地震局科技監(jiān)測司. 地震地形變觀測技術(shù)[M]. 北京: 地殼出版社, 1995.

[3] 中國地震局監(jiān)測預(yù)報司. 地殼形變數(shù)字觀測技術(shù)[M]. 北京: 地震出版社, 2003.

[4] 福建省地震局分析預(yù)報中心GPS課題組.利用GPS技術(shù)開展福建沿海地球動力學(xué)特征應(yīng)用研究,1992.

[5] 張寶紅.美國板塊邊界觀測(PBO)中的鉆孔應(yīng)變觀測設(shè)備[J]. 大地測量與地球動力學(xué),2010, 30(ZⅡ):81-85.

[6] 周碩愚,施順英,吳云,等.強震前后地殼形變場動力學(xué)圖像及參量特征研究[J]. 地震學(xué)報, 1998, 20(1):41-47.

[7] 江在森,馬宗晉,張希,等.GPS初步結(jié)果揭示的中國大陸水平應(yīng)變場與構(gòu)造變形[J]. 地球物理學(xué)報,2003, 46(3):352-358.

[8] WangQi. Current Crustal Movement in Chinese Mainland[J]. ACTA SEISMOLOGICA SINICA, 2003, 16(5).

[9] Morgan T. Page. Susana Custodio , Ralph J. Archuleta and J.M.Carlson, Constraining earthquake Source inversions with GPS date ,Journal of Geophysical Research, 2009,114.

[10] 劉序儼,鄭小菁,王林,季穎鋒. 承壓井水位觀測系統(tǒng)對體應(yīng)變響應(yīng)機制分析[J]. 地球物理學(xué)報, 2009, 52(12):3147-3157.

[11] 李祖寧,劉序儼,吳紹祖,陳光. 臺灣海峽兩岸地殼相對運動分析[J]. 大地測量與地球動力學(xué),2007,27(5):18-21.

[12] 劉序儼,楊軍. 地殼形變分析預(yù)報方法[M]. 北京: 地震出版社,1998.

課題組成員：

1、劉序儼，福建省地震局研究員；

2、鄭小菁，福建省地震局，高工；

3、王 林，福建省地震局工程師；

4、李祖寧，福建省地震局，高工；

5、梁全強，福建省地震局，總工、高工。