基于轨道扰动引力谱分析的方法确定低低观测重力卫星反演重力场的空间分辨率

基于低低卫-卫跟踪重力卫星的轨道特性,从垂直和水平两个方向计算了重力卫星高空扰动引力,并根据其谱特性及星载加速度的测量噪声水平分析了重力卫星能反演重力场的阶数.利用EGM96重力场模型分别计算了400 km、450 km和500 km 轨道高度处重力卫星受到的扰动引力谱及扰动引力谱的平均量级,分析其垂直特性表明:在三个轨道高度处能分别能反演150、140和130阶的重力场模型.利用两颗同轨重力卫星相距220 km的特性,计算了400 km、450 km和500 km 轨道高度处纬度相差2°的两颗卫星纬向扰动引力差,即扰动引力水平分量,分析其谱特性,表明:重力卫星能反演至117阶的地球重力场模型.     

模拟分析低低跟踪模式重力卫星反演地球重力场的精度

本文利用卫星重力反演与模拟软件ANGELS系统(ANalyst of Gravity Estimation with Low-orbit Satellites)对低低跟踪模式的重力卫星的关键载荷精度指标进行了深入分析.模拟结果表明:(1)对短弧长积分法而言,在低低跟踪模式的关键载荷精度指标中,重力场反演精度对星间距离变率精度最为敏感;(2)通过对目前在轨运行GRACE的载荷指标进行分析,发现轨道数据的误差主要影响重力场的低阶部分(约小于25阶),较高阶次部分(约大于26阶)主要受星间距离变率的误差限制;(3)如果下一代低低跟踪模式的重力卫星的目标之一是把重力异常反演精度较GRACE提高约10倍,则在保持轨道高度和GRACE相同的前提下,轨道、星间距离变率和星载加速度计等关键载荷指标需要达到的最低精度分别约为2cm、10nm·s-1和3.0×10-10 m·s-2;(4)轨道精度和混频误差将是影响下一代低低跟踪模式重力卫星重力场恢复能力进一步提高的主要制约因素,距离变率精度和加速度计精度存在盈余.     

基于下一代三向车轮双星编队改善地球重力场空间分辨率研究

由于当前GRACE(Gravity Recovery and Climate Experiment)串行式编队存在"南北向条带误差"等缺陷,因此本文基于星间速度插值法开展了利用下一代三向车轮双星编队ACR(Along-Cross-Radial)-Cartwheel提高地球重力场空间分辨率的可行性研究论证.第一,采用GRACE卫星轨道参数和关键载荷精度,利用三向车轮双星编队ACR-Cartwheel-A/B反演了120阶地球重力场.结果表明:基于ACR-Cartwheel-A/B双星编队反演地球重力场的模拟精度较德国波茨坦地学研究中心(GFZ)公布的EIGEN-GRACE02S地球重力场模型的实测精度平均提高2.6倍,从而检验了基于下一代三向车轮双星编队ACR-Cartwheel-A/B反演地球重力场精度优于当前GRACE串行式双星编队的可行性.第二,通过星间速度插值法,采用卫星轨道参数(初始轨道高度350km、平均星间距离100km、初始轨道倾角89°、初始轨道离心率0.0046)、卫星关键载荷精度指标(星间速度10-7 m·s-1、轨道位置10-3 m、轨道速度10-6 m·s-1、非保守力10-11 m·s-2)、观测时间30天和采样间隔10s,基于经向车轮双星编队Lo-AR(Longitudinal-Along-Radial)-Cartwheel-A/B、纬向车轮双星编队La-AR(Latitudinal-Along-Radial)-Cartwheel-A/B和三向车轮双星编队ACR-Cartwheel-A/B,分别反演了120阶地球重力场;在120阶处,累计大地水准面精度分别为5.115×10-4 m、4.923×10-4 m和3.488×10-4 m.结果表明:(1)由于La-AR-Cartwheel-A/B编队的轨道稳定性优于Lo-AR-Cartwheel-A/B编队,因此基于La-AR-Cartwheel-A/B编队反演重力场精度高于Lo-AR-CartwheelA/B编队;(2)由于ACR-Cartwheel-A/B编队可以同时获得轨向、垂向和径向的重力场信息,卫星观测数据具有各向同性优点,因此ACR-Cartwheel-A/B编队是建立下一代高精度和高空间分辨地球重力场模型的优化选择.     

一种新型重力测量卫星系统确定全球重力场的性能分析

