堆石防波堤不规则波浪反射系数试验研究

结合物理模型试验,分析斜坡坡度、波陡、相对水深、护面类型和破波参数等因素对堆石防波堤不规则波浪反射系数的影响规律。将常用的Van der Meer公式,Seelig公式,Postma公式和Davison公式计算值和实测值进行比较,并结合试验数据,基于有效波高和平均周期定义的Iribarren数,得出堆石防波堤不规则波浪反射系数经验公式。结果表明,该公式能较好地计算不规则波作用下块石和扭王块体护面堆石防波堤波浪反射系数。     

高性能GIS研究进展及评述

互联网技术的发展使地理信息技术得到了前所未有的发展和应用,地理信息计算呈现出计算速度快、运行效率高、应用多样化的发展特征。而随着计算机硬件性能飞速提升,传统的GIS数据处理方式并不能与之匹配,各种缺陷与弊端逐渐显现,亟待更高效的数据处理方式。目前,以并行集群计算技术和分布式网络技术为代表的高性能计算的出现,为这些问题的解决带来了新思路,并逐渐发展形成了新一代的多核并行高性能计算系统。当前,如何利用新型硬件体系结构带来的计算能力,研究新一代高性能GIS计算系统,解决现在所面临的时空数据密集和计算密集问题成为重要挑战。高性能计算是基于一组或几组计算机系统组成的集群,通过网络连接组成超级计算系统以加强数据处理、分析计算性能的一种技术。在实际应用中,逐渐形成Hadoop,Spark和Storm 3大主流分布式高性能计算系统,它们三者各具优缺点。本文从高性能GIS算法、并行GIS计算、内存计算和众核计算4个方面梳理、归纳总结了高性能GIS的技术体系,分析了每类高性能GIS技术特征,综合分析、评述了近年来高性能GIS的研究进展,并对高性能GIS未来发展进行展望,为更完备、高效的高性能GIS体系的建立、发展和应用提供参考。今后,并行GIS计算、高性能计算模式和分布式存储仍然是GIS技术领域发展的重要方向,通过高性能GIS系统可有效地解决时空数据密集、计算密集和网络通讯密集等问题,大大提升GIS地理分析效率。     

沙漠地区地震勘探随机噪声建模及其在噪声压制中的应用

随机噪声是影响地震勘探有效信号的主要因素,其存在大大降低了地震记录的信噪比.在噪声压制方法不断被改进的同时,对随机噪声特性进行研究,了解噪声的产生机制是对其进行压制的先决条件,目前对噪声的研究主要是特性研究以寻找规律性,对其进行定性定量的分析还比较少.本文根据塔里木沙漠地区实际采集环境,考虑到噪声的连续性给计算带来的不便,假设各类噪声源以点源的形式分布在检波器周围,依据相应理论确定各类噪声源的源函数,其激发的噪声经由波动方程传播,将随机噪声作为各类噪声源共同作用的综合波场,建立随机噪声的理论模型.通过分析不同种噪声对地震记录的影响,选取合适的滤波方法对其进行压制,实验结果表明,通过建立沙漠地区随机噪声的理论模型,为选择有效的滤波方法,提高地震记录信噪比起到理论指导作用.     

The utilization of the double focal transformation for sparse data representation and data reconstruction

In many cases, seismic measurements are coarsely sampled in at least one dimension. This leads to aliasing artefacts and therefore to problems in the subsequent processing steps. To avoid this, seismic data reconstruction can be applied in advance. The success and reliability of reconstruction methods are dependent on the assumptions they make on the data. In many cases, wavefields are assumed to (locally) have a linear space–time behaviour. However, field data are usually complex, with strongly curved events. Therefore, in this paper, we propose the double focal transformation as an efficient way for complex data reconstruction. Hereby, wavefield propagation is formulated as a transformation, where one‐way propagation operators are used as its basis functions. These wavefield operators can be based on a macro velocity model, which allows our method to use prior information in order to make the data decomposition more effective. The basic principle of the double focal transformation is to focus seismic energy along source and receiver coordinates simultaneously. The seismic data are represented by a number of localized events in the focal domain, whereas aliasing noise spreads out. By imposing a sparse solution in the focal domain, aliasing noise is suppressed, and data reconstruction beyond aliasing is achieved. To facilitate the process, only a few effective depth levels need to be included, preferably along the major boundaries in the data, from which the propagation operators can be calculated. Results on 2D and 3D synthetic data illustrate the method's virtues. Furthermore, seismic data reconstruction on a 2D field dataset with gaps and aliased source spacing demonstrates the strength of the double focal transformation, particularly for near‐offset reflections with strong curvature and for diffractions.     

