海相火山岩区地质填图方法——以北祁连山东段白银矿田古火山机构厘定为例

20世纪90年代以来,火山机构控矿理论在海相火山岩区块状硫化物矿床研究中得到了越来越多的关注,基于该理论的指导,在北祁连山东段的白银矿田取得了找矿突破性进展。其间,详细的海相火山岩相地质填图工作为白银矿田火山机构与成矿联系的研究及找矿突破奠定了坚实的基础。反映出在海相火山岩区推行合理可行的火山岩相地质填图方法至关重要。以白银矿田海相火山岩区为例,首先,对该区大比例尺火山岩相填图及古火山机构厘定成果进行介绍;其次,通过分析、借鉴此成功案例,总结海相火山岩区填图内容和技术方法;最后,阐述应用该方法在北祁连山中段祁连县一带海相火山岩区进行大比例尺填图的成功尝试。     

碎石土滑坡的因素敏感性计算分析

通过资料搜集整理和分析、现场工程地质调查与勘探和室内外的物理力学试验,采用弹塑性接触有限元强度折减法、不平衡推力法和数理统计分析法,结合工程实例,分析了影响碎石土滑坡各因素的敏感性,结果表明,按因素敏感度从大到小排列,主要影响碎石土滑坡稳定性的因素依次为滑面岩土体内摩擦角、地形坡度、滑体饱水面积比和滑面岩土体的内聚力。地形坡度大、具有低内摩擦角滑面的碎石土滑坡的整体稳定性较差,而连续降雨、强降雨会抬高地下水位增大滑体饱水面积比,降低滑面带岩土体抗剪强度,增加下滑力,对碎石土滑坡的稳定性将产生最不利的影响。同时,在对碎石土滑坡进行加固处理时,必须优先考虑采取提高滑面等效内摩擦角、降低地形坡度结合降低地下水位措施,其次考虑采取提高滑面岩土体内聚力的措施。     

土样比重的测量不确定度评定

土样比重的测量不确定度既判定了自身的可靠性,同时又影响了通过比重及其他一些参数换算得到的孔隙率、饱和度等指标的准确性。选用物理性质较为均匀的原状粉土做样品,基于《土工试验方法标准》(GB/T50123-1999),通过比重瓶法重复测量了样品的比重。根据《测量不确定度评定与表示》(JJF1059-1999)的要求,分析了试验过程中影响因素的主要来源,建立了测量比重不确定度的计算模式,较系统地对其不确定度进行了评定。结果表明:土的比重测量误差主要由称量及不均匀性两个因素引起。试验条件下的土样比重测量结果的扩展不确定度为0.005。     

固体介质中孤立波的传播及演化特征

由于地球介质中广泛存在断裂、微裂缝等地质现象,实际地震资料中会出现形似孤立波的非线性地震现象。因此,对固体介质中孤立波的研究有利于解释这些非线性地震现象的形成机制。本文基于KdV方程,以雷克子波作为初始条件,采用伸缩子机理构建体力模型,利用有限差分的方法模拟了孤立波的演化过程。理论结果表明,非线性地震纵波可以从雷克子波逐渐演化成孤立波,而且地震波的初始振幅和频散系数对模拟结果也有重要影响。通过与实际资料的对比也能说明这种演化的可能性。同时根据方程系数矩阵中元素带状分布的特征,采用稀疏矩阵的存取方法,可以减小计算内存,提高计算效率。     

用于地球物理反演的均匀设计优化算法

将均匀设计、模拟退火和单纯形方法相结合，提出一种全局优化迭代算法．数值模拟表明，它具有搜索速度快、搜索能力强和不需要求导数的特点，可用于地球物理反演研究．     

A generalized transformation approach for simulating steady-state variably-saturated subsurface flow

A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach.     

Estimating open boundary conditions from coastal tidal observations by adjoint approach

