洞庭盆地第四纪沉积环境演变的初步分析

利用水流能量指数，分析了洞庭盆地第四纪沉积环境的演变特征。得出洞庭盆地主要为河流－过流型湖泊相互消长的沉积环境。沉积环境的演变及主要水流方向，受差异性构造沉降的控制。历史时期，人为作用对盆地沉积环境影响很大，并且是近代洞庭湖区洪灾加剧的主要原因。     

槽波地震勘探的数学模型研究

为研究槽波的形成及地震时间剖面上的表现规律，利用非均匀各向同性介质的波动方程，应用有限差分法进行求解并编制了计算程度，给出了一个3层不对称模型的两分量计算结果实例。     

含油水各向异性孔隙介质中地震波传播方程

油气勘探的重点是寻找油气储层的分布 ,而地下油气的存在必然对地震记录产生影响。传统的地震波理论是建立在纯弹性固体基础之上 ,没有考虑固体中所含的流体 ,因此 ,它难以研究含油水岩石的物理性质 ,如岩石的渗透率和孔隙度、油水粘度。从固、液介质系统能量和本构方程出发 ,推导出固体、油相和水相的动力学方程 ,进而建立起含油水两相流体各向异性孔隙介质的地震波传播方程。该方程包括有岩石和油水的物理参数 ,更适合于油田的勘探和开发。     

孝妇河源头地区岩溶地下水资源可持续利用的决策支持研究

基于运筹学和系统工程学的基本理论，建立了地下水最优控制模型，并采用微分动态规划方法进行求解，最后直接根据淄博市孝妇河源头地区地下水资源的可持续利用需求，建立了该区的地下水最优控制模型，并求得了最优解。     

Re-Os同位素体系在陨石研究中的应用

铁陨石中的Re ,Os含量反映其结晶分异历史。通过铁陨石定年修正187Re的衰变常数为 :λ(187Re) =1 6 6 6× 10 -11a-1。ReOs同位素测年法可以直接用于对铁陨石的定年 ,结果表明天然铁陨石大体同时形成 ,但ReOs定年技术已有可能揭示不同化学群铁陨石形成年代的序列 ,但研究尚需深入。这些方法也可以用来探讨铁陨石和石铁陨石的形成源区、冷却历史和后期变化。虽然在石陨石中Re ,Os同位素的浓度很低 ,但也有了探索性研究成果。随着技术的不断发展 ,ReOs同位素体系在天体化学中的作用将愈加明显和重要。     

Plasticity model for sand under small and large cyclic strains: a multiaxial formulation

This paper presents the multiaxial formulation of a plasticity model for sand under cyclic shearing. The model adopts a kinematic hardening circular cone as the yield surface and three non-circular conical surfaces corresponding to the deviatoric stress ratios at phase transformation, peak strength and critical state. The shape of the non-circular surfaces is formulated in accordance with the experimentally established failure criteria, while their size is related to the value of the state parameter ψ. To simulate cyclic response under small and large shear strain amplitudes without a change in model parameters, it was found necessary to introduce: (a) a non-linear hysteretic (Ramberg–Osgood type) formulation for the strain rate of elastic states and (b) an empirical index of the effect of fabric evolution during shearing which scales the plastic modulus. This index is estimated in terms of a macroscopic second-order fabric tensor, which develops as a function of the plastic volumetric strain increment and the loading direction in the deviatoric plane. Comparison of simulations to pertinent data from 27 resonant column, cyclic triaxial and cyclic direct simple shear tests provide a measure for the overall accuracy of the model.     

The effect of ferromagnetism on the equation of state of Fe3C studied by first-principles calculations