本文设计了一种高-低卫星跟踪卫星、低-低卫星跟踪卫星和卫星重力梯度测量相结合的新型重力测量卫星系统,其可在一定程度上发挥卫星重力梯度和低低卫星跟踪卫星两种测量模式各自的优势.基于重力卫星系统指标设计的半解析法,深入分析了不同重力测量卫星系统配置和不同观测量及其不同白噪声水平情况下,新型卫星重力测量模式反演重力场模型的能力.数值模拟分析结果表明：在观测值精度和星间距离相同的条件下,轨道高度是影响重力场反演精度的关键因素;随着星间距离的增大,高频重力场信号反演精度会先提高后降低,轨道高度在200~350 km之间时,星间距离在150~180 km之间时反演精度最优;星间距离变率和卫星重力梯度两类观测值仅在某些精度配置时可达到优势互补,如果某一类观测值精度很高,则另一类观测值在联合解算时贡献非常小或者没有贡献.在300 km轨道高度,若以GRACE和GOCE任务的设计指标1 μm·s^-1/√Hz和5 mE/√Hz来配置新型重力测量卫星系统中星间距离变率和引力梯度观测值的精度,联合两类观测值解算200阶次模型大地水准面的精度比独立解算分别提高1.2倍和2.8倍.如果以实现100 km空间分辨率1~2 cm精度大地水准面为科学目标,考虑卫星在轨寿命,建议轨道高度选择300 km,星间距离变率和卫星重力梯度的精度分别为0.1 μm·s^-1/√Hz和1 mE/√Hz.本文的研究成果可为中国研制自主的重力测量卫星系统提供参考依据.     

基于半解析法有效和快速估计GRACE全球重力场的精度

首先基于半解析法建立了新的GRACE卫星K波段测量系统星间测速、GPS接收机轨道位置和加速度计非保守力误差联合影响累计大地水准面的误差模型;其次,基于各关键载荷精度指标的匹配关系,论证了误差模型的可靠性;最后,基于美国喷气动力实验室（JPL）公布的2006年的GRACE Level 1B实测误差数据,有效和快速地估计了120阶全球重力场的精度,在120阶处累计大地水准面的精度为18.368 cm,其结果和德国地学研究中心（GFZ）公布的EIGEN-GRACE02S全球重力场模型符合较好. 本文的研究为将来国际卫星重力测量计划（如GRACE Follow-On, 360阶）中高阶全球重力场模型精度的有效和快速估计提供了理论基础和计算保证.     

卫星重力梯度恢复地球重力场的频谱分析

根据仪器功率谱密度和重力位系数阶方差的定义,本文建立了卫星重力梯度测量噪声功率谱密度与重力场模型的误差阶方差的直接对应关系,并基于此讨论了重力梯度测量精度、卫星轨道高度以及运行时间对地球重力场恢复精度的影响.相比于传统的基于最小二乘法评估卫星载荷噪声对地球重力场恢复精度的影响而言,本文提出的方法简单、直接,有助于快速设计和确定卫星重力测量计划的有关参数.     

Spectral analysis for recovering the Earth's gravity potential by satellite gravity gradients

根据仪器功率谱密度和重力位系数阶方差的定义,本文建立了卫星重力梯度测量噪声功率谱密度与重力场模型的误差阶方差的直接对应关系,并基于此讨论了重力梯度测量精度、卫星轨道高度以及运行时间对地球重力场恢复精度的影响.相比于传统的基于最小二乘法评估卫星载荷噪声对地球重力场恢复精度的影响而言,本文提出的方法简单、直接,有助于快速设计和确定卫星重力测量计划的有关参数.     

利用解析法有效快速估计将来GRACE Follow-On地球重力场的精度

本文首次利用解析法有效快速估计了将来GRACE（Gravity Recovery and Climate Experiment） Follow-On地球重力场的精度. 第一,基于功率谱原理分别建立了新的GRACE Follow-On卫星激光干涉星间测量系统星间速度、GPS接收机轨道位置和轨道速度以及加速度计非保守力误差影响累计大地水准面的单独和联合解析误差模型. 第二,利用提出的GRACE卫星关键载荷匹配精度指标和美国喷气推进实验室（JPL）公布的GRACE Level 1B实测精度指标的一致性,以及估计的GRACE累计大地水准面精度和德国波兹坦地学研究中心（GFZ）公布的EIGEN-GRACE02S地球重力场模型实测精度的符合性,验证了本文建立的解析误差模型是可靠的. 第三,论证了GRACE Follow-On卫星不同关键载荷匹配精度指标和轨道高度对地球重力场精度的影响. 在360阶处,利用轨道高度250 km、星间距离50 km、星间速度误差1×10^-9m/s、轨道位置误差3×10^-5m、轨道速度误差3×10^-8m/s和非保守力误差3×10^-13m/s²,基于联合解析误差模型估计累计大地水准面的精度为1.231×10^-1 m. 本文的研究不仅为当前GRACE和将来GRACE Follow-On地球重力场精度的有效快速确定提供了理论基础和计算保证,同时对国际将来GRAIL（Gravity Recovery and Interior Laboratory）月球卫星重力测量计划的成功实施具有重要的参考意义.     

