油气勘探中的碳酸盐岩古岩溶研究方法综述

检索了国内外相关文献,重新阐述了古岩溶的概念,同时介绍了岩溶研究进展情况。在此基础上,对国内外学者关于古岩溶研究的方法和技术手段进行总结,分析了几种主要方法的优点和不足之处,并对今后岩溶研究思路和方法提出建议。      

秋季太湖悬浮物高光谱估算模型

2004年10月,在太湖设置67个采样点,现场测量了水体遥感反射比、后向散射系数和辅助参数,实验室分析了水体的悬浮物浓度和水色要素吸收系数.本文在对遥感反射比光谱分析的基础上,比较了几种水体遥感反射比光谱估算悬浮物浓度的方法,结果发现广泛应用的带比值项的算法,虽然可以应用于太湖总悬浮物浓度估算,但是普遍存在相对误差较高的弱点.通过对光谱的分析,确定了750hm单波段算法的参数,并提出了利用近红外812nm波峰高度来估算水体悬浮物的方法.文中还详细阐述了为什么比值算法等在太湖水体悬浮物浓度估算中相对误差比较高,并解释了利用近红外估算精度高的原因,并指明以上算法在遥感应用中的优点和不足,以及现实水体(水质)遥感对现有遥感器配置和遥感算法的需求.     

基于解析信号重构的同步挤压小波变换的时频谱影响因素分析

受测不准原理的制约,小波变换、S变换等时频分析算法无法同时获得高时间分辨率和高频率分辨率。为了满足更高的要求,出现了一种联合小波变换和时频谱重组的新方法—同步挤压小波变换。本文从同步挤压小波变换和基于解析信号重构的同步挤压小波变换的原理出发,通过模型分析算法中参数设置对时频分析结果的影响,包括小波母函数、小波母函数的参数选择和小波阈值等,分析瞬时频率变化率不为零的信号所存在的时频谱模糊现象,并通过控制小波母函数、小波母函数参数以及小波阈值有效地减轻瞬时频率变化率不为零的信号所存在的时频谱模糊现象,时频谱的质量得到一定程度的提高。研究结果对获取高分辨率地震时频谱具有一定的指导意义。      

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.     

Shales in the Qiongzhusi and Wufeng–Longmaxi Formations: a rock-physics model and analysis of the effective pore aspect ratio

The shales of the Qiongzhusi Formation and Wufeng–Longmaxi Formations at Sichuan Basin and surrounding areas are presently the most important stratigraphic horizons for shale gas exploration and development in China. However, the regional characteristics of the seismic elastic properties need to be better determined. The ultrasonic velocities of shale samples were measured under dry conditions and the relations between elastic properties and petrology were systemically analyzed. The results suggest that 1) the effective porosity is positively correlated with clay content but negatively correlated with brittle minerals, 2) the dry shale matrix consists of clays, quartz, feldspars, and carbonates, and 3) organic matter and pyrite are in the pore spaces, weakly coupled with the shale matrix. Thus, by assuming that all connected pores are only present in the clay minerals and using the Gassmann substitution method to calculate the elastic effect of organic matter and pyrite in the pores, a relatively simple rock-physics model was constructed by combining the self-consistent approximation (SCA), the differential effective medium (DEM), and Gassmann’s equation. In addition, the effective pore aspect ratio was adopted from the sample averages or estimated from the carbonate content. The proposed model was used to predict the P-wave velocities and generally matched the ultrasonic measurements very well.     

基于比值均方根的杨氏模量反演方法

杨氏模量和泊松比是表征岩石脆性、评价储层可压性以及描述储层流体特征的重要参数,叠前地震反演是从地震数据中得到此类参数的有效手段.常规的叠前地震反演方法多使用Zeopppritz方程的近似式计算反射系数,近似式的误差影响了反演结果的精度.针对这个问题,推导了比值均方根形式的杨氏模量、泊松比和密度的Zoeppritz方程,并基于广义线性反演构建了杨氏模量和泊松比直接反演方法,有效地避免了使用近似式的局限性,提高了反演的精度.模型测试和实际资料试算结果表明：采用比值均方根形式构建的反演流程稳定,能够从叠前地震数据中获得可信的杨氏模量、泊松比和密度,提供了一种可靠的杨氏模量和泊松比直接反演方法.     

