莺歌海盆地东方区高温超压气藏盖层封盖机制

莺歌海盆地是一个典型的高温超压盆地,盖层能否及如何有效封闭该盆地内中深层高温超压气藏是一个值得探讨的重要问题.为此,以该盆地东方区为例,通过气田实际资料及理论分析盖层的宏观地质特征及微观封闭能力.研究表明:该区高温超压气藏盖层的封盖机制包括低构造倾角、厚层泥岩封闭、浅海泥岩盖层中泥质粉砂岩与粉砂质泥岩频繁互层、超压间接封闭、高温超压引起的高突破压力与排替压力比值;其中前3者是决定该区高温超压气藏盖层对天然气封闭能力的主导因素和先决条件,后两者则是中深层浅海相泥质碎屑岩盖层具备优质封盖条件的基础.进而对比剖析了DF-A、DF-B超压气田及DF-C含气构造等典型实例,认为:DF-B气田气体充满度高的原因在于上述五个要素的最优配置;而高构造倾角的岩性气藏储层能量相对较高,对盖层封闭能力的要求高,在其他盖层条件相似甚至更差的情况下,气藏遭到破坏的风险大,这也是DF-C含气构造尚未取得商业性突破的重要原因.     

重构小波阈值函数在信号去噪中的应用与研究

实际工程应用中,通过物理手段采集到的信号都带有噪声信息,这样就会淹没很多有用信号,传统的“带通、低通、高通”滤波技术往往显得无能为力,因此提取有用的特征信息就需要对原始信号进行去噪。常用的小波阈值去噪（硬、软阈值函数）方法,含有诸多缺点,本文在其基础上重新构造了一个新的阈值函数。在MTALAB（2014a）环境下,用传统的硬、软阈值函数及重构阈值函数对混入高斯白噪声信号进行去噪仿真分析,结果表明,重构的阈值函数处理的信号效果更好,更清晰。     

局部倾角约束最小二乘偏移方法研究

随着石油勘探难度的进一步加大,地震数据往往存在采样不规则、地震道缺失等现象,如果不对其进行处理,会对后续的地震成像产生影响,引入成像噪音.针对这一问题,一般是通过地震道插值或数据规则化对叠前数据进行处理,然后采用常规的偏移方法进行成像,本文则是将地震成像看作最小二乘反演问题,在共成像点道集引入平滑算子,在共偏移距/角度道集引入平面波构造算子(PWC)进行约束,通过预条件共轭梯度法使得反偏移后数据与输入数据之间的误差达到最小,最终得到信噪比更高、振幅属性更为可靠的成像结果.理论模型和实际资料处理表明,本文方法不仅可以有效压制数据不规则对成像产生的噪音,而且具有更高的成像精度.     

全国自动地震速报系统评估

对全国自动地震速报系统进行了初步的评估. 通过2010年1月1日—2011年12月31日的自动定位结果与全国正式目录的对比, 分别得出了3套自动定位系统总体误差以及全国和全球的误差插值图, 并对结果进行了总结. 同时, 还针对3套自动定位系统的漏报和误报情况进行了分析和讨论. 总体来说, 现阶段自动地震速报系统整体上对于首都圈地区M≥3.0, 国内其它地区M≥4.0(部分台站稀疏的西部地区M≥5.0), 国外M≥7.0的地震基本可以达到全面覆盖的程度, 可以作为正式地震速报的参考, 对于建立“自动报—正式报”的两级地震速报制度能提供有力的技术支撑.     

基于结构导向的梯度属性边缘检测技术

地震数据中的不连续性检测对于识别和刻画地下地质体边界具有重要意义.随着勘探目标向岩性油藏和复杂断块油藏的转变,解释人员需要对小断裂、微裂缝等微幅地质体边界信息进行精细刻画.常规边缘检测算子由于算法的局限性不能满足上述需求.本文通过引入能量归一化和沿层检测思想对常规检测算子进行改进,提出一种基于结构导向的梯度属性边缘检测技术.该技术通过结构张量估计得到局部地层区域的倾角和方位角,利用得到的倾角和方位角信息作为约束,使用二维边缘检测算子实现三维数据中的沿层不连续性检测.经过实际资料处理,并与结构导向方差属性对比分析,证实该方法在地层不连续性检测方面具有很好效果,检测得到的不连续性信息在细节上表现更丰富,可以作为地层横向差异和地质体边界识别的有力工具.     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

鲜水河断裂带北西段跨断层基线变化机理探讨

