李江俊

(呂梁學(xué)院汾陽師范分校，山西 汾陽 032200)

時至今日，對重力學(xué)的研究已有400年的歷史[1]。有學(xué)者在16世紀(jì)初就發(fā)現(xiàn)地球存在重力。有學(xué)者在17世紀(jì)發(fā)現(xiàn)，重力在不同的位置存在一定的差別。在上述研究理論基礎(chǔ)下，有學(xué)者提出了地球重力變化產(chǎn)生的原因：地球兩極扁平，促使重力由兩極向赤道增大，從而導(dǎo)致不同位置出現(xiàn)不同程度的重力，并提出水準(zhǔn)面概念[2]。20世紀(jì)初，眾多科學(xué)家投入重力研究中，并不斷完善重力網(wǎng)，采用地球衛(wèi)星采集地球重力數(shù)據(jù)，健全地球重力模型[3]。隨著地球重力模型的完善，有學(xué)者針對地球重力異常部分進(jìn)行了研究，并采用動力學(xué)對重力變化的形變和運動進(jìn)行推斷[4]。由于已知地球地質(zhì)、板塊運行等地球地質(zhì)信息，只需要根據(jù)地球重力場數(shù)據(jù)，通過統(tǒng)一歸算改正，就可以得到重力改正、重力異常等重力場信息，為此提出重力場深部結(jié)構(gòu)與動力學(xué)特征分析這一課題。

1 分析重力場深部結(jié)構(gòu)

此次分析重力場深部結(jié)構(gòu)，將采用反演的方式，分析重力場深部結(jié)構(gòu)，形成重力場深部莫霍面結(jié)構(gòu)圖。

1.1 計算重力場場源邊界

計算重力場邊界，需要計算重力場局部梯度，確定重力場中心和邊緣位置，從而確定重力場場源邊界。因此，假設(shè)重力場數(shù)據(jù)獲取點為x、y、z三點，其重力場場源點則為x0、y0、z0三點，則重力場水平梯度的最大值max{T(x，y)，θ}為

(1)

式中：?x和?y為重力場數(shù)據(jù)獲取點x、y兩個方向上的導(dǎo)數(shù)；θ表示重力場水平梯度方向[5]。此時，根據(jù)式(1)，所提取出的重力場水平梯度的最大值，即可確定重力場場源邊界及梯度變化情況，為重力場深部結(jié)構(gòu)的莫霍面反演奠定基礎(chǔ)。

1.2 重力場深度反演

由于地球的殼幔密度差為常數(shù)，因此在式(1)得到的重力場場源邊界及水平梯度基礎(chǔ)上，采用Parker-Oldenburg密度模型，在空間域中反演重力場深度，得到所有重力場場源參數(shù)，構(gòu)成重力場場源莫霍面結(jié)構(gòu)圖，分析重力場深部結(jié)構(gòu)。

但是在重力場場源反演過程中，考慮地球地殼存在密度分布不均勻、地幔物質(zhì)的質(zhì)量盈虧等問題，這會導(dǎo)致地球重力場因地殼局部場信息的干擾，產(chǎn)生重力場局部重力異常的現(xiàn)象，難以準(zhǔn)確反演重力場深部結(jié)構(gòu)[5]。所以，此次反演重力場深度將抑制地殼淺部的局部場信息，突出重力場深部信息。綜合上述分析結(jié)果，假設(shè)地球重力場中存在的重力異常為Δg；根據(jù)式(1)得到的重力場場源邊界，所分析得出的重力場界面起伏為h；重力場場源深度為Z0；重力場不同區(qū)域中，存在的地球殼幔密度差為Δρ，則反演出的重力場深度F[h(r)]為

(2)

式中：G表示萬有引力；r表示計算點的向徑；n表示空間域?qū)哟?；k表示反演系數(shù)[6]。此時，根據(jù)式(1)和式(2)計算的結(jié)果，以及分析出的重力場場源變化，采用Ueosoft軟件的Euler3D模塊，模擬重力場深部結(jié)構(gòu)。

1.3 形成重力場深部結(jié)構(gòu)圖

經(jīng)過對式(1)和式(2)的計算及分析后，將得到的重力場場源參數(shù)代入Ueosoft軟件的Euler3D模塊，形成重力場深部莫霍面結(jié)構(gòu)圖。為準(zhǔn)確分析重力場深部結(jié)構(gòu)，從而得到動力學(xué)特征，以某區(qū)域的重力場深部結(jié)構(gòu)為例，采用式(1)和式(2)進(jìn)行計算及分析，得到重力場場源參數(shù)，形成如圖1所示的重力場深部莫霍面結(jié)構(gòu)圖。

圖1 重力場深部莫霍面結(jié)構(gòu)圖

從圖1中可以看出，重力場深部結(jié)構(gòu)具有相對的均一性，且重力場場源強(qiáng)度值高，并呈現(xiàn)出自東向西逐漸降低的趨勢。由此可以推斷，重力場場源深度在26～33 km之間，處于穩(wěn)定狀態(tài)，且基本維持在這個區(qū)間之內(nèi)。但是，在此區(qū)間還出現(xiàn)了較大的正值變化和重力正值變化異常區(qū)，所以，在重力場下存在的基底物質(zhì)密度較高，且重力場場源深埋地下，不具備物質(zhì)上涌熔融的條件。而在重力深部莫霍面結(jié)構(gòu)圖中，還存在縱向上的延伸，橫向上的延展變化，形成由深至淺的物質(zhì)交錯斷層。由于重力場深部結(jié)構(gòu)中還存在深度淺且密度大的場源點，因此在形成重力場深部結(jié)構(gòu)的位置，一定會存在構(gòu)造轉(zhuǎn)換，并出現(xiàn)了斷裂的逆沖、山脈的隆升的地質(zhì)運動現(xiàn)象。

