多时相InSAR技术的廊坊地区地下水体积变化研究

基于时序InSAR函数模型,分别建立了单主影像和多主影像时序InSAR误差模型,在理论上完善了时序InSAR数学模型。利用构建的随机误差模型,模拟试验研究了随机误差模型对时序InSAR待估参数精度的影响。结果表明,与等权模型相比较,加权模型(随机误差模型)估计的参数精度有一定提高;但由于加权条件下的参数估计模型复杂、计算效率低,目前利用等权方法进行时序InSAR的参数估计更简便易行。     

利用HJ-1B模拟数据反演地表温度的两种单通道算法

针对HJ-1B热红外波段特点,采用修正型QK & B算法,反演广州市2013-01-14的地表温度（land surface temperature,LST）。建立偏微分方程得出,当辐射率误差为0.01时,引起的LST误差约为0.6 K,LST误差与大气透过率成反比,与大气透过率误差成正比,0.1的透过率误差引起LST误差约1 K。大气水汽含量w误差与LST误差成线性关系,当大气水汽含量误差为0.1 g/cm²时,引起LST误差约为0.2 K。LST反演误差与近地表气温误差和大气平均作用温度误差均成正比,1 K的近地表气温误差引起LST反演误差约1K。总的来说,LST反演误差与区间比值和大气平均作用温度误差和近地表气温误差相关。用算法反演出来的广州市地表温度与MOD11_L2温度产品具有较强的空间一致性,温度差值曲线呈正态分布,主要集中在-0.9~0.9℃区域,选取广州市6个观测点,得出修正型QK & B算法和实测地温平均值相差约为0.31 K,MOD11_L2与实测地温的温度平均值相差0.65 K,误差均小于1 K。通过对修正型QK & B算法偏微分方程的推导,可对HJ-1B/IRS中的LST反演进行更细致和精确的分析,为其他针对环境卫星热红外波段类似反演LST的算法提供一定的借鉴,也为后续提高LST反演精度提供科学依据。     

正规化矩阵正定时半参数估计量的统计性质

对于包含系统误差和粗差的观测数据,本文将混合Cook距离引入到半参数模型中,实现了粗差的定值定位。首先,构造补偿最小二乘函数,利用泰勒展式,根据均值漂移模型与删失模型的等价性,导出了观测数据中剔除第i个观测数据前后参数分量和非参数分量相应估计值之间的关系式,为粗差的定位分析奠定了基础。其次,将混合Cook距离作为诊断统计量,进行粗差定位分析,得到了参数分量和非参数分量的诊断统计量的简洁计算公式,为了提高粗差定位的准确性,给出了混合Cook距离参数 Q 和C的一些常用形式,通过合理选择相应参数,计算参数分量和非参数分量的距离函数,实现粗差的定位,并将系统误差和粗差从观测数据中区别开来。最后通过模拟算例和实测数据验证了本文方法的正确性。     

利用方差分量估计的地震同震滑动分布反演

提出了基于方差分量估计的地震同震滑动分布反演方法,该方法不仅可以确定不同数据集的相对权重,而且还能得到平滑因子的值,突破了方差分量估计方法仅用于各类数据定权的实际应用范围。模拟算例验证了此方法的有效性,从解的概率后验密度分布分析来看,方差分量估计方法得到的结果是线性和非线性混合方法反演结果的一个解。用方差分量估计方法反演了Bam地震同震滑动分布,并与线性非线性混合方法反演的结果进行了对比分析。两种方法得到的Bam地震断层滑动分布之间的空间互相关系数为0.999 9,地震同震的最大滑动量相同,都为3.04 m;其深度仅有微小区别,前者为4.50 km,后者为4.47 km;平均滑动量前者为0.714 m,后者为0.718 m。但前者反演计算用时227 s,后者却达到1.5×10⁶ s,约17 d,前者仅为后者的1/6 608。结果表明,本文方法相对于线性和非线性混合方法具有计算简单、计算量小和计算效率高的优点。     

由GOCE引力梯度张量不变量确定卫星重力模型的半解析法

