塔里木盆地中原探区速度场研究方法

利用自行编程设计的速度场算法,建立了塔里木盆地中原探区的速度场模型,并应用于油气勘探生产之中．在塔中地区,影响平均速度的主要因素是地层埋深和火成岩的岩性．在巴楚地区,影响平均速度变化的主要因素是大的构造运动形成的多角度不整合界面和断距较大的断层．总结出了一套适合地质复杂地区速度场建立的方法;利用该速度场进行变速成图,能够提高成图精度、排除假构造、发现新的圈闭、减小勘探风险,落实发现了一批圈闭,提供探井多口,实现了油气勘探的突破。     

Structural and depositional controls on Plio‐Pleistocene submarine channel geometry (Taranaki Basin,New Zealand)

High‐quality 3D seismic data are used to investigate the effect of the Parihaka Fault on the geometry of submarine channels in Northern Graben of the Taranaki Basin, New Zealand. The Parihaka Fault comprises of four segments (S1–S4) with variable displacements. As part of the Plio‐Pleistocene Giant Foresets Formation, the older Channel Complex Systems 1 and 2 reveal a two‐stage evolution: (a) a syn‐tectonic depositional stage with channels incising the slope during early fault growth (ca. 4.5 Ma) and (b) a stage of sediment bypass (ca. 3 Ma) leading to the infill of hanging‐wall depocentres. The Channel Complex System 3 is syn‐tectonic relative to segment S3 and was formed at ca. 2.5 Ma. We show that the successive generation of new fault segments towards the north controlled the formation of depocentres in the study area. This occurred in association to rotation and uplift of the footwall block of the Parihaka Fault and subsidence of its hanging‐wall block, with fault activity controlling the orientation of channel systems. As a result, we observe three drainage types in the study area: oblique, transverse and parallel to the Parihaka Fault. This work is important as it shows that relay zones separating the Parihaka Fault segments had limited influence on the geometry and location of channel systems. Submarine channels were diverted from their original courses close to the Parihaka Fault and flowed transversally to fault segments instead of running through relay ramps, contrasting to what is often recorded in the literature. A plausible explanation for such a discrepancy relates to rapid progradation of the Giant Foresets Formation during the Plio‐Pleistocene, with channel complexes becoming less confined, favouring footwall incision and basinward deposition of submarine fans.     

利用H-K叠加方法和CCP叠加方法研究中国东北地区地壳结构与泊松比

利用中国东北布设的流动地震台阵（116个）以及国家和区域台网（121个）的宽频带台站记录的824个远震事件,采用P波接收函数CCP叠加和H-K叠加两种方法获得了研究区详尽的地壳厚度图像.研究结果显示,两种方法获得的地壳厚度分布特征具有很好的一致性,中国东北下方地壳厚度存在明显的东西横向差异.重力梯度带西侧和佳木斯地块的台站下方地壳较厚,介于36~41 km之间,而在兴蒙槽地褶带中重力梯度带往东从36 km减薄至34 km左右.松辽盆地北侧、东侧和南侧地壳厚度较薄,为29~34 km,反映了该区复杂的地壳变形过程.CCP剖面显示郯庐断裂深切地壳,敦化—密山断裂下方莫霍面出现错断.H-K叠加得到的地壳平均泊松比显示,东北地区绝大部分台站下方的泊松比值较大,0.24~0.29.长白山、松辽盆地东部、燕山台隆东部和大兴安岭北部,泊松比值达到0.27~0.30,可能有幔源物质上涌,下地壳铁镁组分含量增加.     

基于Chebyshev自褶积组合窗的有限差分算子优化方法

有限差分法广泛应用于地震波数值模拟、成像和波形反演中,差分数值解的精度直接影响着地震成像和反演的效果.因为有限差分算子可以通过截断伪谱法的空间褶积序列得到,而截断窗函数的属性影响有限差分算子逼近微分算子的精度.具体地讲,窗函数的幅值响应的主瓣和旁瓣决定了有限差分算子逼近的精度,主瓣越窄,旁瓣衰减越大,则有限差分算子逼近微分算子的精度越高,更好地压制数值频散.基于此认识,本文提出了一种基于Chebyshev自褶积组合窗截断逼近的有限差分算子优化方法.Chebyshev自褶积组合窗的主瓣较窄,且旁瓣衰减大,其可通过只调节三个参数,更直观和可视化地控制主瓣和旁瓣的形状,改变有限差分算子逼近微分算子的精度;该窗函数截断逼近的有限差分算子不仅有较大的谱覆盖范围,而且精度误差波动较小,这表明低阶的差分算子可以达到高阶算子的精度,且逼近误差更稳定;从经济上来讲,将有效地减少模拟计算花费,提高计算效率.     