2 判斷動力學(xué)特征

根據(jù)上一章節(jié)中形成的重力場深部結(jié)構(gòu)圖像，及其圖像分析結(jié)果，采用動力學(xué)計算、判斷重力場深部結(jié)構(gòu)動力學(xué)特征。

2.1 大地水準(zhǔn)面異常反演動力學(xué)特征

從重力場深部結(jié)構(gòu)分析過程可以發(fā)現(xiàn)，重力場的深部結(jié)構(gòu)會存在一定的異常，從引言提到的眾多學(xué)者的研究成果中可以發(fā)現(xiàn)，重力場深部結(jié)構(gòu)可以影響地質(zhì)信息。所以，重力場深部結(jié)構(gòu)中存在的重力異常，會引起大地水準(zhǔn)面出現(xiàn)異常現(xiàn)象。

為此，模擬大地水準(zhǔn)面在重力場深部結(jié)構(gòu)重力異常影響下出現(xiàn)的異常變化。假設(shè)平均地殼厚度為H，地形高或海深為l，地幔密度為ρ1，地殼密度為ρ2，海水密度為ρ3，則大地水準(zhǔn)面異常N為

(3)

式中：g表示重力加速度。從圖1可以推斷，大地水準(zhǔn)面異常，與重力場深部結(jié)構(gòu)的密度、場源強(qiáng)度值大小等變化有關(guān)，其存在的構(gòu)造轉(zhuǎn)換過程會讓巖石出現(xiàn)向下?lián)锨默F(xiàn)象，產(chǎn)生均衡補(bǔ)償，此時，大地水準(zhǔn)面異常N為

(4)

此外，重力場場源邊界存在的變化也會影響大地水準(zhǔn)面，可能導(dǎo)致大地水準(zhǔn)面出現(xiàn)異常，則有：

(5)

式中：r′表示流動體元的向徑；dv表示積分體元。此時，即可根據(jù)式(2)～(4)得到的大地水準(zhǔn)面異常變化，判斷地球殼幔密度差，從而對大地水準(zhǔn)面的總場進(jìn)行地殼厚度反演。此時，將式(5)帶入，則有：

(6)

式中：t0表示彈性板撓曲界面深度；ZN(k)表示大地水準(zhǔn)面異常反射出的地殼厚度；D表示巖石層的彈性剛度。采用式(6)，即可反演地殼厚度，從而推斷出地殼厚度。

2.2 地球動力學(xué)特征

基于已有的重力場研究成果可以發(fā)現(xiàn)，重力場研究的主要目的是通過重力場推斷出地殼構(gòu)造活動變化，及其存在的動力學(xué)特征。在2.1節(jié)中，分析了大地水準(zhǔn)面異常及其反演后得到的地殼厚度等動力學(xué)特征。本小節(jié)將分析重力場等值線走向及其分布，表現(xiàn)出的地殼構(gòu)造活動的動力學(xué)特征。

從圖1所示的重力場深部構(gòu)造及其分析結(jié)果，以及大地水準(zhǔn)面異常反演動力學(xué)特征推斷出的地殼厚度，可以確定重力場深度的重力空間變化，會受到地質(zhì)塊體、局部構(gòu)造條件、介質(zhì)不均勻、地幔物質(zhì)等作用影響，且與重力異常的空間分布，具有一定的相關(guān)性，反之亦然。因此，基于此次針對重力場深部結(jié)構(gòu)的計算所形成的重力場深部結(jié)構(gòu)圖，其等值線走向和分布特征，反映了地殼構(gòu)造走向的變化規(guī)律，且存在梯度變化異常區(qū)域；梯度帶的范圍和量級在東西兩部也存在較大差異，并與構(gòu)造活動產(chǎn)生的斷裂帶存在密切關(guān)系。

3 結(jié)語

本文從異常角度來分析重力梯度變化帶和量級。由于此次形成的重力深部結(jié)構(gòu)圖中，重力場的重力強(qiáng)度存在較大差異，所以東西重力梯度變化帶和量級也存在較大變化，并存在過渡階段，影響著地殼厚度和梯度變化，出現(xiàn)梯度陡變現(xiàn)象，即較為強(qiáng)烈的物質(zhì)變遷和構(gòu)造變形差異。

在板塊方面，由于各大板塊之間存在相互控制，所產(chǎn)生的制約問題，因此在這過程中，所產(chǎn)生的異常地幔會向上隆升，影響到重力場深部結(jié)構(gòu)，減小重力場深部重力，從而改變地殼構(gòu)造活動，以及其產(chǎn)生的動力學(xué)特征。

此次分析重力場深部結(jié)構(gòu)與動力學(xué)特征的重點是重力場深部結(jié)構(gòu)。根據(jù)重力場深部結(jié)構(gòu)推斷其動力學(xué)特征，推斷由其引起的地質(zhì)變化，為地球地質(zhì)研究提供一定的依據(jù)。但是此次分析未涉及重力場變化與構(gòu)造應(yīng)力場變化存在的關(guān)系，及其對地質(zhì)構(gòu)造造成的影響。因此在今后的研究中，還需對這一領(lǐng)域進(jìn)行深入分析，為地質(zhì)研究提供更多的參考。