《重力正演计算》

合集下载

球体重力异常正演程序报告

球体重力异常正演程序报告球体重力异常正演是地球物理学中的一种重要方法，用于研究地下物质分布和地球内部结构。

本报告将重点介绍球体重力异常正演程序的原理、步骤和应用。

一、原理球体重力异常正演是基于牛顿引力定律和球体模型的数学计算方法。

根据牛顿引力定律，在球体表面上的任意一点，重力加速度可以表示为：g = G * (M / r^2)其中，g为重力加速度，G为引力常数，M为球体的质量，r为球心到该点的距离。

根据球体模型，球体的质量可以表示为：M = (4/3) * π * ρ * R^3其中，ρ为球体的密度，R为球体的半径。

将质量公式代入重力加速度公式，可得到球体表面上的重力加速度公式：g = (4/3) * G * π * ρ * R / r^2二、步骤球体重力异常正演程序的步骤如下：1. 确定观测点的位置和高度，以及球体模型的半径和密度。

2. 计算球体表面上的重力加速度，根据上述公式进行计算。

3. 根据观测点与球心的距离，计算球体表面上的重力加速度的投影值。

4. 重复步骤3，直到计算出所有观测点的重力加速度投影值。

5. 计算观测点的球体重力异常值，即观测点的重力加速度减去球体表面上的重力加速度投影值。

三、应用球体重力异常正演程序在地球物理勘探中有广泛的应用，主要包括以下几个方面：1. 地质勘探：通过球体重力异常正演，可以对地下的岩石密度分布进行推测，从而帮助地质勘探人员确定地质构造和找到潜在的矿产资源。

2. 油气勘探：油气藏通常与地下的密度异常有关，通过球体重力异常正演，可以对潜在的油气藏进行初步判断，指导油气勘探的方向和深度。

3. 地壳构造研究：地球内部的构造和演化与地下岩石的密度分布密切相关，通过球体重力异常正演，可以揭示地壳的变形和演化过程，为地壳构造研究提供重要的参考依据。

4. 火山和地震研究：火山和地震活动通常与地下的岩浆和断层有关，球体重力异常正演可以帮助科学家们理解火山和地震的发生机制，预测可能的灾害风险。

试验二规则二度体重力异常正演计算

握简单条件下规则二度体重力异常的正演计算的方法和形态特征，加深对课堂理论内容的理解。要求学生熟悉计算机常用编程语言，能够独立编制简单的计算程序，并能够使用常用绘图软件对所计算结果绘制剖面和平面图，为后续反演工作打基础。实验内容：编制简单条件下二度体模型的重力正演计算程序，给出不同形状地质体，同时给出不同产状和密度参数，计算出重力异常及其导数的正确结果，绘制相应异常图，观察曲线形态分析总结其规律。

重力异常正演资料

• 若以水平圆柱体的轴线作为Y轴，Z轴垂直向下，在轴线上取一
单位长度， dm d
• 若水平圆柱体有限长，则
密度均匀的水平圆柱体
l
g G
d
l [( x)2 ( y)2]3/2
(x
2Gh0l
h0 )(x2 l2
h02 )3/2
密度均匀的水平圆柱体
• 当 l 时，
g 2Gh0
密度均匀的球体Vg VFra bibliotekzG
v
( z)d dd [( x)2 ( y)2 ( z)2 ]3/2
密度均匀的球体
密度均匀的球体
Vg
GM
[x2
h0 y2
h02 ]3/2
密度均匀的球体
Vg
GMh0 ( x2 h02 )3/2
球体重力异常图
球体重力异常图
利用已知异常计算球体参数
重力异常正演
正问题与反问题
正问题也称为正演计算(Forward Calculation) 已知地质体的形状、产状和剩余密度等，通过理论计算来求得异常的分布和规律。
正问题与反问题
• 反问题也称为反演(Inversion) • 已知异常的分布特征和变化规律，求场源的赋存
状态（如产状、形状和剩余密度等）
正问题与反问题
正演计算是解反问题的基础，解反问题（反演）是目的
正问题与反问题
简单规则几何形体的异常
• 为了简化，假设地质形体孤立存在，密度均匀，地面水平，所取剖面为中心剖面。
• 规则形体：球体、水平圆柱体、垂直台阶、脉状体……
密度均匀的球体(点质量)
• 自然界中，一些近于等轴状的地质体，如矿巢、矿囊、岩株、穹窿构造等，都可以近似当作球体来计算它们的重力异常，特别当地质体的水平尺寸小于它的埋藏深度时，效果更好。

