基于CLG光学流和波场分解的逆时偏移角度道集提取方法研究

角度域共成像点道集是衔接叠前地震数据与储层特征的重要桥梁,对地震偏移成像与储层描述具有重要意义.与克希霍夫偏移和单向波动方程偏移相比,逆时偏移是复杂地区最精确的成像方法.高效稳健地生成逆时偏移角度道集目前仍然是一个挑战.本文主要讨论如何采用光学流方法高效、高质量地提取角度道集.在逆时偏移波场外推过程中,光学流方法可以估计波场传播方向.其中Lucas-Kanade(LK)和Horn-Schunck(HS)方法是光学流方法中两种典型的方法.LK光学流方法是一种局部方法,该方法依赖于局部点的梯度值,但是容易出现奇异现象,HS光学流方法属于全局方法,波场方向估计依赖于整个波场,易受噪声影响,对异常值比较敏感,导致整体波场方向计算精度不高.本文提出采用局部和整体结合(Combining Local and Global,CLG)的光学流方法估计波场传播方向.该方法可以有效地提高波场方向的精度,并且简单高效,便于并行处理.对比HS光学流方法,CLG光学流方法几乎不增加额外的计算量.另外,为了弱化光学流方法无法处理波前重叠问题,本文利用解析波场和方向滤波对波场进行方向分解,仅需波场的空间傅里叶变换即可实现任意波场方向分解,将分解后的波场分别估计波场反向,提取成像结果.进一步地,在估计反射张角和方位角时,本文提出有效的归一化方法和改进的最小二乘除法,提高角度估计的精度和稳定性.最后,理论和实际资料例证了本文提出方法的有效性.     

基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.     

浅海域海底极低频电磁响应特征研究

极低频电磁探测方法(WEM)是一种新型人工源低频电磁探测技术,目前尚未见极低频电磁探测方法的海区研究,缺少相关经验.为测试WEM在海洋环境中的电磁响应特征,本文将"地-电离层"全空间电磁波传播理论应用到含海水介质的"地-电离层"全空间模型,由地表处数值格林函数出发,利用空气层与海水层界面谢昆诺夫势函数满足的边界条件,推导了海底界面的场强表达式.采用高采样密度的Hankel数字滤波方法进行数值模拟计算.为探究WEM在海洋油气资源探测中的适用性,建立简化海洋地电模型,通过分析WEM对不含异常体的背景场和含高阻油气异常体时的电磁响应曲线差异来检测WEM探测海底高阻油气层的可行性.并讨论不同观测方式、不同收发距以及海水层厚度、油气层埋深和油气藏厚度对WEM电磁响应特征的影响.研究结果表明,WEM可在数百千米收发距下探测浅海区域海底下方电性异常体的存在.     

瞬变电磁扩散场到虚拟波场的精细积分变换算法

现有的瞬变电磁解释与反演方法主要对电阻率参数进行解释,较难给出准确的地质构造信息.为了能够较多获得地质构造信息,可将瞬变电磁扩散场信号转换为虚拟波场信号,提高瞬变电磁法的分辨率.本文采用精细积分法,实现了瞬变电磁扩散场到虚拟波场的转换.将一个高度病态的线性方程组转换成求积分的过程,大大降低了解决病态问题的难度.积分步长以2的指数增加,在保证精度的前提下,具有极高的计算效率;设置合理的终止迭代条件,进一步提高本文方法的适用性.对典型地电模型的计算,表明本文方法具有较高的精度;减小两个波峰之间的距离,证明了本文方法具有较好的分辨能力;对含有噪声信号的扩散场进行波场反变换,证明了本文方法具有较好的抗噪性.最后计算三维模型和实测数据的虚拟波场,可以很好地反映出三维地质体的界面信息,证明本文方法稳定、可靠,可以对地下地质目标体的界面进行有效识别.本文方法相比预条件正则化共轭梯度法(PRCG),在计算效率方面提升了 4倍左右;本文方法的相对误差在5％以内,而PRCG的最大相对误差可达40％～50％,精度最高可提升10倍;在分辨率与稳定性方面,也表现出较好的效果.     

复波场的实部和虚部对同时源波形反演的影响

