布格重力异常计算
野外重力测量的计算题
某工区,平均经纬度为E119.6°,N40.1°;σ=2.67g/cm3(中间层校正需要)。近区地改
0.008645mGal,远区地改0.062980mGal。
(1)早基点读数3699.019,时间为上午8:10;仪器高h=60mm。
固体潮校正值为-0.0696mGal,基点坐标为x:4436.440km,y:40465.710km,高程101m。
δg中=-0.0419*σ*Δh=-0.0419*2.67*(95.1970-101)=0.649199019
(3)高度校正
δg高=0.3086*Δh=0.3086*(95.1970-101)=-1.7908058
(4)正常场校正
δg正常场=-0.814*sin(2φ)*D=-0.814*sin(2*40.1°)*(XA-X基点)=-0.814*sin80.2*(4436.305-4436.440)=0.108286474
g''''A=g'''A+δg零A=3753.935478
(4)相对重力值
Δg''''=g''''A-g''''早=3753.935478-3752.816956=1.118521831
3、计算布格重力异常
(1)地形校正
δg地形=近改+远改=0.008645+0.062980=0.071625
(2)中间层校正
格值表
读数
累计重力差
格值因子
3500
3551.02
1.01422
3600
3652.44
1.01423
3700
3753.86
1.01425
3800
3855.29
1.01427
问题:
求格值转换、计算A点的相对重力值、求A的布格重力异常。
ห้องสมุดไป่ตู้解题思路:
1、格值转换:
g'=累计重力差+格值因子*(测点读数-格值表的标准读数)。
δg零晚=-(g'''晚-g'''早)*(t晚-t早)/(t晚-t早)
=-(3752.848080-3752.816956)*(14:23-8:10)/(14:23-8:10)=-0.0311232
-----------------------------------------------
g''''早=g'''早+δg零早=3752.816956
(1)仪器高校正,公式Δg早=0.3086*仪器高(单位m),即Δg早=3.086h
---------------------------------------------
例如Δg早=0.3086*60/1000=0.018516
Δg晚=Δg早=0.018516
ΔgA=0.3086*20/1000=0.006712
(5)布格重力异常
δg布格=Δg''''+δg地形+δg中+δg高+δg正常场=0.156826524
例如早基点:g'早=3652.44+1.01423*(3699.019-3600)=3752.86804
g'A=3753.86+1.01425*(3700.055-3700)=3753.915784
g'晚=3652.44+1.01423*(3698.859-3600)=3752.705764
2、计算A点的相对重力值
g''晚=g'晚+Δg晚=3752.705764+0.018516=3752.724280
--------------------------------------------------------------
(2)固体潮校正。g'''=g''+Δg固体潮
g'''早=g''早+Δg早固体潮=3752.886556-0.0696=3752.816956
-----------------------------------------
g''早=g'早+Δg早=3752.86804+0.018516=3752.886556
g''A=g'A+ΔgA=3753.915784+0.006712=3753.921956(结果为用Excel计算的,相加可能不是这个结果)
(2)A观测点读数3700.055,时间为上午10:23;仪器高h=20mm。固体潮校正值为0.02462mGal,A点坐标x:4436.305km,y:40465.735km,高程95.1970m。
(3)晚基点读数3698.859,时间为下午14:23,仪器高h=60mm,固体潮校正值为0.1238mGal。
g'''A=g''A+ΔgA固体潮=3753.921956+0.02462=3753.946576(结果为用Excel计算的,相加可能不是这个结果)
g'''晚=g''晚+Δg晚固体潮=3752.724280+0.1238=3752.848080
(3)零漂校正
-----------------------------------------------
δg零早=-(g'''晚-g'''早)*(t早-t早)/(t晚-t早)
=-(3752.848080-3752.816956)*(8:10-8:10)/(14:23-8:10)=0
δg零A=-(g'''晚-g'''早)*(tA-t早)/(t晚-t早)
=-(3752.848080-3752.816956)*(10:23-8:10)/(14:23-8:10)=-0.01109755