重力模型公式

重力模型公式重力模型公式，这可是物理学中的一个重要概念呢！咱先来说说啥是重力模型公式。

简单来讲，它就是用来描述两个物体之间引力大小的数学表达式。

就像地球吸引着我们，月亮围绕着地球转，这些现象背后都离不开重力模型公式的作用。

那重力模型公式到底长啥样呢？它是 F = G×(m1×m2) / r²。

这里的F 表示两个物体之间的引力，G 呢是引力常量，m1 和 m2 分别是两个物体的质量，r 是两个物体质心之间的距离。

我还记得有一次，我带着一群小朋友去科技馆玩。

走到一个关于重力的展示区，有个小朋友特别好奇地问我：“老师，为什么苹果会从树上掉下来呀？”我就笑着告诉他：“这就是重力在起作用呀，咱们可以用重力模型公式来计算苹果受到的地球引力有多大呢。

”然后我就给他简单讲了讲这个公式。

这公式在生活中的应用可多啦！比如说，咱们发射卫星的时候，就得好好算算卫星和地球之间的引力，保证卫星能在预定的轨道上运行。

要是算错了，那卫星可就不知道跑到哪里去咯！在天文学里，科学家们也是靠着这个公式去研究天体之间的相互作用。

比如说，通过观测行星的运动轨迹，再结合重力模型公式，就能推测出那些我们看不到的天体的质量和位置。

回到咱们的日常生活，你想想，为什么我们跳起来最后还是会落到地上？这就是因为地球对我们的引力。

要是没有这个引力，那可就好玩啦，咱们说不定就飘在空中，像气球一样到处乱飞！再比如，建筑工人在盖高楼的时候，也得考虑重力的影响。

如果不考虑清楚，房子可能就不稳，说不定哪天就歪了倒了。

还有啊，当我们扔东西的时候，东西飞行的轨迹也和重力有关。

扔出去的球会先往上飞一段，然后在重力的作用下开始下落。

总之，重力模型公式虽然看起来有点复杂，但它真的无处不在，影响着我们生活的方方面面。

通过对重力模型公式的学习和了解，我们能更好地理解这个世界的运行规律。

它就像是一把神奇的钥匙，能帮助我们打开探索宇宙奥秘的大门。

所以啊，同学们，可别小看这个公式，它里面蕴含着无穷的智慧和力量，等着我们去发现和运用！。

简述重力场的正反演问题

简述重力场的正反演问题
重力场的正反演问题涉及重力异常的正演和反演。

正演问题是给定地下某种地质体的形状、产状和剩余密度等，通过理论计算来求得它在地面上产生的异常大小、特征和变化规律，这是正向思维的问题。

反演问题则是依据已获得的异常特征、数值大小、分布情形等并结合物性资料来求解地下地质体的形状和空间位置等，这是逆向思维的问题。

重力正演是指根据地下地质体的形状、大小、密度等物理参数，利用重力场理论计算其在地球表面产生的重力异常。

重力反演则是根据实测的重力异常数据，结合物性资料，推断地下地质体的形状、大小、空间位置等信息。

重力正演是解决正问题的过程，它从地下地质体的物理参数出发，预测其在地球表面产生的重力异常。

重力反演则是解决反问题的过程，它从实测的重力异常数据出发，推断地下地质体的形状、大小、空间位置等信息。

重力场的正反演问题在地球物理学中具有重要的应用价值，例如在矿产资源勘探、地质构造研究、地下水资源调查等领域都有广泛的应用。

通过正反演问题的解决，可以更好地理解地球内部结构和动力学过程，为资源开发和环境保护提供科学依据。

重力正演、反演