A numerical study for estimating the tidal open boundary conditions of a shelf current modrl from tb coastal tidal observations is presented. The method is based on the optimal control/adjoint method. A lrast square fitting of the model state to simulated data is used. Two ideal domains and coastlines are considered. Using the IAP shallow. water model and its adjoint model, some identical twin experiments are carried out to test efficiency and lirnilsd of the method. The results show that the adjoint method can efficiently estimate the open boundary conditions well for gulf/bay like domains. The adjoint method seems to have great potential to improve the accuracy of tide and shelf current modeling in coastal regions. Project supported hy the National Natural Science Fuundation of China (Grant No. 49376256)     

FS方法及其在综合多项震兆中的应用

本文提出了一种模式识别方法--FS方法.该方法的特征取为某些论域上的模糊子集.先以每三个原始特征组成子分类器的特征集,“训练”这些分类器,然后对子分类器进行筛选,最后用筛选出的子分类器的“软”分类结果的加权平均作为判别函数,形成总的分类决策.将FS方法应用于综合多项震兆,并进行了一系列控制试验,还与地震学中几种常用的方法进行了对比试验.结果表明,FS方法是稳定和有效的.     

利用H/V方法研究琼北火山区的浅表层地质变化特征

海南琼北活火山作为华南沿海重要的第四纪火山, 一直受到国内外学者的广泛关注, 目前对该火山的形成演化已开展了大量深部地质结构的研究, 但对其表层地质结构特征仍然缺乏详细了解。水平-垂直频谱比法(简称H/V方法)已被广泛应用于浅表层地质结构的研究, 本文基于35个便携式微功耗地震仪组成的流动台阵观测数据, 采用经典全波场H/V_(full)和瑞利波椭圆率H/V_(Ray)分析方法, 获得了琼北火山区H/V曲线峰值频率(F₀)和放大倍数(A₀)的空间变化特征。琼北火山区H/V曲线峰值频率变化强烈, 以长流-仙沟断裂为界, 断裂西南侧很难从H/V曲线中获取清晰的F₀和A₀, 而断裂东北侧则存在明显的F₀和A₀。结合已有地质资料, 本文认为长流-仙沟断裂西南侧处于凹陷区, 其浅表层火山岩未发生明显的风化作用, 而东北侧位于隆起区域, 地表火山岩受到了比较严重的风化作用。这一研究结果表明, H/V方法可反映表层火山岩的风化状态, 火山岩出露场地H/V曲线的峰值频率和放大倍数随风化作用程度的增强而表现得更加明显。本研究为琼北火山区浅表火山岩风化状态和地震危险性评估提供了重要的信息支撑。     

A synthetic aperture radar sea surface distribution estimation by n-order Bézier curve and its application in ship detection

To dates,most ship detection approaches for single-pol synthetic aperture radar(SAR) imagery try to ensure a constant false-alarm rate(CFAR).A high performance ship detector relies on two key components:an accurate estimation to a sea surface distribution and a fine designed CFAR algorithm.First,a novel nonparametric sea surface distribution estimation method is developed based on n-order Bézier curve.To estimate the sea surface distribution using n-order Bézier curve,an explicit analytical solution is derived based on a least square optimization,and the optimal selection also is presented to two essential parameters,the order n of Bézier curve and the number m of sample points.Next,to validate the ship detection performance of the estimated sea surface distribution,the estimated sea surface distribution by n-order Bézier curve is combined with a cell averaging CFAR(CA-CFAR).To eliminate the possible interfering ship targets in background window,an improved automatic censoring method is applied.Comprehensive experiments prove that in terms of sea surface estimation performance,the proposed method is as good as a traditional nonparametric Parzen window kernel method,and in most cases,outperforms two widely used parametric methods,K and G0 models.In terms of computation speed,a major advantage of the proposed estimation method is the time consuming only depended on the number m of sample points while independent of imagery size,which makes it can achieve a significant speed improvement to the Parzen window kernel method,and in some cases,it is even faster than two parametric methods.In terms of ship detection performance,the experiments show that the ship detector which constructed by the proposed sea surface distribution model and the given CA-CFAR algorithm has wide adaptability to different SAR sensors,resolutions and sea surface homogeneities and obtains a leading performance on the test dataset.