南海东北部洋陆过渡区域地震折射波衰减特征

地震波衰减研究是了解岩石圈构造特征的有效方法.本文基于OBS2016-2测线的海底地震仪(OBS)数据, 首次对南海东北部洋陆过渡区域地震波衰减特征进行研究, 通过正演模拟获得该区域二维纵波衰减(Q_P)结构.结果表明, 下陆坡的上地壳存在一个宽约40 km, 厚度约为4~5 km的高衰减区, 其特征为低纵波和横波速度(V_P为5.5~6.3 km·s^-1和V_S为3.1~3.6 km·s^-1)以及较低的波速比(V_P/V_S为1.72~1.80), 对应较低的Q_P(280~410), 推测与断裂发育有关, 且受到火山活动的影响.洋陆过渡及洋壳区域的上地壳高衰减区具有低Q_P(300~400)和高V_P/V_S(1.90~1.96)特征, 可能对应较多的断裂发育及流体运移.洋陆过渡区域下地壳高速异常体表现为相对低的Q_P(550~600), 对应较高的V_P(7.0~7.8 km·s^-1)和V_S(3.5~3.8 km·s^-1)以及较高的V_P/V_S(1.85~1.96), 推测与蛇纹石化作用有关.蛇纹石化可能进一步增加岩石的孔隙度并导致更多的流体运移, 使得洋陆过渡及洋壳区域存在较高的地震波衰减.Q_P结构有助于我们分析南海大陆边缘的地震波衰减特征, 结合V_P、V_S以及V_P/V_S, 可以更好地了解该区域地质结构和岩石属性, 对进一步挖掘OBS数据信息有重要参考价值.

     

角度域广义Radon变换多参数非迭代线性反演

基于Born近似的线性反演是当前流行的定量化地震成像与反演方法, 然而, 现有线性反演方法存在多种不可避免的问题, 比如迭代计算成本过高, 背景模型依赖, 以及多参数提取不确定/不稳定.为避免这些问题, 本文提出了角度域广义Radon变换(AD-GRT)多参数非迭代反演方法(2D声波双参数介质).AD-GRT是一种地震数据与角度域模型之间的积分变换, 该角度域模型能够保留地震散射数据的完整信息.基于该角度域模型, 地震数据可被精确重建, 即使背景模型是错误的.介质双扰动参数可在角度域模型中准确且稳定提取.为解决传统角度域离散方法存在的幅值震荡和不连续问题, 设计了两种有效的角度域离散单元分裂方法, 实现了平滑连续的角度域幅值提取.为提高AD-GRT反演的分辨率, 考虑了震源子波反滤波.利用三个数值模型验证了AD-GRT在角度域模型反演、数据重建、多参数提取方面的有效性.

     

基于SVD和BRD的二维核磁共振测井正则化反演算法研究

二维核磁共振能从"弛豫-扩散"两个维度上展现流体性质,在识别稠油储层方面具有理论优势,是当前核磁共振测井技术的研究热点.本文深入研究二维核磁共振测井原理,系统分析CPMG-DE脉冲序列测量扩散系数与弛豫时间的方法,结合核磁共振二维谱数理模型,提出一种基于SVD和BRD的正则化反演算法.该算法通过SVD压缩数据,采用带非负约束的Tikhonov正则化方法求解流体"弛豫-扩散"分布,并基于BRD算法迭代确定最佳正则化因子.模拟实验与数值分析表明,该算法无需先验信息、运算效率高、相对误差小,在原始数据信噪比低至50时,仍可有效获取流体(T2,D)二维分布.在二维核磁共振测井数据实时解释应用中,该方法较传统反演算法(如TSVD)具有较大优势.同时,在自主研发的核磁共振测井仪测量CuSO4溶液(T2,D)分布的实验显示,本文设计算法对弛豫时间和扩散系数的反演误差分别仅为2%和4%,较TSVD算法有较大改善.     

融合降水实况分析产品在四川地区的适用性评估

利用四川省地面自动站2018年6月—2019年5月的逐小时降水观测资料,在邻近插值和双线性插值对比分析的基础上,从晴雨准确率、降水时空特征、降水分量级检验等多个方面,对国家气象信息中心研制的融合降水实况分析产品在四川地区的适用性进行评估分析.评估结果表明:(1)邻近插值和双线性插值对评估结果影响小.(2)融合降水实况分...     

基于加速度传感器数据融合的隧道施工地面沉降监测技术研究

地面沉降的监测和控制一直是制约隧道施工安全的关键问题,目前地表沉降量一般利用水准仪等传统设备、通过人工定期测量获取,无法进行实时动态观测。随着无线传感器网络和微机电系统的高速发展,通过在地表布置无线传感器节点,即可对隧道施工引起的地面沉降进行实时监测。本文对处于盾构施工期的某隧道开展了基于加速度传感器数据融合的地面沉降监测技术研究。首先,在隧道上方和两侧布置了内嵌两套加速度传感器的节点盒,获取原始监测数据;然后,利用旋转矩阵求解其空间法向量,并通过投影得到地表节点在隧道轴向上的倾斜角;最终,通过集中式卡尔曼滤波对上述两类传感器的倾斜角进行融合估计,推导出传感器倾斜角随时间的变化情况。基于该方法,通过布置传感器节点监测网,即可估算隧道施工引起的地表绝对沉降量。     

基于数据规则化和稀疏反演的三维表面多次波压制方法