2)当σ>o时，极大值一侧对应着上升盘，极小值一侧对应着下降盘，在极小值十分清晰且大干极大值的绝对值时，属正断层类型，反之则属逆断层类型。
二度铅垂柱体对于沿水平方向延伸较长而横截面近于矩形的矿脉，可以当成二度铅垂柱体来研究。在正演它的异常时，坐标系及有关参数的选取见图,用 (x+α)与(x一α)分别代替铅垂台阶各公式中的 x，并将结果相减，即获得这一形体的重力异常及各阶导数异常的公式:
当柱体的下底 H→+∞ 时，便可获得底部无限延伸的铅垂脉的相应公式Δg→∞
( x − a) 2 + h 2 V xz = Gσ ln ( x + a) 2 + h 2 h h 2ah V zz = 2Gσ (tg −1 − tg −1 ) = 2Gσtg −1 2 x−a x+a x + h2 − a2 ⎡ ⎤ x+a x−a 2a ( a 2 + h 2 − x 2 ) V zzz = 2Gσ ⎢ = 2Gσ 2 − 2 2 2 2 ⎥ ( x + a) + h ⎦ ( x + a 2 + h 2 ) 2 − 4a 2 x 2 ⎣ ( x + a) + h
GM GMD = 2 2 nD ( x1 / n + D 2 ) 3 / 2
x 1/n = ± D n 2 / 3 − 1
取n=2，得x1/2=0.766D(X正半轴)和x’1/2=－0.766 D (X负半轴)，说明异常半极值点的横坐标为球心深的0.766倍
4、当D不变，使M加大m倍时，异常也同样加大
[( x + a ) 2 + H 2 ][( x − a ) 2 + h 2 ] V xz = Gσ ln [( x + a ) 2 + h 2 ][( x − a ) 2 + H 2 ] H h H h ) − tg −1 − tg −1 + tg −1 V zz = 2Gσ (tg −1 x+a x+a x−a x−a ⎡ ⎤ x+a x+a x−a x−a − + − V zzz = 2Gσ ⎢ ⎥ 2 2 ( x + a) 2 + H 2 ( x − a) 2 + h 2 ( x − a) 2 + H 2 ⎦ ⎣ ( x + a) + h

两种新的长方体重力异常正演公式及其理论推导

+(
z) arctan ( z)R | 2 | 2 | 2 ( x)( y) 1 1 1
g(x, y, z) = G ||| ( x) ln{( y) + R} + ( y) ln{( x) + R}
( z) arctan ( x)( y) | 2 | 2 | 2 ( z)R 1 1 1
g(x, y, z) = G ||| ( x) ln{( y) + R} + ( y) ln{( x) + R}

两种新的长方体重力异常正演公式及其理论推导
骆遥 1, 2
1 中国科学院地质与地球物理研究所，北京（100029） 2 中国科学院研究生院，北京（100049） E-mail：geo@
摘要：在前人推导长方体重力场、磁场正演理论表达式工作的基础上，重新对长方体重力场正演理论表达式进行理论推导，提出了两种全新的长方体重力异常正演公式形式，并给出了全部的理论推导过程，对比模型正演计算结果表明，新导出长方体重力场正演理论表达式的正确。关键词：长方体，重力场，正演，积分中图分类号：P631
线数据单位为 g.u.
Fig2. The cubic model gravity contour map
5. 结论
综合前人对长方体重力场正演理论表达式的推导过程，并借鉴长方体磁场及其梯度场理论表达式的推导，推导出了两种新的长方体重力场正演理论表达式（11）式和（13）式，对比模型正演计算结果表明，新导出长方体重力场正演理论表达式是完全正确的。
似积分的推导[15~17]，对 2 的推导有：
2= 2( 1
z)2 ( {(
y) {( x) +
R}2

第四节地质体参数的计算重力勘探5-正反演

i ) ln
2 i 1 i2

2 i 1

2 i

(i1

i
)
tg
1
i i
tg1
i1 i1

（二）任意形状三度体
1、线元法
➢用一组垂直于y轴的平面
和一组垂直于X轴的平面分别切割地质体，则任意两个平面的交线包合在地质体之内的部分形成一个线元。
x 时， g Gf h
1
P(x,0)
●x
h2
h 1 △σ △h
2
△σ △h
主剖面异常曲线单调变化，断层正上方梯度最大；平面异常等值线呈条带状分布，与断层线平行。
在前述三个特征点上，异常值与埋深无关；异常形态与埋深有关，埋藏越浅，水平梯度越大。
等值线为一系列平行台阶走向的直线，在断面附近等值线最密，称为“重力梯级带”，且异常向台阶延伸方向单调增大。
第四节地质体参数的计算
正演与反演
正问题也称正演，是指给定地质异常体的形状、产状和剩余密度分布，通过计算得出重力异常的大小、特征和变化规律等。
反问题也称反演，是指根据重力异常的数值大小、变化规律等场的特征，结合已知的地质资料和地质体的物性参数，求解地质体的形状和空间位置等。
正问题从给定地球物理模型，通过数值计算或物理模拟，得出相应地球物理场的过程，目的是认识和掌握地球物理场的特征与场源之间的对应关系；
当α=90°（垂直断层）时，重力异常极大值与极小值绝对值相等，曲线以原点O为中心对称
当α<90°（正断层）时，下降盘一侧异常极小值明显
当α>90°（逆断层）时，上升盘一侧异常极大值明显

重力的概念及其计算方法

重力的概念及其计算方法重力是自然界中普遍存在的力之一，它对物体产生作用，使得物体始终朝向地球的中心运动。

本文将对重力的概念进行介绍，并详细阐述计算重力的方法。

一、重力的概念重力是地球或其他天体之间产生的相互吸引力，它是由于物体之间的质量而产生的。

根据普遍引力定律，重力的大小与两个物体的质量成正比，与其距离的平方成反比。

也就是说，物体的质量越大，重力越强；物体之间的距离越近，重力也会增加。

二、计算重力的方法计算重力的方法主要基于牛顿第二定律和万有引力定律。

1. 根据牛顿第二定律计算重力牛顿第二定律表明，物体所受到的力与其质量和加速度成正比。

根据这个定律，我们可以通过以下公式来计算重力：重力 = 质量 ×加速度在地球表面上，加速度可以被近似地视为9.8米每秒的平方。

因此，重力的计算公式可以简化为：重力 = 质量 × 9.8其中，质量的单位可以是千克或克，重力的单位是牛顿（N）。

例如，一个质量为50千克的物体在地球表面上的重力为50千克 × 9.8 = 490N。

2. 根据万有引力定律计算重力万有引力定律是描述两个物体之间引力的定律，它可以用来计算物体之间的重力。

根据万有引力定律，重力的大小与物体的质量和距离的平方成正比。

具体计算公式如下：重力 =（G ×质量1 ×质量2）/ 距离的平方其中，G是一个常量，被称为普适引力常量，其值约为6.67430 ×10^-11 N·(m/kg)^2。

质量1和质量2分别是相互作用的两个物体的质量，距离是它们之间的距离。

计算得到的重力单位也是牛顿（N）。

需要注意的是，计算重力时，质量的单位需要一致，距离的单位也需要一致。

在国际单位制（SI）中，质量的单位是千克，距离的单位是米。

三、总结重力是地球或其他天体之间相互吸引的力，它可以通过牛顿第二定律和万有引力定律来计算。

根据牛顿第二定律，重力等于物体的质量乘以加速度；根据万有引力定律，重力等于普适引力常量乘以两个物体质量的乘积，除以距离的平方。

地球重力场模型正演计算软件系统的设计

地球重力场模型正演计算软件系统的设计孙腾科;赵杏杏【期刊名称】《现代测绘》【年(卷),期】2014(000)003【摘要】This paper introduces the earth gravity model forward calculation (calculation of height anomaly ,gravity anomaly , vertical deviation etc .)-EGMCAL ,and presents the software framework of the overall design and function .The software is proved correct by using GFZ developed the related software ,and using EGMCAL to process the eastern China bridge vertical control network data ,the calculation results of three kinds of high precision model in this area .%本文介绍了自主开发的地球重力场模型正演（计算高程异常、重力异常、垂线偏差等）计算软件-EGMCAL ，给出了该软件的整体设计框架和功能。

利用GFZ（德国地学中心）研制的相关软件验证了该软件的正确性，并利用 EGMCAL对我国东部某大桥高程控制网数据进行了处理，得出三种精度最高的模型在此区域的计算结果。

【总页数】4页(P19-21,24)【作者】孙腾科;赵杏杏【作者单位】河海大学地球科学与工程学院，江苏南京210098;河海大学地球科学与工程学院，江苏南京210098【正文语种】中文【中图分类】P225【相关文献】1.计算机辅助堤防设计计算软件系统 [J], 徐泽平;徐建忠;梁建辉2.基于可扩展计算机网络设计软件系统的开发设计 [J], 洪辉3.可扩展计算机网络设计软件系统的开发设计分析 [J], 武装4.地球重力场参数计算程序的设计与应用 [J], 祖滢;田林亚5.可扩展计算机网络设计软件系统的开发设计分析 [J], 陈荣高因版权原因，仅展示原文概要，查看原文内容请购买。

重力的概念与计算

重力的概念与计算重力是自然界中的一种基本力量，它对我们的生活和宇宙的演化都起着重要的作用。

重力概念的理解和计算方法的掌握对于我们深入研究物理学和解决实际问题都至关重要。

本文将介绍重力的概念和重力计算的相关知识。

一、重力的概念重力是指物体之间的相互吸引力。

根据万有引力定律，任何两个物体之间都存在着引力，而且这个引力与它们的质量和它们之间的距离有关。

重力是由质量产生的，而不是物体的形状或大小。

在地球上，重力对我们的生活产生了明显的影响。

我们之所以能够站立在地面上，是因为地球对我们施加了一个向下的重力。

重力还决定了物体下落的速度和轨迹，是许多自然现象的根本原因，如水流的形成、天体的运动等。

二、重力的计算1. 重力的计算公式根据牛顿的万有引力定律，两个物体之间的引力可以通过下面的公式进行计算：F =G × (m1 × m2) / r^2其中，F表示两个物体之间的引力，G为万有引力常量，m1和m2分别代表两个物体的质量，r表示两个物体之间的距离。

2. 重力的单位根据国际单位制的规定，重力的单位是牛顿（N）。

万有引力常量G的单位是N·m^2/kg^2，质量的单位是千克（kg），距离的单位是米（m）。

3. 重力的计算实例以地球对物体的重力为例，地球的质量约为5.97 × 10^24 kg，半径约为6.37 × 10^6 m。

假设一个物体的质量为100 kg，距离地球表面为1000 m，则可以计算出地球对该物体的重力：F =G × (m1 × m2) / r^2≈ (6.67 × 10^-11 N·m^2/kg^2) × ((5.97 × 10^24 kg) × (100 kg)) / (6.37 × 10^6 m + 1000 m)^2通过计算，可以得出地球对该物体的重力约为980 N。

16重力勘探-重力异常正反演

] 3Gm
mh
(x h ) h
2
3 2 2
] G m 2h 2 x 2 (x h )
2 5 2 2
g zz
g x
g z
[Gm g z g z g zz z h
2h 2 x 2 ( x2 h ) h
5 2 2
] 3Gmh 2h 2 3 x 2 (x h )
x1
h1 h2 h 2 )
h x 1
4
1 1 4 2tg 2 h x1 1 1 4 2tg 2 h x1 1 1 4 tg 4 h x1 1 1 1 4 tg tg ( 1) 4 h h x1

