地球内部合力强度方程建立
李务伦 李相通
摘要：地球是一个处于外部引力场，且存在自转的星体。地球内的任意一点无时不受自身引力场引力强度、外部引力场引力强度、自转强度以及浮力或沉力强度的作用，这些作用影响着地球构造，是地球构造运动的驱动动力；这些强度首次统一于一个方程之中，是构造运动具有周期性、方向性等的根源。
引言
地球内部存在自身引力场，不停地自转，同时还处在外部引力场之中，内部的热运动形成浮力或沉力，这些构成了地球的构造动力。对目前地球已发现的地球构造周期性、地球韵律[13]等，用这一构造动力进行解释，可能是不错的选择！
1.地球为球态的原因及内部合力强度方程
地球是球态、月亮是球态，所有星球都是球态。宇航员在太空中将相当量的液体从储液盒中挤出，如若在地面将摊成薄薄一层，而在太空，由开始的不规则，经不太长时间后，就变成球态。这是为什么？要回答这一个问题离不开万有引力定律。
万有引力定律表达式如下：
F=Gm1m2/r2………………（1）
G万有引力常数；m1、m2物体质量；r质量为m1、m2之间质心间距离；F引力，引力方向在m1、m2质心连线上
根据式（1），对于任意质量为m质点，质点在距离为r的一点p所产生的引力场强度[1]（或简称引力场，即一单位质点在p点所受引力）E为：
………………（2）
如图1所示，实心均质球外一点p的引力强度E可以通过以下步骤求出。设球半径为R且均质，密度为ρ，质量为M，OO1=y，图示球外P点距球心距离为r。圆平面O1垂直y轴，与平行于圆平面O1的另一圆平面相距dy，平行于Z轴的AO1B，O1A=O1B=k，相同质量质点A和B在P点的引力强度分别为dE＇和dE，方向如图示。对称的A、B两点在平行Z轴方向上合引力强度为零，所以半径为k的圆在P点平行平面OXZ的引力强度为零，圆平面O1内，A点或B点对P点平行Y轴引力强度为：，圆平面O1内半径为k的圆环在P点平行y轴的引力强度
由于球的对称性，所以球外任意方向一点P，根据Ey又可写为
（r＞R）………………（3）
高斯定理：引力场中引力通量用Φ表示，高斯定理可表述为：通过一闭合曲面的引力通量等于该曲面所包围的所有物质质量代数和的-4πG倍，即
。
图1，球内部任意一点的引力强度，既可积分方法求出，也可根据高斯定理得出。球内任意一点到球心距离为r的引力强度为：
（r＜R）………（4）
3）球体内的引力位
如图1，球外任一点P距球心距离为r的引力位V，可由式（3）积分得出，根据场论推知为：
………………（5）
半径为R球内部引力位，可由式（3、4）积分得出，根据场论推知为：
（r＜R）………（6）
4）球体内的压力
图2均匀球体，密度为ρ，半径为R，球内压力可根据式（4）计算，下面叙述这一过程：如图2所示半径为r的球上任取一微元ds，微元的高度为dr，因此体积为dv=dsdr，微体积的质量dm=dvρ，一单位质点为1=dvρ/dm，将其带入（4）式，于是有：
上式两边同乘dm于是得：
…（7）
（7）式右侧是两质点的相互吸引力，所以：
…（8）
图2
（8）式两边同除以就得微元对球面的压力dP，所以：
…（9）
根据力的可加性，对（9）式从r到R积分，于是得到：
…（10）
根据矢量线所满足的微分方程[2]15dx/Ax=dy/Ay=dz/Az，（3）和（4）引力线方程如下。
引力线方程为：
…（11）
从方程（11）得出，引力线均为经过球心的直线方程，所以所有引力线均垂直于半径为r的球面，而等引力强度面也为同样的球形，所以等引力强度面方程为：
x2+y2+z2=r2…（12）
方程（11）与（12）的关系可用图3表示。根据上面的分析，球内压力由（4）求出，所以压力方向线也为（11），等压力面方程为（12）；另引力位方程由（3）、和、(4)求出，所以等引力位面方程也为（12）。
根据上述叙述可得出单一物质形成的球内具有以下性质：①引力线为直线、压力方向均指向球心；②等引力强度面、等压力面、等引力位面均为球形；③引力线与等引力强度面、等压力面、等引力位面垂直；④球内部等引力强度面上，引力强度值处处相等；⑤球内等压力面上，压力值处处相等；⑥等引力位面上，引力位值处处相等。这些性质决定了太空中液体很快收缩为球形，推而广之，星体为球态的原因也源于此。
1.2均匀旋转球体内合力强度
如图3所示，球均质刚性，旋转角速度为ω。A为球内任意一点，距球心距离为r，该点的引力强度为式（4），A点处单位质点离心力为：
…（13）
单位质点离心力，也可以看做是离心力强度，这样一来与引力强度就统一起来了，进而可以将它们投影到图示的坐标轴上：
…（14）
…（15）
…（16）
…（17）
椭球为（a＞b）的极坐标方程为：
…………………（18）
对比（17）、（18）两式可以发现合力强度方程（17）也为椭球方程。如果图3球不具刚性，内部质点将依合力强度方程（17）运移，最后变为椭球形，因此地球为椭球形。根据方程（14、15、16）也可求出类似（11）强度线方程，以及由（17）求出等压力面和等位面方程，而这些方程只不过球形复杂，于要解决问题关系不大，因此对这些不做详细讨论，以后也如此。
为了使方程（17）更具有一般性和普遍性，将具有密度ρ=ρ(r)（r小于R）随半径变化，这时球内引力强度：
…（19）
…（20）
地球在赤道引力强度与离心力强度之比Er/Eω=4πGρ/3/ω2≈292。如此小的离心力强度，使得地球变为椭球形。
1.3外力场对旋转球体内合力强度的影响
地球处于太阳系中，太阳系位于银河系猎户臂上，月亮是离地球最近的星体。因此地球处在其它星球形成的合引力场中。月球在地球位置引力场强度有远、近和平均之分，其值分别如下：E远=2.98117×10-5N、E近=3.71384×10-5N、E平=3.31727×10-5N。如此小的引力强度引起潮起潮落及固体潮。太阳的质量占太阳系质量的98%，下面忽略其它行星作用，只算太阳在地球处引力强度。同样太阳与地球有远、近和平均之分，其值分别如下：E远=5.7349×10-3N、E近=6.1318×10-3N、E平=5.9282×10-3N。太阳系距银心2.6万光年，银河系质量为4.1771×1041kg（目前有的认为银河系质量为5.9673×1042kg），九成质量位于银盘上。很难推出从太阳到银河系中心所有物质形成的引力强度在太阳处合引力强度公式，但可根据“一均匀球体或均匀球层在其外一点所产生的引力强度等于将其全部质量集中于球心所产生的引力强度”和高斯定理可粗估其强度。根据高斯定理，2.6万光年以外的物质对太阳系的引力强度为零。2.6万光年以内物质保守估计不应少于银河系总质量的四分之一，即约1041kg（1042kg）。因此可粗估引力强度为E银=1.1×10-4N（或1.1×10-3N）。根据月球在地球处产生的引力强度两引力线与地球相切的最大夹角不足2o，可以认为这三者在地球处引力线是平行线，合引力强度引力线也平行，如图4所示，因此可设在地球处存在如下引力强度方程：
……（21）
式中：t时间，r位于球内距球心距离，σ与Z轴夹角，τ在XY平面内与X轴夹角
在图4的A点，球自身、自转、Ew合力分量为：
…（22）
…（23）
…（24）
合强度为：
…（25）
…（26）
月球在地球的引力强度量级前面已叙述，如此小的量级则引起地球引力潮，太阳和银河系的引力量级也会在构造中留有痕迹。
1.4合力强度方程成立证据
太阳系绕银心一周约2.5Ga，8次穿越银道面，穿越银道面时间间隔是30～35Ma左右。在地球穿越银道面时，陨击撞击地球，与地球穿越银道面时间相合[3]；还有学者对地质事件与穿越银道面做了研究，研究表明地球各圈层主要地质事件的旋回周期如表1和图5所示[4]35-36。从表1和图5不难发现银河系对地球的引力强度，影响了地球的方方面面。地球的韵律分级如表2所示，从表2更能看出地外引力强度对地球的作用；地球存在经向和纬向构造带，以及54种朝东定向的构造现象[13]42-43，它们是地球自转导致的，因此可以断定方程（25）或（26）是成立的。
表1 地球各圈层主要地质
事件的旋回
地质事件
旋回周期/Ma
二级海平面旋回
35±3
海底扩张不整合
34±2
火山活动增强期
34±2
构造运动加强期
33±3
古气候旋回
33±3
碳酸盐侵入
34±1
生物集群绝灭
30±4
天体撞击事件
32±2
暗色溢流玄武岩
32±1
金伯利岩侵入
35±1
地磁极倒转
30±2
图5 地史上35Ma左右的重要地质事件集中期与穿越银道面时间对比
表2 地球的韵律分级
分段
时间尺度/a
地、气、水、天的韵律性变动代表事件
天体的代
表性效应
科学意义
分级
地壳
大气
表水
长韵律
109～106
2～1Ga
构造宙-壳、幔、核岩浆等的演化阶段
气演变宙-成分演变(O2、CO)、巨冰
海演变宙
银河效应：太阳-银心距、旋臂结构、麦哲伦星云
地球和太阳系的演化、起源
1～0.25Ga
构造期-造山期、成矿期、岩浆期
气演变期-大冰期
海演变期
成矿阶段、期的划分、成矿预测
250～2Ma
构造幕-褶皱幕、洋底变动幕
气演变幕
海演变幕-海水进退、海积大旋回
地球均态动力学
中韵律
106～103
2～0.1Ma
地文期-夷平面、洋脊张裂幕(磁极性分期)
气文期冰期、雨期、森林-草原交替期
湖文期-湖泊盛衰分期
米兰科维奇效应：地球轨道偏心率、黄赤交角、岁差
基础自然环境学
100～10ka
地文幕-大阶地、洋脊张裂幕(磁极性韵律)
气文慕-冰碛线进退
湖文幕-湖积旋回
地文期与近代生物演变
10～1ka
地文阶-活断层错动分
气文阶
湖文阶-湖积韵律、洪积扇叠置
地球次匀态动力学
段(河系错位、断层陡坎变动)
短韵律
103～1
1～0.1ka
地震世-古地震、古火山、古土壤、古文化层
气候世-小冰期、结冰温度变化
水面世-海蚀阶地
行星会聚效应：太阳黑子视效应，地球自转效应
自然灾害的长期预报
100～10a
地震期-地震盛衰分期、火山活动分期
气候期-温度世纪周、树轮
水面期-贝壳堤
自然环境与社会
10～1a
地震幕-地震活动幕
气候幕-环流形式变换、台风多少
水面幕-潮线
地球暫态动力学
微韵律
1～10-3
1～0.5a
地形年变-地面升降、伸缩、倾斜
天气年变-气温、气压、降水
水面年变-年升降
太阳辐射效应，月、日潮汐效应
自然灾害的短期预报
30～15d
地形月变-月潮
天气月变-天气过程的更换
水面月变-月潮
多种自然参数的相互关系
1～0.5d
地形日变-日潮
天气日变-气温
水面日变-日期
地球实态动力学
2物体的浮沉
图6（a）中，左侧长方体。在外力F作用下侵于液体中，其上作用力F=V(ρ2-ρ1)E,式中E为引力场强度，V为长方体体积，当图示該力不存在时，此时的-F为长方体向上运动力即浮力，而长方体质量为m=Vρ1，那么单位质量的浮力为：
……（27）
由于E为引力强度，（ρ2-ρ1）/ρ1为常数，所以E浮可以称为浮力强度，其方向与引力强度方向相反；使得引力场在长方体所处位置引力强度减少了E浮，因而物体上浮。
当长方体的力F去掉，长方体将如图6（b）示的上浮。长方体的向上运动，E浮的存在引起长方体上部压力增加量为△P＞0和质量相对增多，图中等压力线L1因此有水平如钟形上展，变为图示的L1。长方体的上升，长方体底面对接触底面的液体，因E浮的存在引起长方体底部压力增加量为△P＜0和物质量相对减少，使得图示的等压力线L2有水平如钟形下延为L2。长方体上部物质将如图示的方向运动，而长方体周边物质在黏滞力的作用下如图示向上运移，从而形成图示的流体旋转。
对于图6（a）右侧长方体，当上提力F去掉后将下沉，运动情况与上浮相反，如图6（c）所示，沉力强度为：
…（28）
图7 重力剖面图
在图7[1]中自由空气异常极小，表明洋中脊是均衡，但中部对应洋中脊形态，由洋中脊最高处向两侧，自由空气异常还是由大向小变化的趋势，说明物质有富余；布格异常是负的，表明洋脊下面质量短缺，这成了一对矛盾。大洋中脊为海底扩张的地方，热的物质上涌而使海底扩张形成，因而具有物体上浮的特征，反映在等压力面上就具有6（b）的形态，这种形态正是洋脊存在质量即富余又短缺的情况，这样质量富余和短缺的矛盾得以解决。
3热结构中合力强度方程建立
地球于星云中吸积而成，聚首的星云物质存在力、化学势、物质展布的三种不平衡，物质并非严格依照方程（11），（12）排列，也非时下的地球具有圈层。以往的研究认为，重组和运移是在热结构中进行的[5]、[6]、[7]、[8]，热结构如8所示。物质的重组和运移一分为三，小密度物分异后到热结构的图8的左右上角，大密度物聚于加热结构的底部。
星云物质聚集后，因引力位能而生热，热使得物质重组迁移，并使粒子间距离增大，从而使得物质上浮，同时使得其内部大比重粒子下沉聚集，这其中就包括放射性重元素U、Tu等，这样位能转为热能和衰变热（链式核裂变热）[9]形成上部的热源供给地，上部物质如图示持续上升，到最高后返回到出发地，热结构形成。对上部供热，是热源；热源促进热结构中粒子上升，在热结构底部形成负压，既然是源，就符合拉普拉斯方程∆u=0，因此就有了图8下部热源、负压源。
中间上升部位是因与周边存在密度差，也就是存在式（27）的浮力强度。设该浮力强度位于在图4中的A点，浮力强度与Z轴夹角为λ，在XY平面内与X轴夹角为μ，因此浮力强度在各轴的分量如下：
…（29）
…（30）
…（31）
将上述三式与式（22、23、24）对应相加，其A点合强度为：
…（32）
…（33）
负压是由上浮形成，因此负压也存在负压强度，它等于负E浮。在图8物质上升到最高点，开始向四周运动并降温，在图8的沟处回返，形成俯冲。因此在俯冲带内存在下沉和负压，所以在图4中A点，沉力强度与Z轴夹角为θ，在XY平面内与X轴夹角为δ，因此沉力强度在各轴的分量如下：
…（34）
…（35）
…（36）
设在图中A点，负力强度与Z轴夹角为γ，在XY平面内与X轴夹角为ψ，因此负力强度在各轴的分量如下：
…（37）
…（38）
…（39）
上述从（34）到（39）与（22、23、24）对应相加，A点强度为：
…（40）
…（41）
方程（32、33、40、41）均为椭球方程，是多因素方程。
4小密度物与热结构间的平衡分析
图9为图8左上的一部分。小密度右下与热结构接触处，对小密度任意一质点，存在图9示的受力：热结构对小密度的支撑力N支，热结构对小密度物摩擦力F摩，摩擦力使得左侧小密度物产生图9示的反抗力F抗，及上覆小密度物对该质点重力F引。在图示的坐标系中F抗与X，F引与Y夹角均为α。在X、Y轴的合力分别为：
FX=F摩+F引sinα-F抗cosα…(42)
FY=N支-F引cosα-F抗sinα…(43)
当FX=0，沟的位置稳定；FX＞0，沟的位置向图示左下移，反映热结构扩大，小密度物漂移或消减；FX＜0，沟的位置向图示右上移，反映热结构减弱，沟后退。正是（32、33、40、41、42、43）这些方程的存在，使得地球构造具有多样性。
5热结构存在验证
图10
帕皮提是一火山岛，离汤加弧不远，可以先假定其下存在图8式的点源热结构。用测线验证是否存在图8的海底地形。图10是谷歌地球在帕皮提和汤加弧间的截图，图中Ⅰ为帕皮提到海沟一条测线，在测线上等距离连续取得海底高程，并对应的做出海底测线地形图。在海底测线地形图上，除A处存在凸起外，其它处虽有波动，但不及A处，对这一测点线，做一手工拟合，图中红色线既是。从手工拟合线上可以看出：帕皮提到海沟的海底地形，基本上与图8 S到G相同。测线Ⅱ测线Ⅲ，基本和测线Ⅰ相同，可见热结构的存在。
当然热结构的存在，不止上面的证据，如海沟吞噬海水[14]太平洋下、大西洋及非洲下存在两个“大型低速切变区”，以及图7的重力图的解释等，均可佐证热结构的存在。在文献[8]中还对转换断层，洋脊链，三联点做了理论上的分析，有对热结构形成再了解的同行可参考文献[8]。
6结语
通过对引力场的分析，自转等的分析，在地球内部建立了由各种力强度的合力强度方程。从合力强度方程看，各家学者在同样的地质体中，提出不同的动力学理论是有道理的，但如何用合力方程解决地质问题，以后还将进行探讨！然而，随着洋底“消失的古陆”的被发现[10]；及“现今地球上裂谷的伸展并不一样，已经存在的洋中脊分布方向也不统一，为什么会形成这样分布的大洋中脊系统？…”问题提出[11]；“板块自驱动模式”[12]等的出现。也许上述系统性方程的给出，对问题的最后解决,若错误可提供可避免的覆辙，若正确或许是一可借鉴的思路，因为本文给出的合力强度适用于地球内部任何一点。也许这种探讨不对，望批评指正！
参考文献
[1]傅承义 陈运泰 祁贵仲 地球物理学基础[M] 科学出版社 1985年
[2]谢树艺 矢量分析与场论[M] 人民教育出版社 1978年12月
[3]万天丰 尹延鸿 全球岩石圈板块为什么会运移？[J] 自然杂志 2019第41期
[4]巫建华 刘帅 大地构造学概论与中国大地构造学纲要[M] 地质出版社 2008年10月第一版
[5]李务伦 地球成长及构造形成的热力学分析[J] 吉林省教育学院学报 2018年第7期 179-183
[6]李务伦 引力与热力的博弈--大地构造的二力作用[J] 山西科技 2019年第4期 23-29
[7]李务伦 大地构造运动的动力系统架构及应用[J] 山西科技 2019年第6期 37-44
[8]李务伦 李相通 引力热力在地球演化中的作用探析[J] 山西科技 2020年第4期 48-59
[9]马学昌 地壳运动驱动力的探讨--核能与地球演化[J] 地质学报 2016年1月 24-35
[10]任继舜 徐芹芹等 寻找消失的大陆[J] 地质论评 2015年9月 969-989
[11]赵文津 大陆漂移，板块构造，地质力学[J] 地球学报 2009年12月 717-713
[12]梁光河 大陆漂移的源动力[J] 地球科学前沿 2013年3月 86-96
[13]巫建华 刘帅 大地构造学概论与中国大地构造学纲要[M] 地质出版社 2008年10月第一版
[14]Chen cai,Douglas A.Wiens,Weisen Shen & Melody Eimer Water input into the Mariana subduction zone estimated from ocean-bottom seismic data Nnture Published 14 November 2018 563,389-392(2018)