基于引力位系数相对权重的卫星重力场测量分析

卫星重力场测量已成为最有效的全球重力场测量手段.本文结合典型的重力卫星和重力卫星研究计划,分析了卫星重力测量的三种原理,并基于各阶位系数的相对权重讨论了各种原理的应用优势.分析可知,卫星受摄轨道适用于恢复长波重力场,低轨星间距离变化率适用于恢复中长波重力场,重力梯度适用于恢复中短波重力场.针对中长波高精度重力场测量的需要,设计了综合获取低轨星间距离变化率与受摄轨道的重力卫星方案,该方案由两组内编队组成星星跟踪复合编队,轨道高度为250km,星间距离为50～100km.     

基于下一代四星转轮式编队系统精确和快速反演FSCF地球重力场

第一,由于重力卫星编队轨道的稳定性设计是建立下一代高精度和高空间分辨率地球重力场模型的关键,因此为保证下一代四星转轮式编队系统的稳定性,轨道根数的最优设计如下:(1)轨道半长轴a、轨道偏心率e、轨道倾角i和升交点赤经Ω保持不变;(2)每对卫星的近地点幅角ω和平近点角M分别相差180°;(3)初始近地点辐角ω设置于赤道处,初始平近点角M设计于极点处;(4)卫星编队系统椭圆轨道的半长轴和半短轴之比为2:1. 第二,基于下一代四星转轮式编队系统,利用星间速度插值法,通过相关系数(激光干涉测量系统的星间速度0.85、GPS接收机的轨道位置和轨道速度0.95、星载加速度计的非保守力0.90)、观测时间30天和采样间隔10 s,反演了120阶FSCF-1/2/3/4(Four-Satellite Cartwheel Formation)地球重力场,在120阶处累计大地水准面精度为1.162×10^-4 m,较目前GRACE地球重力场精度至少提高一个数量级. 第三,下一代四星转轮式编队系统具有低轨道高度、高精度测量、全张量观测、弱混频效应和强时变信号的优点.     

测井数据岩石物理参数敏感性的应用研究

利用测井资料所获得的基本参数进行岩石物理参数敏感性分析,可以较大提高含烃岩石的识别能力。本文讨论了对地震属性敏感性高的参数的选择及其使用方法,以便提高流体和岩性的识别能力。为了更好的与地震数据解释相结合,我们在敏感性参数空间构造了研究区的岩石物理模板,用以指导实际地震数据定量解释。     

各向异性TI介质qP反射波走时层析成像

地震走时层析成像是反演地层各向异性参数分布的有效方法,但是关于地震各向异性介质走时层析成像的研究并不多,其技术远远没有达到成熟的阶段.在野外数据采集时,地表反射波观测方式相对井间和垂直地震剖面观测方式的成本更低,利用qP反射波走时反演各向异性参数具有更加广泛的实用价值.本文实现的TI介质地震走时层析成像方法结合了TI介质反射波射线追踪算法、走时扰动方程和非线性共轭梯度算法,它可以对任意强度的TI介质模型进行反演,文中尝试利用qP反射波走时重建TI介质模型的参数图像.利用qP反射波对层状介质模型和块状异常体模型进行走时反演,由于qP波相速度对弹性模量参数和Thomsen参数的偏微分不同,所以可以分别反演弹性模量参数和Thomsen参数.数值模拟结果表明:利用qP反射波可以反演出TI介质模型的弹性模量参数与Thomsen参数,不同模型的走时迭代反演达到了较好的收敛效果,与各向同性介质走时反演结果相比较,各向异性介质走时反演结果具有较好的识别能力.     

上海及邻近地区介质参数的研究