n
2G n 1
（3）反演圆柱体半径
R
R
2
• 如果剩余密度已知，圆柱体半径：
（4）导数异常
g z max h2 x2 2G 2 (x h 2 ) 2
x 0
h 2 2G 2G 4 2 h h
△gz
h2 x2 g z 0 2G 2 h x0 2 2 (x h )
g ( x x) g ( x x)) 测量水平梯度原理 g x 2x
△g • 反演： • （1）等值线的球心（x0,y0,0）→异常体中心在地表的投影 A（x,0,0） • （2 ）反演深度
1 g 1 g max n n

△g
F r
h
g 1 G
n n
m h
h h2 h1
h h h h1 h2
2
O P(x,0,0) 1 2 Nhomakorabeal 2 l1 2 1

重力异常正演计算知识讲解

重力异常
• 1、重力异常的概念
• 地下物质密度分布不均匀引起重力随空间位置的变化。
• 在重力勘探中，将由于岩石、矿物分布不均匀所引起的重力变化，或地质体与围岩密度的差异引起的重力变化，称为重力异常。
重力异常
重力异常
• 在观测的重力值中，包含了重力正常值及重力异常值两部分。
• 用实测重力值减去该点的正常值，也能够得到重力异常。
VXZ4G S [( (x )2x )(( zz))2]2dd
密度均匀的水平圆柱体
• 对于剩余密度均匀的无限长水平圆柱体，可视为质量集中在轴线上的物质线
d d S （S 是水平圆柱体的横截面积）
S
g2G(x)h 20(h z0z)2
密度均匀的水平圆柱体
VXZ4G[( (x )2x)((hh00zz))2]2
VZZ2G[((h0x)z2)2(h(0zx))22]2
VZZZ4G 3([( x )2x()h20 (zh)0 (zh)02 ]3z)3
密度均匀的水平圆柱体
• 若以水平圆柱体的轴线作为Y轴，Z轴垂直向下，在轴线上取一
单位长度， dmd
• 若水平圆柱体有限长，则
重力场的等价性
• 重力场的等价性： • 地下不同深度、形状、密度的地质体
在地表面可引起同样的重力异常。 • 重力场的等价性给重力异常的解释带来一
定的困难。
简单规则几何形体的异常
• 为了简化，假设地质形体孤立存在，密度均匀，地面水平，所取剖面为中心剖面。
• 规则形体：球体、水平圆柱体、垂直台阶、脉状体……
V
(z)d d d
g z V zG v[(x)2 (y)2 (z)2]32
g V z V zG v[(x) 2( ( z)y d )2d (d z)2]32