基于压缩感知的高分辨率平面波分解方法研究

平面波分解方法是地震资料处理中的一项关键技术,它广泛地运用在平面波偏移、各类Beam偏移等成像方法中.平面波分解不仅可以提高偏移成像的效率,而且可以压制地震资料中的随机噪音,提高地震数据的信噪比.线性Radon变换（LRT）是一种常见的实现平面波分解的方法,但常规的LRT存在以下两个缺点：（1）分辨率受测不准原理限制;（2）变换结果存在很多噪音和空间假频.为了克服LRT的上述缺点,本文提出一种基于压缩感知的高分辨率平面波分解方法,并利用加权匹配追踪（WMP）技术实现了该方法.该方法将LRT视为一个参数估计问题,并将LRT结果的稀疏性作为约束条件,在压缩感知理论的指导下利用WMP方法得到高分辨率、高信噪比的平面波分解结果.另外,该方法还可以用于提取地震数据的线性信号、压制随机噪音、实现高维地震数据规则化等地震资料处理技术.数值实验结果证明：WMP方法可以有效地提取地震数据中的线性信号,提高LRT的分辨率和信噪比,从而改善平面波分解的质量.     

基于地理加权回归的漫湾库区景观破碎化及影响因子分析

应用地理加权回归模型分析漫湾库区景观破碎化指数——有效筛网大小与相关因子之间的空间关系。选取的解释变量分别是距道路的距离、距乡村的距离、距河流的距离、坡度。结果表明：大坝修建后4种解释变量与有效筛网大小呈现较显著的正相关性。与线性回归模型相比,地理加权回归模型的拟合效果显著提高。1974~1988年,有效筛网大小对各影响因子最敏感的区域面积呈现显著的时空变化这为确定水电站建设及其他因素对景观破碎化影响的大小,并进一步改善库区景观破碎化的现状提供了依据。     

基于导航雷达海面回波图像去噪研究

根据海面X波段雷达回波图像数据的时空特点,提出了一种新的去噪方法。该方法利用线检测对雷达图像中的射线状噪声进行分离,然后针对这种特殊的噪声形式提出一种阈值分割的方法,最后在极坐标下针对噪声点进行插值修复。实验表明,与传统的均值滤波和中值滤波等方法相比,新方法无论在主观视觉还是信噪比方面,都有了一定改进,图像保留了更多的细节信息,具有较好的去噪效果。     

Fast Geographically Weighted Regression (FastGWR): a scalable algorithm to investigate spatial process heterogeneity in millions of observations

Geographically Weighted Regression (GWR) is a widely used tool for exploring spatial heterogeneity of processes over geographic space. GWR computes location-specific parameter estimates, which makes its calibration process computationally intensive. The maximum number of data points that can be handled by current open-source GWR software is approximately 15,000 observations on a standard desktop. In the era of big data, this places a severe limitation on the use of GWR. To overcome this limitation, we propose a highly scalable, open-source FastGWR implementation based on Python and the Message Passing Interface (MPI) that scales to the order of millions of observations. FastGWR optimizes memory usage along with parallelization to boost performance significantly. To illustrate the performance of FastGWR, a hedonic house price model is calibrated on approximately 1.3 million single-family residential properties from a Zillow dataset for the city of Los Angeles, which is the first effort to apply GWR to a dataset of this size. The results show that FastGWR scales linearly as the number of cores within the High-Performance Computing (HPC) environment increases. It also outperforms currently available open-sourced GWR software packages with drastic speed reductions – up to thousands of times faster – on a standard desktop.     

山区煤矿三维地震勘探叠前新抽道集方法

目前，山区地震资料处理，仍采用常规的抽道集后进行水平叠加，但是山区激发点与接收点高差大，且煤层埋深相对较浅，按常的抽道集方法不能实现共反射点叠加，从而影响了解释小陷落柱、小断层的精度。本文提出一种新的抽道集方法，同时解决了因地形高差和目的层起伏引反射点在X、Y方向偏离中心点的归位问题，无是二维还是三维勘均可实现共反射点叠加，从而所得资料可准的反映地下构造，该技术在目的层埋藏浅的开展地震勘探有重要的指导作用和应用价值。