在频率域里研究了上海及邻近地区的Q值、小震震源参数,初步得到如下结果:①以上海佘山为中心,平均直达P波Q值呈方向性分布。在相同震中距下,上海西面溧阳方向(?)值较高,上海北部南黄海方向(?)值较低;溧阳和南黄海地震活跃区域的(?)值比其邻近地震相对不活跃区域的偏高,②相同震级的地震相比之下,上海—溧阳方向的地震破裂尺度小,上海—南黄海方向的地震破裂尺度大;并且震级在2.5—4.0级之间的地震,破裂尺度的数值溧阳的比较集中,南黄海的则比较离散。     

spartina in Mariager Fjord,Denmark: The effect on sediment parameters

The analysis of 3 subsamples from each of 12 sediment cores from a micro-tidal flat shows that spreading of Spartina effects sediment parameters. The effect diminishes with distance to vegetation.     

储层弹性与物性参数地震叠前同步反演的确定性优化方法

储层弹性与物性参数可直接应用于储层岩性预测和流体识别,是储层综合评价和油气藏精细描述的基本要素之一.现有的储层弹性与物性参数地震同步反演方法大都基于Gassmann方程,使用地震叠前数据,通过随机优化方法反演储层弹性与物性参数;或基于Wyllie方程,使用地震叠后数据,通过确定性优化方法反演储层弹性与物性参数.本文提出一种基于Gassmann方程、通过确定性优化方法开展储层弹性和物性参数地震叠前反演的方法,该方法利用Gassmann方程建立储层物性参数与叠前地震观测数据之间的联系,在贝叶斯反演框架下以储层弹性与物性参数的联合后验概率为目标函数,通过将目标函数的梯度用泰勒公式展开得到储层弹性与物性参数联合的方程组,其中储层弹性参数对物性参数的梯度用差分形式表示,最后通过共轭梯度算法迭代求解得到储层弹性与物性参数的最优解.理论试算与实际资料反演结果证明了方法的可行性.     

基于阶跃响应在时域解析测定地震计参数

提出基于阶跃响应在时域解析测定地震计参数的方法，并对合理确定响应的时间零点提供了一种方法。对JC-V100 3D型便携式短周期地震计及与其配套的EDAS-3M型16位地震数据采集记录仪组成的系统进行了测试。结果表明：对JC-V1003D短周期地震计的传递函数进行测试是可行的，对地震计的阻尼系数Dc和自振周期Tc的测试结果是令人满意的。     

2021年12月中国大陆地区M≥4.0地震震源机制解测定

为了填补测震台网地震震源参数目录常规数据产出方面的空白,本文利用中国地震台网记录的宽频带波形资料,采用近震全波形反演方法对2021年12月1—31日发生在中国大陆地区的M≥4.0共6次地震进行震源机制反演。结果显示逆断型3次,走滑型2次,未知型1次。      

2022年1月中国大陆地区MS≥4.0地震震源机制解测定

本文利用中国地震台网记录的宽频带波形资料,采用近震全波形反演方法得到2022年1月1—31日发生在中国大陆地区的M_S≥4.0共13次地震的震源机制解。结果显示逆断型1次,走滑型9次,未知型3次。      

渤海地区地震参数的修定

本文介绍了一种横向不均匀介质中地震参数的修定方法。该方法以渤海及其邻区地壳、上地幔结构的最新成果为基础构造速度模型,集中起可利用的渤海周边地区的山东、河北,北京和辽宁等台网的资料,统一修定了近年来渤海及其邻区598个地震的震源参数。所有地震都给出了震源深度并且讨论了地震定位的精度。文中给出了近年来渤海内M_L≥4.0级的地震震源参数修定结果,讨论了该区震源深度分布的特征。     

浙江珊溪水库地震震源参数研究

利用浙江省区域数字地震台网记录的珊溪水库地震波资料,采用Brune(1970)模式,通过几何扩散校正、介质衰减校正、仪器校正等将速度记录谱归算为震源位移谱,使用遗传算法计算拐角频率及零频极限,进而计算震源半径、地震矩等小震震源参数,并系统分析了各种参数之间的关系。研究中使用发生于2002～2006年间ML2.0～4.6珊溪水库库区地震180个,分析结果为:地震矩范围为1011～1015N.m,震源破裂半径为74～727m,地震矩、拐角频率、矩震级及体波震级之间表现出一定的对数线性或半对数线性关系,并可分别表示为:lgM0=0.0061R 11.71,lgM0=1.11ML 10.11,MW=0.7447ML 0.7。地震应力降范围为0.08～4.95MPa,应力降与震源半径之间表现出一定的相关关系。