重力数据处理解释方法

重力数据处理解释方法重力数据处理是指对地球或其他天体的重力测量数据进行处理和分析的方法。

通过重力数据处理，人们可以了解地球内部的物质分布和结构变化，研究地质构造、地壳运动和岩石的物理性质等。

下面将介绍几种常见的重力数据处理方法。

1.重力测量数据的收集与处理在进行重力数据处理之前，首先需要进行重力测量数据的收集。

常用的重力测量仪器有落体仪和重力仪等。

测量数据包括重力值、测量点的经纬度和高程等。

然后，对采集到的重力数据进行预处理，如去除仪器漂移、大气压力和海洋潮汐等影响因素，得到相对重力值。

2.重力异常的计算与分析重力异常是指实际测量值与参考重力值之间的差异，它反映了地下物质分布和地下结构的变化。

通过对重力数据进行异常的计算与分析，可以揭示地下构造和地质过程。

常用的重力异常计算方法有布格(Bouguer)异常和自由残差(Free-air)异常等。

布格异常是将测量值减去由海平面到测量点计算得到的理论重力值，同时考虑海平面以下的地壳质量；自由残差异常则是将布格异常再减去由海平面到一些参考高度计算得到的理论重力值。

3.重力数据的噪声处理重力数据中可能存在各种噪声，如仪器误差、大气压力和海洋潮汐等。

