如何观测和拍摄日环食

今年5月21日上午将发生日环食,在我国广东和福建等地可目睹到。如果错失这次良机,就要等到2020年6月21日在我国才能再次看到。这次日环食,最大食分为0．9439。环食带从我国北部湾开始,经过我国东南部、东海、日本南部、太平洋、美国西部,在美国得克萨斯州结束。日环食从5月21日清晨6时9分开始,至9时37结束,历时3小时28分。日环食持续时间最长为5分42秒。太阳最大高度为61度。     

宇宙的命运与Ia型超新星（中）

宇宙膨胀的发现 1923年,俄罗斯宇宙学家亚历山大·弗里德曼首先用∧=0的爱因斯坦场方程建立了宇宙学膨胀模型（∧为宇宙学常数）,他发表论文阐述了膨胀宇宙的思想,即曲率分别为正、负、零时的三种情况,称为弗里德曼宇宙模型。     

基于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。佳木斯—萝北段以右旋走滑为主、拉张为辅;五常—佳木斯段以拉张为主、右旋走滑为辅。     

汶川地震频次场异常特征分析

地震频次场是描述地震发生频次时空特征的一种数学方法。将2008年5月12日汶川8.0级地震震中附近区域（30.0°~33.0° N、101.5°~105.5° E）作为研究对象,以自然正交函数展开方法分析频次场典型场时间因子的时间变化特征。当取前8个特征值对应的典型场时,拟合精度可达0.936 8;其中6个典型场显示有异常变化,占总场比重的0.691 1;异常出现时间最早在2004年9月,即震前3年7个月,最迟在震前1个多月,表现出短临异常特征。研究结果表明利用地震频次场方法能够较为理想地提取汶川8.0级地震震前异常。     

时变重力场球面模型反演算法和模拟实验

时变重力场是研究地球内部介质物性变化的重要手段。本文提出了一种适用于地面流动重力测量获得的时变重力信号的场源反演方法,该方法采用球坐标系下的六面体单元来模拟场源介质,适合大尺度地震流动重力测量数据的等效源模型构建。通过引入重力时变信号的一阶光滑先验条件,压制了时变重力信号中的短周期高频分量,可用于提取与地震孕育相关的长周期信号。通过理论和模型实验证明了本文算法的可靠性和稳定性,并使用南北地震带南段2014—2017年的流动重力实测数据进行了反演解释,获得了地壳内部等效场源的视密度时变信号,变化量级在正常地壳密度的±0.7‰之间,其空间形态受川滇菱形块体边界控制。研究成果可用于时变重力场模型解释和深部场源特征提取,可为地震重力前兆信号分析和相关研究提供完备的方法保障。