First-principles calculations have been used to determine the equation of state of Fe₃C in both its low-pressure magnetically ordered and high-pressure non-magnetically ordered states; at 0 K the ferromagnetic transition was found to occur at about 60 GPa. In the high pressure, non-magnetically ordered regime at 0 K the material may be described by a Birch-Murnaghan third-order equation of state with V₀=8.968(7) ų per atom, K₀=316.62(2) GPa and K′=4.30(2). At atmospheric pressure the ferromagnetic phase transition in Fe₃C occurs at ∼483 K; preliminary measurements of the thermal expansion by powder neutron diffraction show that this transition produces a large effect on thermoelastic properties. The volumetric thermal expansion coefficient in the paramagnetic phase was found to be 4.34×10⁻⁵ K⁻¹ at T∼550 K. By applying a thermal expansion correction to the calculated equation of state at 0 K, predicted values for the density and adiabatic incompressibility of this material at core pressures and temperatures were obtained. These results appear to be sufficiently different from seismological data so as to preclude Fe₃C as the major inner core-forming phase.     

Equation of state of gold and its application to the phase boundaries near 660 km depth in Earth’s mantle

The pressure-volume-temperature equation of state (EOS) of gold is fundamental to high-pressure science because of its widespread use as an internal pressure standard. In particular, the EOS of gold has been used in recent in situ multi-anvil press studies for determination of phase boundaries related to the 660-km seismic discontinuity. These studies show that the boundaries are lower by 2 GPa than expected from the depth of the 660-km discontinuity. Here we report a new P-V-T EOS of gold based on the inversion of quasi-hydrostatic compression and shock wave data using the Mie-Grüneisen relation and the Birch-Murnaghan-Debye equation. The previously poorly constrained pressure derivative of isothermal bulk modulus and the volume dependence of Grüneisen parameter (q=d lnγ/d ln V) are determined by including both phonon and electron effects implicitly: K′_0T=5.0±0.2 and q=1.0±0.1. This combined with other accurately measured parameters enables us to calculate pressure at a given volume and temperature. At 660-km depth conditions, this new EOS yields 1.0±0.2 GPa higher pressure than Anderson et al.’s EOS which has been used in the multi-anvil experiments. However, after the correction, there still exists a 1.5-GPa discrepancy between the post-spinel boundary measured by multi-anvil studies and the 660-km discontinuity. Other potential error sources, such as thermocouple emf dependence on pressure or systematic errors in spectroradiometry, should be investigated. Theoretical and experimental studies to better understand electronic and anharmonic effects in gold at high P-T are also needed.     

在构造和气候因素制约下雅鲁藏布江的演化

雅鲁藏布江位于印度和欧亚大陆汇聚带内，其形成受到冈底期山和喜马拉雅山差异性抬升的控制。冈底期山抬升在先，发生在中生代晚期至新生代早期。一系列起源于冈底期山和青藏高原的水系向南先是流主特提斯海。在特提斯海关闭后流入印度次大陆。喜马拉雅山构造抬升要晚于冈底斯山，大规模抬升发生在中新世早期，其抬升阻断了这些河流的通道，水流开始汇聚在这两个造山带之间，牙鲁藏布江由此形成。在雅鲁藏布江大拐弯地区，在海拔4500m处存在一个平坦的侵蚀面，并构成雅鲁藏布江大峡谷最高的一级谷肩，这表明雅鲁藏布江在下切前就在该面上流动，而且流速不大。在大拐弯以南，雅鲁藏布江的下游－布拉马普特拉河位于印度洋热带季风带内，其下切和源侵蚀速率很大。印度洋热带季风形成于6－9MaB.P.。因此，该河流很可能形成于该时期，要比雅鲁藏布江年轻，它在向北的溯源侵蚀的过程中袭夺了雅鲁藏布江，袭夺处可能就大拐弯的北端，因此大拐弯是袭夺成因。     

水气二相流方程的一种离散数值解法

提出了水气二相流方程的一种数值解法.在利用有限元方法求解水气二相流方程时,引入了离散Newton迭代方法,用于非线性有限元方程组的线性化处理,将这一步计算的收敛阶由原有研究的线性收敛提高到平方收敛,并避免了直接应用Newton迭代方法给编程带来的不便.同时在求解两相的有限元方程组时,采用两相方程组并行迭代的方法,与联立计算相比节省了大量的内存空间.