乳山金矿煌斑岩及流体和氢、氧、锶同位素研究

乳山金矿床是胶东近年来探明的一个赋存在花岗岩中的石英脉型金矿床。矿出露的煌斑岩金含量较低，不可能是金矿化的重要物源者。     

A mixed finite element method for Hele-Shaw cell equations

A new numerical method to solve the system of equations describing two phase flow in a Hele-Shaw cell is presented. It combines a mixed finite element method, the method of subtraction of the singularity and a front tracking grid in a single computational strategy. This choice of discretization techniques is well motivated by the difficulties present in the system of equations and the physics of the problem. The new method was tested against analytical solutions and also by solving the Saffman–Taylor viscous fingering problem for finite and infinite mobility ratios. In both cases convergence under mesh refinement is achieved for the fingers developed from an initial sinusoidal interface. Finger splitting is observed for low values of the surface tension and high mobility ratio. Different explanations, based in our results, are provided for this phenomenon. This revised version was published online in July 2006 with corrections to the Cover Date.     

基于动态图像变形的PP与PS波层位直接匹配

多分量地震资料的矢量偏移、多波地震资料的联合解释与反演均需要估算纵横波的速度比,实现纵波与转换横波在时间或深度域的匹配.基于DTW,本文实现了一种适用于PP与PS波直接匹配的动态图像变形算法.该算法分为三个部分：首先,使用二阶对称动态规划算法逐样点递归计算PP与PS波走时或深度的误差累积和;其次,在以误差累积和为目标函数的回溯阶段设定变形窗,并在纵横波速比约束的变形窗内递归回溯搜索匹配路径;最后,根据最大相关系数判定准则在匹配路径中确定最佳匹配路径,获得使PP与PS波匹配的拉伸或压缩时移量.利用所获得的拉伸压缩时移量计算纵横波速度比就可以实现PP与PS波之间的匹配.模型与实际陆上多分量地震资料测试结果表明：该方法具有较高的匹配精度,且对于信噪比、相似度较低的多分量地震资料,该方法也能产生较好的匹配效果.

     

Deposition Rates of Terrestrial and Marine Organic Carbon in the Osaka Bay, Seto Inland Sea, Japan, Determined Using Carbon and Nitrogen Stable Isotope Ratios in the Sediment

Carbon and nitrogen stable isotope ratios (¹³C and ¹⁵N) of surface sediments were measured within Osaka Bay, in the Seto Inland Sea in Japan, in order to better understand the sedimentation processes operating on both terrestrial and marine organic matter in the Bay. The ¹³C and ¹⁵N of surface sediments in the estuary of the Yodo River were less than –23 and 5 respectively, but increased in the area up to about 10 km from the river mouth. At greater distances they became constant (giving ¹³C of about –20 and ¹⁵N about 6). It can be concluded that large amounts of terrestrial organic matter exist near the mouth of the Yodo River. Stable isotope ratios in the estuary of the Yodo River within 10 km of the river mouth were useful indicators allowing study of the movement of terrestrial organic matter. Deposition rates for total organic carbon (TOC) and total nitrogen (TN) over the whole of the Bay were estimated to be 63,100 ton C/year and 7,590 ton N/year, respectively. The deposition rate of terrestrial organic carbon was estimated to be 13,200 (range 2,000–21,500) ton C/year for the whole of Osaka Bay, and terrestrial organic carbon was estimated to be about 21% (range 3–34) of the TOC deposition rate. The ratio of the deposition rate of terrestrial organic carbon to the rate inflow of riverine TOC and particulate organic carbon (POC) were estimated to be 19% (range 3–31) and 76% (range 12–100), respectively.     

Hydraulic parameters in channels with wall vegetation and gravel bed

In this experimental study,field observations and laboratory experiments have been carried out to assess the impacts of the vegetated channel walls and aspect ratio on flow velocity profiles,shear stress distribution and roughness coefficient of channel.Results show that the presence of vegetation cover on channel wall causes deviation of the Reynolds stress distribution from the linear one under uniform flow condition.It is also noticed that the Reynolds stress distribution is influenced by the aspect rati...     

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.

     

绝对重力仪激光器光纤对准调节方法概述

