黃勁松 傅 瑜
(中國運載火箭技術研究院,北京 100076)
g
與一些因素的集合有關,這些因素包括所研究點的地理位置和高度、地球、太陽和月球的重力、地球圍繞自身軸和沿繞太陽軌道的轉(zhuǎn)動、地球軸相對黃道面的傾角、一年和一晝夜的時間等,即它是對不同位置和時間取不同值的變量。準確地求出以m/s、cm/s、伽(1伽=1cm/s)、毫伽(1毫伽=0.001cm/s)計量的自由落體絕對加速度,對宇航學、導航、大地測量、估計全局地球物理氣象過程以及軍事應用等都具有非常重要的意義。g
=981274±3毫伽的形式。γ
值的計算進行較多研究,他們在理論探索和在波茨坦坐標系進行大量實際測量試驗結(jié)果分析的基礎上,提出了一系列在精確球體條件(H
=0)下確定γ
值的解析關系式。1930年于斯德哥爾摩市舉行的大地測量會議上認定下面的Cassinis公式為國際公式(γ
值的計量單位是毫伽)γ
=978049.
0·(1+0.
0052884sinφ
-0.
0000059·sin2φ
)(1)
但之后在詳細完成不同地理點的多次重力測量與衛(wèi)星觀測數(shù)據(jù)分析的基礎上確認,波茨坦系把加速度值提高了大約14毫伽。因此在1971年于莫斯科市舉行的國際大地與地球物理測量會議上曾確立了新的確定標準的重力加速度公式,該公式考慮了觀測點在海平面以上高度H
的影響,其表達形式如下:γ
=978031.
85·(1+0.
005302·sinφ
-0.
0000059·sin2φ
)-0.
3086·H
(2)
公式(2)給出了考慮物體緯度和高度位置時重力加速度的標準值(近似值),為了更精確地確定絕對值g
,必須附加地考慮由各種因素作用造成的異常值(簡化),為此本文介紹一種考慮太陽和月球影響的方法重力加速度值計算方法。Δg
與所研究時刻月球在軌道上的位置有關,特別是該軌道被地球-月球地心軸與地球表面交點的坐標(φ
,λ
)確定。于是,考慮上述因素的重力加速度g
的廣義量由公式(3)確定:g
=γ
+Δg
+Δg
(3)
式中:γ
——考慮到所研究點的緯度和高度位置,按公式(2)確定的自由落體標準加速度;Δg
,Δg
——分別為由太陽和月球引力以及地球圍繞太陽運動形成的引力加速度值的異常。(4)
式中:M
——太陽的質(zhì)量1.
989·10kg
;G
——重力常數(shù)6.
6720·10Nm
/
(kg
·s
);α
——地球軸相對黃道平面的傾角23.
5°;-90°≤φ
≤90°——所研究點的地理緯度。它與地心緯度之間的差別很小(<0.
1°),對北半球、南半球和赤道的取值范圍分別是φ
>0,φ
<0,φ
=0;R
——從所研究的海平面上高度H
處的點到地球中心的距離,按公式(5)計算R
=R
+H
(5)
式中:R
——球形地球的矢徑;R
——地球的赤道半徑,R
——地球的極半徑,其值分別為R
=6378.164·10m,R
=(1-ε
)R
=6356.799·10m;ε
——地球扁率,0.
00335;按公式(6)計算(6)
式中:η
——考慮一年時間的系數(shù);ξ
——考慮一晝夜時間的系數(shù),S
——數(shù)值上等于天數(shù)的以度為單位的角度:對于第一個半年從所考慮的日期到6月22日,而對于第2個半年從6月22日到所考慮的日期(太陽系半年);t
——當?shù)匾粫円箷r間內(nèi)的一個數(shù)(以小時記);按公式(7)計算(7)
式中:V
——地球在黃道平面內(nèi)沿黃道軌道運動的線速度當前值,V
——地球繞太陽軌道運動在遠日點的線速度,29.27×10(m/s),V
——地球繞太陽軌道運動在近日點的線速度,30.27×10(m/s)。按公式(8)計算(8)
式中:L
——所研究的日期從地球到太陽的日心距,L
,L
分別表示在遠日點(6月22日)和近日點(12月22日)上從地球到太陽的距離L
=152098.233·10m,L
=147098.291·10m,e
——地球軌道的偏心距,0.0167112。φ
,λ
),按公式(9)計算月球重力F
對重力加速度的影響Δg
=f
(F
)(10毫伽)(9)
對地球朝向月球方向該公式取負值“-”,而對地球的背對月球方向取正值“+”。
式中:M
=7.35·10kg為月球的質(zhì)量;L
——從月球到所研究的地球表面上坐標為(φ
,λ
)的點的距離,其中-180°≤λ
≤180°為所研究點的地理經(jīng)度(λ
>0,λ
<0,λ
=0相應地對應于東經(jīng)、西經(jīng)和零子午線),按公式(10)計算L
=L
-R
·cosσ
(10)
式中:L
——地球到月球的當前距離,對于月球的橢圓軌道按公式(11)計算(11)
式中:e
——月球繞地球軌道的偏心率,0.
055;a
=384.40·10m,月球軌道的長半軸(沿拱點線);θ
表示月球軌道平面內(nèi)地球-月球地心軸相對拱點線(相對近地點)的偏角(頂點在地球中心),根據(jù)計算日期月球的位置求出;σ
——頂點在地球中心的地球-月球地心軸坐標(φ
,λ
)和所研究點(φ
,λ
)的地理位置之間的夾角,按公式(12)計算(12)
式中:Δφ
=|φ
-φ
|;Δλ
=|λ
-λ
|,當Δλ
>180°時應該按照Δλ
=360°-|λ
-λ
|計算。角σ
的變化范圍是:對于地球朝向月球的那一面0≤σ
≤90°;對于地球背對月球的那一面90°≤σ
≤180°。根據(jù)專門計算所研究日期上月球在繞地球軌道上的位置計算出點(φ
,λ
)的坐標值和角度θ
,在這種情況下緯度坐標的變換范圍-18.
36°≤φ
≤18.
36°,當月球沿拱點線位于軌道上時,近地點和遠地點上達到極值(最小和最大值)。零值φ
=0對應于月球處于節(jié)線上(在升交點和降交點附近),這時θ
=90°。圖1 太陽引力及離心力相關加速度
圖2 太陽引力及離心力相關加速度和ΔgS
在早晨和晚上(t
=6,18),因為作用力朝向地球表面的切線方向,所以太陽對重力加速度不產(chǎn)生影響。除此之外,對地球繞自身軸和沿繞太陽軌道運動的分析得出,對于躺臥狀態(tài)的人體組織來說,最有利的狀態(tài)是沿子午線躺臥。當身體轉(zhuǎn)變到頭朝西或者是頭朝東時,夜間和早晨單方向的切向和切向-離心加速度使血液向頭部或腳部流動。
對于近日點夜間與白天(t
=0;12,S
=180°)地球上不同點處考慮高度H
的按公式(2)計算出的最大的標準加速度γ
和按公式(3)~(5)計算出的廣義加速度g
值的比較結(jié)果在表1中給出,括號中給出標準加速度γ
和廣義加速度g
值之間的差值。表1 近日點時各地晝夜的加速度
從表1中可以看出,與一年中和一晝夜內(nèi)的時間有關,廣義加速度g
取與標準加速度γ
不一樣的變化的值,相應地達到30.07和-28.00毫伽。在近日點(P)北極(t
=0;12)加速度g
達到最大值983229.81毫伽,在欽博拉索火山頂峰(t
=12)達到最小值976073.25毫伽和最小值。在兩極的加速度的大小有不同的值——北極的大,南極的小。對于地球上不同點,月球引力F
對重力加速度的影響Δg
=f
(F
),在2019年12月26日(S
=179°,t
=12)新加坡市全日食期間(H
=15m
,φ
=1.
35°,λ
=103.
85°),在地球-月球地心線與節(jié)點線上地球表面的交點處,即θ
=90°廣義加速度g
的計算結(jié)果在表2中給出。從表2中可以看出,月球引力F
對重力加速度的影響取變符號的值——在地球朝向月球的這面取負值,在地球背對月球的面取正值。對于所研究的情況,加速度Δg
=f
(F
)的最大值(絕對的)等于-3.29毫伽,在地球-月球地心線與節(jié)點線上地球表面的交點處。表2 新加坡全日食期間各地的加速度
在自由落體標準加速度國際公式的基礎上,得到了與地球、太陽和月球引力、海平面以上的高度、一年和一晝夜間的時間有關的地球表面任意一點的廣義加速度的計算公式,并得出在兩極的加速度的大小有不同的值——北極大南極小。
地球在近日點時,加速度在北極(子夜、中午)和欽博拉索火山頂峰(子夜)處,分別達到最大值983229.81毫伽和最小值976073.25毫伽。在一晝夜期間由重力和離心力造成的變號加速度分量沒有被完全“平衡”,例如,對欽博拉索火山在近日點子夜分別是-568.57毫伽和+540.56毫伽。
在近日點太陽對地球上重力加速度施加的影響最大,一晝夜間加速度波動幅值達55.34毫伽(對欽博拉索火山夜間加速度增加達27.34毫伽,而白天減少達-28.00)。在早晨和晚間太陽對加速度的影響不大,因為重力沿地球表面的切線方向。
來自月球引力的加速度對地球朝向月球的表面取負值而對背對月球的那一面取正值,并且其絕對值相對較小(達3~4毫伽)。
由太陽和月球影響帶來的一晝夜和一年間加速度的不對稱波動對地球運動的集成作用最終被地球繞自身軸旋轉(zhuǎn)和沿軌道運動的特點所“吸收”。