2013年芦山M_S7.0地震发生后,鲜水河断裂北西段的侏倭、虚墟等跨断层测点基线测项均出现不同程度的异常变化,异常核实结果认为基线异常属实.该异常变化是否与芦山地震的发生存在一定关联性,值得深入研究.本文首先分析芦山地震发生前后鲜水河断裂北西段不同基线长度的变化特征.基于鲜水河断裂带及邻区三维非线性黏弹性有限元模型,以GPS观测资料、同震静态滑移量作为约束,通过开展多组数值模拟实验,探讨下地壳不同流变特征下,芦山地震发生前后不同阶段,不同跨断层基线长度的时序动态变化特征.初步研究结果表明,1）侏倭、虚墟基线测点原始观测资料均表现出地震发生时基线长度减小,发生后短时间内反向快速恢复及增加的协同变化;2）在以GPS观测数据为约束的块体间差异性运动的动力学边界条件下,鲜水河断裂表现为左旋走滑运动;3）模拟结果显示,地震发生时测点间距离迅速减小,表现为断层右旋的特征,这与实际观测资料反映的结果一致;4）鲜水河断裂带两侧地块的下地壳黏滞系数分别取10¹⁸Pa·s、10¹⁹ Pa·s时,由于黏弹性松弛效应的影响,造成基线长度在地震发生后短时间尺度内快速增加.对比分析认为,地震后短时间内实测资料反映的基线长度快速增加的特征可能是黏弹性松弛效应与构造作用共同作用的结果.

     

BP人工神经网络在MM5预报福建沿海大风中的释用

利用变分方法建立预报场和预报倾向场这一预报场组合与模式预报非系统性误差之间的映射关系,来估计GRAPES（Global/Regional Assimilation and PrEdiction System）模式的非系统性误差,从而对预报做出修正。采用两种不同的历史样本建立这一映射关系,其中,利用相同时刻历史样本建立映射关系的方法称为DEM方法;通过相似面积比选取“相似样本”来建立上述映射关系的方法称为SEM方法。以FNL分析资料作为评判预报误差的依据,根据2002—2010年7月GRAPES模式500 hPa高度场48 h预报的回报资料,利用两种不同的方案进行非系统性误差的估计及预报订正试验。对279个检验样本的试验结果表明：SEM方法和DEM方法都对非系统性误差有一定的估算能力,二者估算的非系统性误差空间分布和量级与模式非系统性误差较一致,SEM方法的修订效果略优于DEM方法,但并不明显。对预报做出系统性误差和非系统性误差两步订正后,DEM方法和SEM方法的订正有效率分别为98.566%和100%,可明显提高预报的准确性。

     

上海TWP3型边界层风廓线雷达探测性能评估

基于2018年7月—2019年6月上海宝山移动式风廓线雷达（mobile wind profiler radar,MWPR）资料,从获取率、水平风误差和上海6类典型天气条件下测风数据的可靠性三方面,对其探测性能进行了评估与分析。结果发现：（1）获取率受湍流强度影响最大,受大气温度、湿度影响次之。获取率具有较明显的日变化和季节变化特征,边界层内夏秋季获取率高于冬春季,且有午后获取率高、深夜和清晨低的特点;自由大气中冬春季获取率高于夏秋季。以获取率80%为其业务准入标准,边界层内MWPR的获取率达标。（2）相比L波段探空测风数据,MWPR整层的水平风速均偏小,低模态时有34.4%的数据位于偏小2.0~6.0 m·s^-1区间内;风向偏差超过±15°以上的水平风向数据占比为48.9%。（3）定性分析结果显示,MWPR探测到的垂直速度的指向性是合理的。雨强较大时,MWPR探测的垂直速度实际上是降水粒子的下降速度与大气湍流速度之和。MWPR的水平风数据在台风天、雪天等大气均一性强的时段内可靠性高。总体上,上海宝山MWPR的探测性能基本能满足日常监测服务之需,但其测风精度与业务准入标准相比尚有一定差距。

     

均方误差意义下的正则化参数二次优化方法

Tikhonov正则化法是大地测量中应用最为广泛的病态问题解算方法之一。影响正则化法解算效果的重要因素是正则化参数,然而,最优正则化参数的确定一直是正则化解算的难题,如L曲线法确定的正则化参数具有稳定性好、可靠性高的优点,但存在过度平滑问题,导致正则化法对模型参数估值精度改善较小。本文从均方误差角度分析了正则化参数对模型参数估计质量的影响。基于奇异值分解技术,提出了由模型参数投影值分块计算均方误差的方法,避免了均方误差迭代计算,并基于均方误差最小准则给出了正则化参数优化方法,实现了对L曲线正则化参数的优化。数值模拟试验与PolInSAR植被高反演试验结果表明,正则化参数优化方法有效改善了正则化法解算效果,提高了模型参数估计精度。