卫星和雷达资料在暴雨数值预报中的应用研究

采用美国风暴分析预报中心开发的中尺度模式（Advanced Regional Prediction System,ARPS）的资料分析系统ADAS（ARPS Data Analysis System）,把中国CINRAD/SA、SB多普勒雷达反射率、风云-2E卫星的红外和可见光资料引入模式初始场,利用高分辨率中尺度模式WRF对2012年7月21日中国华北地区暴雨过程进行数值模拟。结果表明：加入雷达和卫星资料的WRF数值模拟,由于修正了中低层水汽,对降水系统预报有一定的调整,预报的降水落区位置和降水量中心值更接近实况;同时加入卫星和雷达资料同化后,明显改善了WRF降水系统预报偏慢的情况。雷达与卫星资料同化对WRF模式降水预报影响不同,雷达资料影响的时间短,只在6 h内预报起作用,而卫星资料影响时间较长,能够改善24 h内的预报。加入卫星和雷达资料,修正了水汽场,使降水预报在开始时间、量级和范围上均有改善,更接近实况。     

不同围压和流体饱和状态下致密砂岩弹性各向异性特征

为研究致密砂岩声波速度及其各向异性随围压的变化规律以及不同流体饱和状态下的弹性各向异性特征,钻取了不同方向的岩心并在实验室超声波频率下对致密砂岩的声学特性进行了测量,分别给出干燥和饱和水状态下,不同方向样品纵横波速度、刚性系数以及各向异性系数随围压的变化规律,并对实验结果进行了分析讨论.实验结果表明致密砂岩纵横波速度、纵横波速度比以及刚性系数均随围压增加而增加,但其在不同饱和状态下的变化率却截然不同;纵横波速度比、各向异性系数在饱和水状态下变化规律不明显,表明孔隙流体的存在对于岩石物理性质有着非常重要的影响.这方面的实验工作不但对于考察不同流体性质对致密岩石弹性各向异性影响是必要的,而且有助于致密砂岩油水和气层的识别.     

平均入射角道集PP波与PS波联合反演

在界面两侧地层的弹性参数弱反差的假设难以成立的情况下,本文提出用平均入射角道集进行PP波与PS波的联合反演.首先,在PP波与PS波AVA(amplitude versus angle,振幅随入射角变化)道集的基础上,分别选择小入射角范围与大入射角范围的AVA道集进行局部加权叠加,以获得由两个角度组成的平均入射角道集,并作为后续反演的输入数据.然后,再通过最小二乘原理建立了PP波与PS波联合反演目标函数,推导了模型修改量的向量公式,建立了平均入射角道集联合反演的流程.模型数据与实际数据的测试结果表明:在信噪比较低、地层弹性参数反差较大、层厚较薄的情况下,该反演方法的精度在很大程度上超过了基于近似反射系数的反演方法,为复杂油气藏勘探提供了新的思路.     

Comparison of P‐wave attenuation models of wave‐induced flow

Wave‐induced oscillatory fluid flow in the vicinity of inclusions embedded in porous rocks is one of the main causes for P‐wave dispersion and attenuation at seismic frequencies. Hence, the P‐wave velocity depends on wave frequency, porosity, saturation, and other rock parameters. Several analytical models quantify this wave‐induced flow attenuation and result in characteristic velocity–saturation relations. Here, we compare some of these models by analyzing their low‐ and high‐frequency asymptotic behaviours and by applying them to measured velocity–saturation relations. Specifically, the Biot–Rayleigh model considering spherical inclusions embedded in an isotropic rock matrix is compared with White's and Johnson's models of patchy saturation. The modeling of laboratory data for tight sandstone and limestone indicates that, by selecting appropriate inclusion size, the Biot‐Rayleigh predictions are close to the measured values, particularly for intermediate and high water saturations.     

three-dimensional sensitivity coefficients of apparent resistivity and preliminary application

We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.     

小台站的大数据和云计算设想

地球科学是一门观测科学,相关科学研究及公益服务,如气象、地震等,都依赖于野外台站的观测。目前,由于多种因素的作用,导致台站观测和后续分析应用严重脱节。本文基于大数据和云计算的思想,将台站工作和目前后续研究中心的工作进行有机结合,提出了一套新的工作思路,即将数据计算分析和区域监测预报的工作由台站来完成,然后由一个综合的数据中心负责存储和管理台站提供的"大数据"信息,并根据台站需要和反馈情况研发相关软件,为台站的分析计算提供云计算的平台。