基于萤火虫和蝙蝠群智能算法的瑞雷波频散曲线反演

反演瑞雷波频散曲线能有效获取地层横波速度和厚度.但由于其高度的非线性、多参数、多极值等特点,传统的全局搜索方法易出现收敛速度慢、早熟收敛及搜索精度低的问题.鉴于此,本文提出并测试了基于萤火虫优化算法（FA）和带惯性权重的蝙蝠优化算法（WBA）的新的瑞雷波频散曲线反演策略.在瑞雷波频散曲线反演中,FA全局搜索能力强,但后期搜索精度低,而WBA局部搜索能力强,搜索精度高,但易出现早熟收敛.故本文将二者结合,提出了一种新的优化策略,称其为WFBA,即在反演前期使用FA,后期使用WBA,很好地解决了FA后期搜索精度低及WBA早熟收敛的问题.本文首先反演了三个典型理论模型的无噪声、含噪声的数据,验证了WFBA对瑞雷波数据反演的有效性与稳定性.然后将WFBA与WBA、FA单独反演以及不含惯性权重的FBA和粒子群优化算法（PSO）反演的结果进行了对比,说明了WFBA相对于WBA、FA、FBA和PSO具有更稳定、收敛速度更快、求解精度更高等优点.最后,反演了来自美国怀俄明地区的实测资料,检验了WFBA对瑞雷波数据反演的实用性.理论模型试算和实测资料分析表明,WFBA很适用于瑞雷波频散曲线的定量解释,具有很高的实用性价值.     

文莱-沙巴盆地深水褶皱冲断带构造变形特征及成因机制

为揭示活动陆缘深水褶皱冲断带的特征及成因,本文利用地震和区域地质资料的综合分析,系统阐述了文莱—沙巴盆地深水褶皱冲断带的构造变形特征,并结合盆地演化动力学特点,探讨其构造变形机制及其对深水区油气成藏的影响.研究结果表明,文莱—沙巴盆地深水褶皱冲断带具有"垂向分期、平面分段"的特点,垂向上,以中中新统底界面为界可划分为下部（始新世-早中新世）和上部（中中新世-现今）两套逆冲褶皱冲断体系,其中下部逆冲褶皱冲断带的形成与古南海的俯冲作用密切相关,上部逆冲褶皱冲断带是中中新世以来三角洲前缘重力滑动与苏禄海扩张造成的区域挤压应力远程效应共同作用的结果,且苏禄海扩张造成的远程挤压效应主控平面上南北段褶皱冲断带变形的差异性,导致北段褶皱变形强度大于南段,具有背斜褶皱数量多、褶皱间距离短、逆冲断层倾角陡的特点,南段反之;且晚上新世以来北段深水区地层缩短量大于陆架区伸展量,两者之差为2~6 km,而南段两者相当,仅受三角洲前缘重力滑动影响.整个褶皱冲断带发育断弯、断展、断滑褶皱等3种断层相关褶皱以及叠瓦扇和冲起构造2种逆冲构造组合,是多期NW向挤压应力作用下形成的大型逆冲推覆构造,以前展式向盆地扩展.此外,由于中中新世以来逆冲断层的持续活动,研究区深水褶皱冲断带发育众多构造圈闭,油气成藏条件优越,且南段优于北段,靠近陆坡的近端优于远端,可作为勘探部署重点.     

地磁活动对气候要素影响的研究进展

地磁活动是太阳爆发现象引起地球近地空间磁场扰动的重要空间天气过程之一.地球磁场的变化具有多种时间尺度,其中从数十年到数世纪的长时间地磁场变化主要是由地核磁场引起的,而从数秒到数年的短时间地磁变化与太阳活动有关.近年来,越来越多的统计研究表明,地磁活动与太阳活动和地球气候变化之间存在着显著的相关性.地球磁场和地球大气系统的耦合现象驱动着人们探索地磁活动对地球天气和气候系统影响的研究.本文的目的就是综述国内外地磁变化对气候影响的研究进展,介绍我们最新的研究成果,探索地磁活动对气候要素的影响特征和可能机理过程,为深入研究地磁活动对地球天气和气候的影响提供基础和依据,以期对地磁活动和气候要素关系有进一步的认识.     

面向目标自适应三维大地电磁正演模拟

本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.