提出了利用半解析方法解算重力场与稳态海洋环流探测器(gravity field and steady-state ocean circulation explorer,GOCE)引力梯度张量不变量观测值来恢复全球卫星重力模型的方法,该方法比最小二乘方法效率高,同样可给出模型的验后方差。推导了半解析法快速解算引力梯度张量不变量的理论公式,给出了该方法利用卫星重力梯度观测数据快速求解重力场模型的基本步骤。利用89.5°倾角圆形严格重复轨道上的模拟梯度观测值验证了半解析方法用于张量不变量解算的理论严密性,并利用GOCE任务61 d 梯度仪坐标系(gradiometer reference frame,GRF)下无噪声和有噪声的I₂和V_zz模拟观测值进行了模型反演和结果分析。结果表明,利用I₂观测值求解模型的精度略优于仅用V_zz分量解算模型的精度。     

基于DEnKF方法的考虑次网格变异性的MODIS雪盖同化

基于通用陆面模型(CoLM)和确定性集合卡尔曼滤波算法发展了一个考虑模型次网格变异性的MODIS雪盖同化方案,提高雪深模拟的估计精度。利用北疆阿勒泰地区5个气象站点2007年11月至2008年4月逐日雪深观测数据对同化结果进行了验证。结果表明,该同化方案不需要对MODIS雪盖观测数据进行扰动,能明显提高雪深模拟的精度。另外,雪深同化结果与地面观测雪深具有一致的时间变化趋势,能准确地反映积雪深度在各个不同时段的变化特性。     

庙西南凸起馆陶组油气差异成藏规律

针对庙西南凸起馆陶组主要成藏砂体类型与油气富集程度的差异,分析油气运聚条件及其对油气成藏的影响,总结差异成藏规律。研究表明,油田区为多断阶“脊—断”接力式油气运聚模式,油气先沿馆陶组底部区域砂岩输导脊向高部位断阶区横向运移,再由“通脊”断裂将油气分配至浅层聚集成藏。增压应力的分布与强度是影响油气运移与保存的关键,控制主要成藏砂体类型与油气富集程度。增压区利于油气保存而不利于油气运移,以厚砂体成藏为主;弱增压区、亚增压区利于油气运移但不利于油气保存,主要为薄砂体成藏。构造低部位弱增压与高部位增压型断块运移与保存条件最优,油气最为富集。同一断块内部,“脊—断”条件控制不同井区的油气丰度,“通脊”断裂切至输导脊的构造位置越高、在输导脊上累计错动面积越大,该区域的油气丰度也越高。这种油气差异成藏规律有效指导了庙西南凸起馆陶组的高效评价,同时为类似地区的油气勘探提供了重要参考价值。     

南黄海中部浅地层剖面数据处理新进展

1原始资料分析南黄海中部主要包括南黄海中部平原、黄海槽浅谷、济州岛西北部沙脊地形、苏北岸外舌状地形体系、侵蚀堆积斜坡等地形单元,地形结构较为复杂,同时由于采集环境和采集设备等因素的影响,导致采集的浅地层剖面数据信噪比较低,多次波干扰严重,多次波的存在影响地震成像的真实性和可靠性,干扰地震资料的解释。     

印尼苏拉威西地区波尼盆地形成及其石油地质特征

波尼盆地是一个新生代裂谷盆地,勘探程度极低,其含油气系统及主要勘探层系不明。采用类比分析方法,通过对周边构造演化和沉积特征相似盆地的烃源岩、储层特征及含油气系统分析,进一步推测波尼盆地发育始新统潮坪-潟湖相烃源岩,品质较好,为Ⅱ—Ⅲ型干酪根,普遍处于成熟—高成熟阶段,以生油为主。该盆地共发育3套潜在储盖组合,最主要为中、上中新统—上新统Tacipi组生物礁储层及上新统区域性海相泥岩储盖组合,油气勘探前景较好。     

时间域精细速度建模的策略与试验

速度在地震勘探中十分重要,压制多次波、偏移成像和地震反演都需要准确的速度模型。为了高效和准确的获取速度模型,基于时间域精细速度建模方法,利用疑似含气区的二维多道地震资料,针对初始模型、层位约束和基于人工速度拾取的速度反演3个方面开展了试验,探讨实际资料中时间域精细速度建模方法的关键点和效果。试验表明：无需人工速度拾取,时间域精细速度建模能够实现浅地层速度的准确反演,进而高效地完成浅地层速度异常的识别;层位约束能提高速度建模的效率和精度;在地层结构复杂区域,人工速度拾取有利于时间域精细速度建模。