高精度的成像需要准确的速度场信息,波形反演被认为是目前具有最高分辨率的速度反演方法之一.计算效率是目前全波形反演需要考虑的一个主要问题.为了很好地解决计算效率的问题,本文引入了一种高效的无串扰同时源反演方法,并详细介绍了其原理与计算流程.同时,基于此多震源同时反演方法,本文拓展出实波场反演、虚波场反演及复波场反演的三种反演策略,进而分析复波场的实部和虚部对全波形反演的影响.实验表明,相比于复波场全波形反演,无论是实波场反演还是虚波场反演,反演分辨率有所降低;相比于实波场反演,虚波场反演对初始模型的依赖性较小,目标函数的非线性较弱;最后,通过使用组合反演策略,即初始阶段采用虚波场反演,中后期阶段采用复波场反演,不仅可以降低反演的非线性,而且能够保证高精度建模.     

依兰—伊通断裂带黑龙江段构造运动特征

为研究依兰—伊通断裂带黑龙江段构造运动特征,基于2016—2019年GPS和地质资料,解算了该断裂的三维速度场,通过构建断层模型反演了滑动速率。结果显示:依兰—伊通断裂带黑龙江段总体呈下沉趋势,沉降速率在1~2 mm/a,断裂呈右旋走滑态势,闭锁层15 km以下走滑速率为(1.7±0.4)mm/a。佳木斯—萝北段以右旋走滑为主、拉张为辅;五常—佳木斯段以拉张为主、右旋走滑为辅。     

三维井间地震侧面波传播特征及其分离研究

井间地震与常规地震相比具有较高的分辨率,同时波场也更加复杂.利用视速度特征的差异很容易地分离井间上下行波场,根据纵横波偏振特征的不同可将纵横波完全分离,但是针对井间地震侧面波的研究并不多.本文构建高速透镜体和断层模型进行波场数值模拟,通过对三维地震波场的传播特征进行分析可知,井间地震接收的波场既有岩性分界面的响应,还会有侧向地质不均匀体产生的干扰,严重影响成像的质量和精度.经偏振旋转后的T分量记录对物性横向变化较为敏感,综合分析X、Y、Z分量和T分量上不同连井剖面的波场传播特征,可对侧向传播的干扰波进行有效识别,利用线性Radon变换可以在一定程度上将侧面波进行压制,为后续井间地震资料处理提供高信噪比的数据.本文通过理论模型进行试算研究,验证了该方法的有效性,并对理论记录取得了较好的分离效果.     

关于地平式反射望远镜的CCD平场(英文)

现有标准设计的反光望远镜在拍摄CCD平场时,绝大多数都受残留散射光的影响.但是赤道式与地平式反光镜所受影响是不同的.在赤道式上做时间序列较差测光时,只要待测星永远位于CCD的固定像元上,不太准确的CCD平场也能得到高精度的测光结果.当需要0.1%～0.3%精度的平场时,则可以采用夜天平场.地平式的特点是,它的CCD相机必须置于旋转器上,在跟踪天体时不停地旋转,以抵消地球自转的影响.上述用于赤道式的方法失效,因此,在CCD平场时,消除散射光的影响比赤道式更为重要.一个典型的地平式反光镜的例子是NAOC兴隆天文台的EOS 1米镜.虽然该台已附加了防散射光的装置,但是对所有B、V、R、I滤光片,在不同旋转器位置拍摄的CCD平场,仍然有2%～3%的差别(主要是梯度).该文给出了改进的建议,必须满足下面两个条件:Cij=C(r);旋转器的中心与反光镜的光学中心重合.此问题的解决对所有地平式反光望远镜都有普遍意义.     

1．2m地平式望远镜视场旋转研究与消旋

简要地讨论了1.2m地平式望远镜的运动特征,定量地给出物方视场旋转的公式,并对像方视场旋转的量和方向给予确定.通过对三种消旋方式的比较得出物理消旋更适合1.2m地平式望远镜视场消旋的结论.     

异地异构天文数据资源的统一访问

天文观测数据资源具有时间跨度大,数据量大,存储、管理分散,管理工具驳杂等特点.如何提供给天文学家一个统一访问这些异地存放的异构数据资源的方案,是虚拟天文台的一个重要研究课题.本文提出了一种基于网格的解决方案.我们基于OGSA-DAI技术实现了对天文星表数据,图像数据和光谱数据的封装(DataNode).我们设计的虚拟天文台数据访问服务(VO-DAS)实现了对DataNode的自动资源发现以及异地异构数据库的联合访问.这将使天体的多波段交叉证认成为可能.VO-DAS支持国际虚拟天文台联盟(IVOA)的各项相关标准,使得它具有良好的互操作性.它的对外接口简单实用,可以针对不同需求的天文数据用户发展出多种网格应用产品.通过在两个科学范例中应用VO-DAS服务,证明了它的可行性.