为了提高重力数据的质量和准确性，需要对噪声进行处理。

常用的噪声处理方法有滤波、平滑和插值等。

其中，滤波是通过将数据在频域进行变换，并去除高频成分来降低噪声影响；平滑则是通过对数据进行平均或加权平均来降低噪声的波动；插值是指通过已知数据点之间的关系来估算未知数据点的值。

4.重力数据的反演与解释通过对重力数据进行处理和分析，可以推断地下的物质分布和结构变化。

重力数据反演方法主要包括正演和反演两个过程。

正演是指根据已知的地下模型，通过数值计算得到理论重力数据；反演则是根据测量的重力数据，通过数值计算反推出地下的物质分布和结构。

常用的反演方法有二维反演、三维反演和重力异常分解等。

反演结果的解释需要结合地质、地球物理等其他数据，如地震资料和磁力资料等，以确定地下结构的精确性和可靠性。

正常重力值计算公式

正常重力值计算公式正常重力值的计算公式在地球物理学和大地测量学中具有重要的地位。

它可不是一个随随便便就能弄明白的简单式子，这里面藏着不少的学问呢！咱先来说说啥是正常重力值。

想象一下地球是一个特别大的圆球，而且这个球的物质分布特别均匀。

在这种理想情况下，每个地方的重力大小都有一个标准的值，这就是正常重力值啦。

正常重力值的计算公式那可是经过了好多科学家的努力才搞出来的。

比较常用的一个公式是 1901 年赫尔默特提出的。

这个公式长这样：γ = γ₀(1 + βsin²Φ - β/2sin²2Φ) 。

这里面的γ就是正常重力值，γ₀是赤道上的正常重力值，Φ是纬度，β是和地球有关的一个常数。

听起来是不是有点复杂？别担心，我给您举个例子您就明白啦。

有一次，我带着学生们去做实地测量。

我们来到了一个小山坡上，准备测量当地的重力值。

学生们都特别兴奋，一个个拿着仪器，跃跃欲试。

其中有个叫小明的孩子，特别积极，一直在问我各种问题。

我就跟他说：“小明啊，咱们要测的这个重力值，和咱们所处的位置，也就是纬度有关系。

你看这个正常重力值的公式，就能算出来咱们在这个地方理论上应该是多大的重力值。

”小明一脸迷茫地看着我，我笑了笑，接着说：“比如说，咱们现在在这个纬度是 30 度的地方，赤道上的正常重力值咱们假设是 9.78 米每秒平方，那个β常数假设是 0.0034，那咱们就能算出来这里的正常重力值啦。

”我带着他们一步步地把数据代入公式，算出来之后，再和实际测量的值进行对比。

学生们这才恍然大悟，原来这些公式不是在书本上干巴巴的数字和符号，而是能实实在在用来解决问题的工具。

在实际应用中，这个正常重力值的计算公式用处可大了。

比如说在地质勘探的时候，通过计算正常重力值和实际测量值的差异，就能发现地下是不是有异常的地质结构。

