海域航空重力快速构建区域大地水准面

我国在海域开展了大规模的航空重力勘探,这些资料对构建高精度大地水准面具有重要价值.基于此,本文提出一种利用海域航空重力测量数据快速构建大地水准面的方法.该方法基于移去-恢复法思想,利用位场最小曲率方法对航空重力数据进行高精度向下延拓并获取相应的扰动位,实现航空重力测量快速构建海域大地水准面.与斯托克斯积分计算相比,采用了处理效率更高的频率域位场转换,解决了向下延拓及垂向积分时航空重力异常数据空白及扩边问题,具有较高的位场转换精度.本文应用EGM2008模拟航空重力数据进行模型验证,计算结果与其给出的水准面的精度相当;同时,也选取GRAV-D计划的航空重力数据进行实际验证,计算结果与xGEOID18B水准面模型精度基本一致.模型验证和实际应用验证了本方法的实用性.     

Numerical Simulation to Evaluate the Effects of Upward Lightning Discharges on Thunderstorm Electrical Parameters

A theoretical discussion of the discharge effects of upward lightning simulated with a fine-resolution 2D thunderstorm model is performed in this paper,and the results reveal that the estimates of the total induced charge on the upward lightning discharge channels range from 0.67 to 118.8 C,and the average value is 19.0 C,while the ratio of the induced charge on the leader channels to the total opposite-polarity charge in the discharge region ranges from 5.9%to 47.3%,with an average value of 14.7%.Moreover,the average value of the space electrostatic energy consumed by upward lightning is 1.06×10^9 J.The above values are lower than those related to intracloud lightning discharges.The density of the deposited opposite-polarity charge is comparable in magnitude to that of the preexisting charge in the discharge area,and the deposition of these opposite-polarity charges rapidly destroys the original space potential well in the discharge area and greatly reduces the space electric field strength.In addition,these opposite-polarity charges are redistributed with the development of thunderstorms.The space charge redistribution caused by lightning discharges partly accounts for the complexity of the charge structures in a thunderstorm,and the complexity gradually decreases with the charge neutralization process.     

基于泰勒级数迭代的重力归一化总梯度

针对重力归一化总梯度常规计算方法中向下延拓计算的不稳定性,提出了基于泰勒级数迭代法进行快速稳定向下延拓的重力归一化总梯度场( GH 场) 计算模式。水平圆柱体模型试验表明: 该方法的GH 场表现为单一高值特征,而且极大值位置与模型体质心位置相吻合。非均匀( 含油气) 球冠模型表明: 在油气含量较低时,Fourier 级数法和Fourier 变换法无法识别出“两高夹一低”的储油特征,而本文方法可以清晰地识别。这充分说明基于泰勒级数迭代的重力归一化总梯度计算方法不仅提高了计算结果的稳定性,而且还提高了异常的分辨率。     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

High primary productivity and f-ratio in summer in the Ulleung basin of the East/Japan Sea

To better understand the cause of high summer primary productivity in the Ulleung Basin located in the southwest part of the East/Japan Sea, the spatial dynamics of primary, new, and regenerated productivities (PP, NP, and RP) were examined along the path of the Tsushima Warm Current system in summer 2008. We compared hydrographic and chemical parameters in the Ulleung Basin with those of the Kuroshio Current in the Western Pacific Ocean and the East China Sea. In summer, integrated primary productivity (IPP, 0.37–0.96 g C m⁻² d⁻¹) and integrated new productivity (INP, 26–221 mg N m⁻² d⁻¹) within the euphotic zone in the Ulleung Basin were higher than those in the East China Sea and the Western Pacific Ocean (0.17–0.28 g C m⁻² d⁻¹, 2−5 mg N m⁻² d⁻¹, respectively). In contrast, there was no pronounced spatial variation in integrated regenerated productivity (IRP, 43–824 mg N m⁻² d⁻¹). Strong positive correlations between IPP and INP (also the f-ratio), and between nitrate uptake rate in the mixed layer and nitrate upward flux through the top of pycnocline in summer in the Ulleung Basin imply that the high IPP was mainly supported by supply of nitrate from the underlying water in the euphotic zone. Shallowing of the pycnocline depth as the current enters the East/Japan Sea facilitates nitrate supply from the nutrient-replete cold water immediately below the pycnocline through nitrate upward flux. A subsurface maximum in PP at or above the pycnocline and a high f-ratio further support the importance of this source of nitrate for maintaining the high summer PP in the Ulleung Basin. In comparison, the high PP layer was observed at the surface in the following fall and spring in the Ulleung Basin. Our results demonstrate the importance of hydrographic features in enhancing PP in this oligotrophic Tsushima Warm Current system.     