表面多次波是海洋地震勘探中的主要问题.目前,二维数据驱动的表面多次波压制技术(SRME)已经比较成熟,并且已经成为工业界压制海洋表面多次波的主流方法.但是由于二维SRME算法没有考虑横测线方向上多次波的贡献,导致在处理实际三维海洋资料时存在比较大的误差.将二维SRME算法扩展到三维空间后可以得到三维SRME算法,但是由于目前实际采集的三维海洋资料的观测系统存在拖缆漂移,而且横测线方向采样过于稀疏,直接应用三维SRME算法无法准确预测表面多次波.本文提出的通过数据规则化配合稀疏反演的三维表面多次波压制方法能够解决这种实际资料和三维SRME算法之间的矛盾.本文通过研究数据规则化与反规则化技术,使得数据分布满足三维SRME的要求;通过研究稀疏反演技术,有效解决了横测线方向采样稀疏对于多次波预测的影响,三维实际海洋资料的应用结果验证了方法的有效性和可行性.     

四川盆地深部地壳结构——深地震反射剖面探测

四川盆地位于扬子地块的西北部,被褶皱构造带所围绕,受周缘构造带的侧向挤压作用,盆地卷入了多期次和多边界的构造变形,为开展盆山耦合作用及多边界、多期次构造叠加与复合关系的研究提供了不可多得的理想野外实验室.为揭示四川盆地地壳结构,本文通过对3条不同时间采集的深地震反射剖面数据进行拼接联线处理,获得跨越四川盆地的330 km深地震反射偏移成果剖面,揭示了四川盆地地壳上地幔细结构：沉积层从西北向东南逐渐变薄,在龙门山前沉积层厚度超过15 km,在华蓥山下沉积层减薄到~8 km,且褶皱变形形成华蓥山薄皮褶皱冲断带;莫霍面出现在13~15 s（双程走时）,埋深约40~45 km;并发现从下地壳延伸至地幔的东南向的倾斜反射,从13 s向下延伸至18 s,结合四川盆地及其周边地区其他地球物理和地球化学花岗岩同位素年龄等资料,我们认为这些倾斜反射层是扬子克拉通地台西北缘发生的新元古代俯冲的遗迹.

     

一种数字三分量地震微测井仪器设计

针对当前地震微测井仪器在施工和数据采集方面的难题,从提高数据采集质量及仪器功能性方面提出了几种改进方法.主要利用MEMS加速度传感器结合多路并行处理采集方法设计了数字检波器,并集成姿态检测技术直接实现水平分量数据的校正功能,保证系统采集的质量;提出了Ethernet-RS485分布式总线自动控制技术,实现各级联探管的随机组合,同时采用电动方式控制探管操作,提高系统的集成化程度.文中详细分析了利用姿态检测技术实现水平分量数据校正的方法,解决了级联探管随机组合控制和震源激发方位的判定问题.仪器单节探管支持三分量并行采集,有效分辨率达18位,参数可远程配置,同时具备快速显示探管姿态的功能,且支持远程控制,非常适于复杂环境下的工程勘探,整个系统实现全数字化操作且仅需四芯电缆互连,从而缩小了设备的体积和重量,降低施工布线的复杂性.

     

六种不同煤阶煤的品质因子特征

通过6种不同煤阶煤的弹性测试,分析了煤岩变质程度（以镜质组最大反射率表示）与煤岩品质因子之间的关系;在此基础上,进一步分析了煤岩速度与品质因子之间的关系及纵波与横波品质因子的各向异性特征.研究发现：煤岩镜质组最大反射率与纵波、横波品质因子均存在良好的线性相关性;垂直层理方向的波速与品质因子存在良好的线性关系,好于走向与倾向方向;垂直层理方向的纵波与横波品质因子线性相关,相关系数可达93.5%;煤岩的三方向品质因子存在各向异性,且强于速度的各向异性.通过横波品质因子与纵波品质因子间相互换算的理论公式计算结果与实验室实测结果的对比,证明了该理论公式的可靠性,误差小于10%,为多分量转换波的吸收衰减补偿提供了理论依据.     

基于压缩感知的Curvelet域联合迭代地震数据重建

由于野外采集环境的限制,常常无法采集得到完整规则的野外地震数据,为了后续地震处理、解释工作的顺利进行,地震数据重建工作被广泛的研究.自压缩感知理论的提出,相继出现了基于该理论的多种迭代阈值方法,如CRSI方法（Curvelet Recovery by Sparsity-promoting Inversion method）、Bregman迭代阈值算法（the linearized Bregman method）等.CSRI方法利用地震波形在Curvelet的稀疏特性,通过一种基于最速下降的迭代算法在Curvelet变换域恢复出高信噪比地震数据,该迭代算法稳定,收敛,但其收敛速度慢.Bregman迭代阈值法与CRSI最大区别在于每次迭代时把上一次恢复结果中的阈值前所有能量都保留到本次恢复结果中,从而加快了收敛速度,但随着迭代的进行重构数据中噪声干扰越来越严重,导致最终恢复出的数据信噪比低.综合两种经典方法的优缺点,本文构造了一种新的联合迭代算法框架,在每次迭代中将CRSI和Bregman的恢复量加权并同时加回本次迭代结果中,从而加快了迭代初期的收敛速度,又避免了迭代后期噪声干扰的影响.合成数据和实际数据试算结果表明,我们提出的新方法不仅迭代快速收敛稳定,且能得到高信噪比的重建结果.