在航海和航空领域，也能帮助确定位置和导航。

总之，正常重力值的计算公式虽然看起来有点复杂，但只要咱们用心去理解，多结合实际的例子，就能发现它的妙处。

重力异常正演总结

一、正方体重力异常正演推导过程（蒋）直立长方体外任意一点的重力异常公式推导如下：p(x,y,z)：长方体外任一点ρ：长方体密度θ(ε,η,τ): 长方体内任一点坐标222111222321112221()(,,)[()()()]a b h a b h a b h a b h GM V G d d r V z g x y z G d d d Z x y z ρεηττρεητεητ==∂-∴∆==∂-+-+-⎰⎰⎰⎰⎰⎰设：r =；ln[()]u y r η=-+；则1.()xdu d y r r εεη-=-+要求原式：222111ln[()h b a a b hG y d ρηε--+⎰在此先求：ln[()y d ηε-+⎰22.()11.()()...()()().ln[()].[()].u d x u x x d y r rx x y r d y r r εεεεηεεηεη=-=----+-=--+--+⎰⎰⎰上式第二项：2().[()].x d y r r εεη--+⎰222222222().[()].[()()].()()().()()...()()[()()].x r y d x z rz y y z d x d d d x z r x z rεηεεττηητεεεεετετ---=-+-----=---+-+--+-⎰⎰⎰⎰⎰上式最后一项：222().().[()()].y z d x z r ητεετ---+-⎰2222222222222222222().().()2().()222().().[()()]..[()()].[()()]().().[()]..[().().()]().()..()().().y z x r x y ry z y z d r x z y z y z r x d r x y r z y z d r z x y ητεεηητητεετητητεεεητητετεη--------+-=-+--+-----=--+----=+---=-⎰⎰⎰212222().[.()]()..[().().()]().()rz x r z d x y r z y x τετεεητηε------+---⎰212().[.()]()..().()().()rz x r z r d d y x y x τετεηεηε⎡⎤----=⎢⎥----⎣⎦⎰2222.()().()1().().().()1().().arctan ().()r z x y z r z d y x z r z y x τεηττηεττηε---⎡⎤-=--⎢⎥--+⎣⎦⎡⎤-=--⎢⎥--⎣⎦⎰由上述推导过程，可以得出：ln[()y d ηε-+⎰().ln[()]().arctanx x y r z z εεηεττ-=--+-+--().().ln[()]().arctan ().()z r y x r z y x τηετηε⎡⎤-+--++-⎢⎥--⎣⎦ 222111(,,)ln[()]h b a a b hg x y z G y d ρηε∴∆=--+⎰222111().().ln[()]().arctan ().ln[()]().arctan ().()h b a ab hxz r G x y r z y x r z z y x ετρεηετηεττηε⎡⎤--=---+-+-+--++-⎢⎥---⎣⎦将上式完全展开后，().arctan xz z ετετ----可以消除。

球体重力异常正演程序报告

《应用地球物理学》课程作业基于MATLAB的球体重力正演程序实验报告1一程序简介本程序基于MATLAB软件的GUI模块编写，旨在实现球体重力正演结果的可视化分析。

MATLAB是一个高级的编程语言，其矩阵思想方便了地球物理的编程工作。

随着该语言和相应软件的发展，其内部也集成了许多模块，如该实验用到的GUI模块。

在该模块中，可以通过窗口、按键和赋值框等基本元素的组合，编写出可视化的应用程序，再配合MATLAB强大的作图功能，可以实现正演结果的展示与分析。

该程序应包含以下内容：1.可以自由输入参数，如球体半径，埋深和剩余密度。

2.可以计算出Δg、V ZZ、V XZ和V ZZZ这四种重力异常及其导数的对应值。

3.可以绘制剖面图及平面图两种图像。

二源程序由于GUI程序的头文件均大同小异，这里只列出赋值框及绘图按键的程序代码。

% --- Executes on selection change in popupmenu1.function popupmenu1_Callback(hObject, eventdata, handles)% hObject handle to popupmenu1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: contents = cellstr(get(hObject,'String')) returns popupmenu1contents as cell array% contents{get(hObject,'Value')} returns selected item from popupmenu1s=get(hObject,'value');handles.s = s;guidata(hObject, handles);function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text% str2double(get(hObject,'String')) returns contents of edit1 as a doubled = str2double(get(hObject,'string'));handles.d = d;guidata(hObject, handles);function edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text% str2double(get(hObject,'String')) returns contents of edit2 as a doubler = str2double(get(hObject,'string'));handles.r = r;guidata(hObject, handles);function edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text% str2double(get(hObject,'String')) returns contents of edit3 as a doublero = str2double(get(hObject,'string'));handles.ro = ro;guidata(hObject, handles);% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ro = handles.ro;r = handles.r;d = handles.d;x=-2*d:2*d;G=6.67e-11;pi=3.14159;switch handles.scase 2z=4*ro*pi*G*r^3*d./(3*(x.^2+d^2).^1.5);plot(x,z.*1e6);xlabel('X/m');ylabel('\Deltag/g.u.');case 3z=4*ro*pi*G*r^3*(2*d^2-x.^2)./(3*(x.^2+d^2).^2.5);plot(x,z.*1e9);xlabel('X/m');ylabel('Vzz/E');case 4z=-4*ro*pi*G*r^3*(d.*x)./(x.^2+d^2).^2.5;plot(x,z.*1e9);xlabel('X/m');ylabel('Vxz/E');case 5z=4*ro*pi*G*r^3*(2*d^2-3.*x.^2)./(x.^2+d^2).^3.5;plot(x,z.*1e9);xlabel('X/m');ylabel('Vzzz/nMKS');end% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ro = handles.ro;r = handles.r;d = handles.d;x=-2*d:2*d;y=-2*d:2*d;[X,Y]=meshgrid(x,y);G=6.67e-11;pi=3.14159;switch handles.scase 2c=4*ro*pi*G*r^3*d./(3*(X.^2+Y.^2+d^2).^1.5);contour(X,Y,c.*1e6,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 3c=4*ro*pi*G*r^3*(2*d^2-X.^2-Y.^2)./(3*(X.^2+Y.^2+d^2).^2.5);contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 4c=-4*ro*pi*G*r^3*(d.*X.*Y)./(X.^2+Y.^2+d^2).^2.5;contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 5c=4*ro*pi*G*r^3*(2*d^2-3.*X.^2-Y.^2)./(X.^2+Y.^2+d^2).^3.5;contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');end三程序运行结果假设埋深为100m，球体半径为10m，剩余密度为2kg/m3.Δg的剖面图和平面图如下：V zz的剖面图和平面图如下：V XZ的剖面图和平面图如下：V ZZZ的剖面图和平面图如下：。

重力怎么求物理公式

重力怎么求物理公式
牛顿在1687年出版的著作“原理”中，提出了重力场的定义。

法则性质：任何两个物体之间产生的引力是受他们之间距离的平方成
反比，且受物体质量的乘积成正比，从而得到重力定律的物理公式：F
＝Gm₁m₂/R²，其中G为万有引力常数，m1和m2为两物体的某一质量，R 为它们之间的距离，F为引力。

该公式表达的是任何两个物体之间产生的引力的大小，其中万有
引力常数G的值由实验测定而来，常数G是一种自由参数，可以被用
来计算物体之间的引力。

m1和m2为两个物体的质量，可用kg（公斤）单位表示。

R为两物体之间的距离，量纲为米（m）。

根据以上物理公式，可以知道物体之间的万有引力取决于它们的距离和质量，这就是
重力的定义。

例如，如果有一块重物重达20kg，另一块物体重达25kg，并且
它们之间的距离为5m，并假设G取为6.67×109Nm2/kg2，则它们之间
引力F可以用物理公式计算而得：F＝Gm₁m₂/R²＝
(6.67×109)(20)(25)/(5)²＝2.67×1010N。

总而言之，牛顿在17世纪提出的重力定律为我们提供了一个有
效的物理公式，即：F＝Gm₁m₂/R²，它告诉我们：物体之间的引力取决
于它们的距离和质量，它们的引力具有反比的性质，物体之间的距离
越远，引力越小。

《重力正演计算》

球体重力异常正演程序报告

试验二规则二度体重力异常正演计算

重力异常正演资料

重力模型公式

简述重力场的正反演问题

重力正演、反演

两种新的长方体重力异常正演公式及其理论推导

第四节 地质体参数的计算 重力勘探5-正反演

重力的概念及其计算方法

地球重力场模型正演计算软件系统的设计

重力的概念与计算

16重力勘探-重力异常正反演

重力异常正演计算知识讲解

重力数据处理解释方法

正常重力值计算公式

重力异常正演总结

球体重力异常正演程序报告

重力怎么求物理公式

第四节地质体参数的计算重力勘探5-正反演