A geoid solution for airborne gravity data

Airborne gravity data is usually attached with satellite positioning of data points, which allow for the direct determination of the gravity disturbance at flight level. Assuming a suitable gridding of such data, Hotine’s modified integral formula can be combined with an Earth Gravity Model for the computation of the disturbing potential (T) at flight level. Based on T and the gravity disturbance data, we directly downward continue T to the geoid, and we present the final solution for the geoid height, including topographic corrections. It can be proved that the Taylor expansion of T converges if the flight level is at least twice the height of the topography, and the terrain potential will not contribute to the topographic correction. Hence, the simple topographic bias of the Bouguer shell yields the only topographic correction. Some numerical results demonstrate the technique used for downward continuation and topographic correction.     

大巴山弧形断裂（镇巴-高川段）的磁性特征及构造意义

秦巴山区的大巴山断裂（镇巴-高川段）及周边地区实施了高精度磁力勘探,经过向上解析延拓、垂向不同阶导数及欧拉反褶积等新方法处理及其分析,结果表明研究区磁场特征具有明显的规律性：横向上分别显示出沉积区、挤压构造变形区及侵入岩区的磁性特征,这是研究区多样化构造特征的具体反映;纵向上,浅部磁场呈现方向性磁场,深部呈现稳定基底的特征,这种磁性差异分界面构成了北大巴山地区逆冲推覆构造的滑脱面.而大巴山主断裂在不同区段表现出不同磁性结构的多样化特征,也反映了不同挤压应力作用下变质程度的差异.     

不同地质条件下隧道洞口仰坡地震破坏特性研究

针对洞口段均质围岩仰坡、含软弱夹层仰坡和桁架梁加固仰坡围岩3种工况,开展大型振动台模型试验,分析隧道洞口段仰坡模型土破坏形态。试验结果表明,均质边坡洞口段模型土在动力作用下坡肩土体先出现张拉裂缝,随着激振加速度增加,坡肩土体局部出现倾倒崩塌,最后沿坡面滑落堆积;含软弱夹层边坡在地震力作用下,仰坡坡脚部位土体挤压破碎,坡顶表面沿软弱夹层位置出现张拉裂缝,上覆土体沿软弱夹层滑动,最后土体大规模崩塌、滑落;洞口段加固边坡在动力作用下基本保持整体稳定,只部分梁格出现了局部的掉块,模型土顶部出现沿隧道轴向的细微裂缝。分析了隧道洞口段衬砌结构破坏形态,试验结果表明均质边坡洞口段AB两段衬砌模型裂缝分布较CD段复杂,认为洞口段地震影响深度为AB两段衬砌的长度,对应于实际工程中40 m;受软弱夹层影响,跨软弱夹层部衬砌模型裂缝形态较复杂;仰坡加固后,洞口段衬砌模型受力改善。研究结果可为山岭隧道洞口段边坡抗减震研究和设计提供参考。     

基于Broyden方法的不对称与对称周期解延拓算法与应用

考虑周期解的数值延拓问题并提出基于Broyden拟牛顿法来延拓周期解的一种有效算法,先后以布鲁塞尔振子、平面圆型限制性三体问题(Planar Circular Restricted Three-Body Problem, PCRTBP)的周期解为例进行了验证.这里的Broyden方法包含线性搜索、正交三角分解求线性方程组的步骤.对一般的周期解,周期性条件方程组中含有周期作为待延拓参数,可用周期来决定积分时长,将解代入周期性条件得到积分型的非线性方程组,利用Broyden方法迭代延拓直至初值收敛.根据两次垂直通过一个超平面的轨道是对称周期轨道的性质,可采用插值的方法求得再次抵达超平面的解分量,得到周期性条件方程组,再用Broyden方法求解.结合哈密顿系统的对称性和PCRTBP周期轨道的一些分类,对2/1、3/1的内共振周期解族进行了数值研究.最后,对算法和计算结果做了总结和讨论.