绝对重力仪激光器光纤由于应力变化、形状改变、震动或耦合器维护等原因破坏了输出激光的偏振态,使其在光纤内的传输变得非常不稳定,导致光纤输出激光的消光比值小于100,无法满足绝对重力观测条件。针对该问题,介绍了光纤对准调节的方法,通过设计消光比检测装置检验消光比值是否大于100。对#181和#250两台WEO100型He-Ne激光器的光纤进行对准调节实验,并通过消光比检测装置实测,两台激光器输出的光束消光比值满足绝对重力测量的条件,表明该对准调节方法科学有效。     

川东北地区三叠系飞仙关组白云化流体的Ca/Mg比定量估算

作为川东北地区三叠系飞仙关组天然气气藏最重要储集岩的粒屑白云岩和结晶白云岩的成因始终是沉积和石油地质学家高度关注的研究热点之一.在总结川东北地区三叠系飞仙关组不同类型白云岩形成温度的基础上,通过已有的白云化流体Ca/Mg值和温度之间关系式获得了川东北地区三叠系飞仙关组白云化流体的Ca/Mg值的主要分布范围.结果表明:不同类型白云岩对应的白云化流体具有不同的Ca/Mg值,同时证实了三叠纪海水或与之有关的海源流体、三叠纪海水与大气淡水混合流体的Ca/Mg值都符合在40～130℃区间、甚至更高温度下白云化(白云石沉淀)作用的要求,微晶白云岩和具原始结构的粒屑白云岩可能与地表—近地表的低温、小Ca/Mg值混合流体中相对快速的白云石沉淀作用有关,并形成了相对细的他形白云石;结晶白云岩(包括结构类似的过渡岩石类型、白云石胶结物)可能与埋藏的高温、大Ca/Mg值三叠纪海水或与之有关的海源流体中相对缓慢的白云石沉淀作用有关,并形成了相对粗的自形白云石.白云化流体的Ca/Mg值估算,可望为川东北地区三叠系飞仙关组优质储层白云岩的成因研究提供一些有价值的基础资料.     

帕米尔东北缘地壳结构的P波接收函数研究

利用位于新疆帕米尔东北缘地带12个固定数字地震台和天山动力学Ⅱ期10个流动宽频带数字地震台记录的高质量远震波形数据,应用接收函数H-κ叠加方法研究了帕米尔东北缘的地壳厚度-泊松比特征和部分台站下方的壳内界面深度.研究发现:(1) 帕米尔东北缘的Moho面起伏变化剧烈,其总体分布呈现东薄西厚、南厚北薄的特征,由塔里木盆地向天山延伸,地壳厚度约从45 km加深到55 km,从塔里木盆地向西昆仑山延伸,地壳厚度约从45 km加深到69 km;(2)沿着天山动力学Ⅱ期剖面,位于塔里木盆地北缘台站的壳内间断面的深度约为13~16 km,向北进入天山南麓加深到20 km左右,继续向北进入南天山山区壳内间断面不明显,可能暗示塔里木盆地基底向北俯冲,俯冲距离可能到达南天山的山前;(3)研究区地壳泊松比变化复杂(约从0.20到0.31),显示地壳物质组成的复杂性和显著的不均匀构造;(4)整个研究区的地壳厚度和泊松比之间没有明显的相关性,但天山动力学Ⅱ期剖面的结果表明,从塔里木盆地北缘到西南天山,地壳厚度和泊松比之间存在反相关关系,意味着天山地壳的增厚可能主要是通过以长英质岩石为主要组成成分的上地壳叠置而成;(5) 研究区全部地震台地壳厚度与海拔高程的线性回归方程表明地壳厚度与海拔的相关性相对较弱(相关系数为0.66),天山动力学Ⅱ期10个台站的地壳厚度与海拔具有很好的相关性(相关系数为0.85),可能表明沿该剖面地壳整体上处于相对均衡的状态.     

A simple critical state interface model and its application in prediction of shaft resistance of non-displacement piles in sand

A simple semi-hyperbolic state-dependent constitutive model for sand-structure interfaces is proposed. The model formulation is consistent with critical state soil mechanics since void ratio evolves continuously with shear strain from initial state towards asymptotic critical state at extremely large shear strains. The model takes into account influence of normal stiffness on volume change and stress path. The proposed interface model is implemented in a pile segment analysis scheme for simulation of shaft resistance mobilization in non-displacement piles. Results reveal that the proposed pile segment analysis can well predict shaft resistance of model piles embedded in different sands.