基于神经网络方法获得最优化月球内部结构模型

由于观测手段有限,目前对月球内部结构的认识还存在很大的不确定性,至今仍没有一个被广泛认可的内部结构模型,且现有对月球内部结构模型的研究几乎很少关注观测值对观测精度的影响.本研究采用混合密度神经网络方法得到了月球内部结构模型的后验概率密度分布,获得了平均月球内部结构模型(Mean模型)、最大后验概率对应的月球内部结构模型(MAP模型)以及满足1-σ准则的月球内部结构模型(1-σ模型),其中MAP模型即为本文给出的最优化月球内部结构模型.此外,研究结果表明月球低速区S波波速低于月幔S波波速,因此本文结果支持月幔底部存在一个低速区的观点.不同观测值观测精度对模型影响的研究结果表明,勒夫数k₂存在一个约为0.0220的下边界,且其观测精度对月球内部结构模型的影响显著大于平均密度和平均转动惯量.

     

Steady- and transient-state inversion in hydrogeology by successive flux estimation

A calibration method to solve the groundwater inverse problem under steady- and transient-state conditions is presented. The method compares kriged and numerical head field gradients to modify hydraulic conductivity without the use of non-linear optimization techniques. The process is repeated iteratively until a close match with piezometric data is reached. The approach includes a damping factor to avoid divergence and oscillation of the solution in areas of low hydraulic gradient and a weighting factor to account for temporal head variation in transient simulations. The efficiency of the method in terms of computing time and calibration results is demonstrated with a synthetic field. It is shown that the proposed method provides parameter fields that reproduce both hydraulic conductivity and piezometric data in few forward model solutions. Stochastic numerical experiments are conducted to evaluate the sensitivity of the method to the damping function and to the head field estimation errors.     

反演地层电阻率的单调同伦法

本文从电法测井的数学模型出发,给出了一种反演地层参数的单调同伦法,把非线性方程组的求解转化成常微分方程初值问题的数值求解,从而给出一种大范围收敛的反演方法.最后的数值结果表明,本文给出的算法是十分有效的.     

柔度法纤维模型在方钢管混凝土柱滞回仿真中的应用

本文编制了基于二阶柔度法纤维模型梁柱单元的矩形钢管混凝土框架结构非线性分析程序,采用虚拟阶步法和改进的滞回路径追踪方法消除了其中截面切线刚度矩阵奇异时无法求逆的隐患,提高了算法的稳定性.运用此程序对矩形钢管混凝土柱进行滞回仿真,并与试验结果对比,结果表明:柔度法纤维模型梁柱单元可以有效解决框架结构的几何非线性问题和压弯耦合的材料非线性问题,使用较少的单元数量就可以得到较准确的计算结果;本文所采用的改进措施切实有效,保证了程序的稳定性.     

Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.     

O(T2)精度下椭球界面Dirichlet边值问题的积分解

研究了边界是参考椭球面的Laplace方程Dirichlet边值问题的求解,在O(ε4·T)精 度下给出了参考椭球界面上扰动重力位Dirichlet外问题的积分解式. 该结果理论上优于目 前常用的球近似下的积分解式,从而为研究物理大地测量中边值问题的求解提供了新的依据     

微分求积法在结构动力分析中的应用

微分求积法（DQM）是1种求解微分方程初（边）值问题的数值方法,通常以较小的计算工作量即可获得较高的数值精度。这种方法应用于工程领域时多用来解决梁、板等结构的静力分析或结构特征值分析等问题,即对边值问题的微分方程的求解。结构动力分析属于初值问题,荷载和结构反应都具有特殊性,直接套用DQM求解边值问题并不能获得问题的解。本文尝试利用微分求积原理建立求解结构动力反应的具体方法。借鉴单元法的思想,将荷载持时划分为若干个时步,在每个时步内对动态荷载和结构反应进行离散,然后用DQM对时步逐个进行求解,得到体系在整个时域内的反应过程。通过对3种不同自振周期的线弹性单自由度体系在不同频率简谐激励下反应的计算,阐释了本文方法的可行性以及高精度、高效率的特点,通过数值试验确定了时步内相对较优的节点数,并为时步长度的选取提供了建议。     

Definition and interests of reciprocity and reciprocity gap principles for groundwater flow problems

We introduce the reciprocity and reciprocity gap principles for flow problems in hydrogeology and illustrate their interest in addressing identification problems. The reciprocity principle is derived from mechanics and establishes for flow problems a relationship between different sets of forcing terms, including sources, sinks and boundary conditions, and the resulting head fields. The reciprocity gap principle compares different head fields resulting from the same forcing terms applied to different structures. We give general 2D expressions of the reciprocity and reciprocity gap principles for transient flow problems and give two examples of applications for the identification of transmissivity values and interfaces between different transmissivities. Identification capacities of the reciprocity and reciprocity gap principles yielding direct inversion methods could be used as initial guesses for more advanced inverse problem methodologies.     

3-D rheologic model of earthquake preparation (II): Strain field and its applications

Introduction Based on the elastic theory of the hard inclusion (Dobrovolskii, 1991), we developed an inclusion theory of rheologic medium, and applied the results of bulk-strain field of a rheologic inclusion model to explain the spatial-temporal evolution process of earthquake precursors (SONG, et al, 2000). In the former paper (SONG, et al, 2003), we derived the viscoelastic displacement field of the rheologic inclusion model on the basis of the analytic expression of displacement field o…     

气候模式模拟中海量数据的并行抽取与可视化分析诊断

在地球动力学和气候模拟等领域, 数值模拟产生的数据规模达到Tb至Pb量级。 实现这些海量数据的实时可视化和实时诊断分析面临很大的困难, 时空多尺度数据抽取可以解决这一瓶颈。 高精度数据的可视化结果展示需要高分辨率的显示设备, 并行的大屏幕显示技术是解决这一问题的有效手段。 地球科学数值模拟、 并行数据抽取和高分辨率显示都需要搭建高性能计算机集群。 本文在基于Lagrange插值的多维度、 多尺度、 多分辨率并行数据抽取算法的同时, 利用并行计算节点及LCD显示器, 基于Rocks cluster系统搭建起一个176核, 4×10×1024×1280分辨率的高性能计算模拟、 数据抽取和并行显示输出的集成平台, 并将该平台成功应用到气候模式模拟产生的海量数据的并行抽取和并行显示。