Assessing model of highway slope stability based on optimized SVM

Considering the geological hazards attributed to the highway slope, using a common simple model cannot accurately assess the stability of the slope. First, principal component analysis (PCA) was conducted to extract the principal components of six factors (namely, bulk density, cohesion, internal friction angle, slope angle, slope height, and pore water pressure ratio) affecting the slope stability. Second, four principal components were adopted as input variables of the support vector machine (SVM) model optimized by genetic algorithm (GA). The output variable was slope stability. Lastly, the assessing model of highway slope stability based on PCA-GA-SVM is established. The maximal absolute error of the model is 0.0921 and the maximal relative error is 9.21% by comparing the assessment value and the practical value of the test sample. The above studies are conducive to enrich the assessing model of highway slope stability and provide some reference for highway slope engineering treatment.     

山东及邻区地震重新定位及其与活动构造的关系

Earthquakes taking place from 1975 to 2010 in and around Shandong Province are relocated using double-difference （HypoDD） and Hypoinvers 2000 （Hypo2000） methods, after correction of the onset times of seismic phases. The results show that the relocated seismicity is clearly associated with regional tectonics in space, and is also in agreement with the existence of deep faults imaged by wide-angle and deep seismic reflection profiling ; most of the focal depths are in the range of 5 - 25km, and there are clearly two predominant depths： 10km and 16km, which are inferred to be on the bottom of the upper crust and in the middle crust, respectively. The pattern of seismic activity indicates that moderate and strong earthquakes are likely to occur in the brittle-ductile transition zone between the upper and the lower crust, as the outcome of the deep tectonic dynamic process and the movement and deformation of faults in the upper and shallow crust under the regional stress field.     

沿107.6°E南北向剖面鄂尔多斯地块及周缘地区地壳结构

鄂尔多斯地块紧邻青藏高原东北缘,位于华北克拉通的西部,在我国中生代、新生代以来东部地区的构造活动中起到了重要作用.对鄂尔多斯及其周缘地区的研究可以提供有关华北克拉通的形成、演化和破坏过程的重要信息.本文选取了纵贯鄂尔多斯的107.6°E附近南北剖面上的44个流动地震台站进行分析,采用接收函数方法,进行Kirchhoff偏移成像,并且结合在该区域内前人的地震面波频散进行联合反演,获得剖面下方的地壳内部精细结构.研究结果显示:(1)莫霍面在鄂尔多斯北部较平缓,约45km深;在鄂尔多斯南部有所加深,达到50km;其北边的河套盆地的地壳厚度约为50km;南边的渭河盆地到秦岭地区及四川盆地的地壳厚度从约为40km增厚到47～50km.(2)河套盆地下方存在大规模的低速异常,最深可达25km,反映了其显著的拉张构造和沉积历史.(3)秦岭造山带下方的低速异常对应于其主要为长英质的地壳组分,可能是由于中生代的拆沉作用导致的地壳下部基性岩石层的缺失.(4)以38°N为界的鄂尔多斯地块,南北部地壳速度结构存在差异,可能表明了这两部分经历的构造历史不同.     

空-地-井重力异常正则化协同密度反演方法

围绕综合利用空-地-井多源重力异常联合反演提升精度的目标,提出正则化协同密度反演方法,从而有效利用多维数据的横向和纵向变化特征提高反演精度,且无需先验信息的约束.此外,还提出利用奇异值谱和深度分辨率工具评价航空数据观测高度和层数对反演分辨率的影响.通过理论模型试验本文方法在空-地、地-井等不同情况下的应用效果,结果表明空-地-井重力异常正则化协同密度反演方法能获得更高分辨率、高精度的反演结果,且证明不同的钻孔位置对反演结果有较大的影响,从而可指导实际的空-地-井联合勘探.最后,利用文顿盐丘实际勘探数据进行空地和空地井重力数据协同反演,反演结果垂向分辨率明显优于地面观测数据反演结果,盐盖的顶面及中心埋深与前人研究和地质资料解释相吻合,验证了方法正确性及实用性,为推进我国空-地-井立体重磁勘探提供了重要的技术手段.     

Incorporating known petrophysical model in the seismic full-waveform inversion using the Gramian constraint

In this paper, we develop a general approach to integrating petrophysical models in three-dimensional seismic full-waveform inversion based on the Gramian constraints. In the framework of this approach, we present an example of the frequency-domain P-wave velocity inversion guided by an electrical conductivity model. In order to introduce a coupling between the two models, we minimize the corresponding Gramian functional, which is included in the Tikhonov parametric functional. We demonstrate that in the case of a single-physics inversion guided by a model of different physical type, the general expressions of the Gramian functional and its gradients become simple and easy to program. We also prove that the Gramian functional has a non-negative quadratic form, so it can be easily incorporated in a standard gradient-based minimization scheme. The developed new approach of seismic inversion guided by the known petrophysical model has been validated by three-dimensional inversion of synthetic seismic data generated for a realistic three-dimensional model of the subsurface.     

Multiscale phase inversion for 3D ocean-bottom cable data

We invert three-dimensional seismic data by a multiscale phase inversion scheme, a modified version of full waveform inversion, which applies higher order integrations to the input signal to produce low-boost signals. These low-boost signals are used as the input data for the early iterations, and lower order integrations are computed at the later iterations. The advantages of multiscale phase inversion are that it (1) is less dependent on the initial model compared to full waveform inversion, (2) is less sensitive to incorrectly modelled magnitudes and (3) employs a simple and natural frequency shaping filtering. For a layered model with a three-dimensional velocity anomaly, results with synthetic data show that multiscale phase inversion can sometimes provide a noticeably more accurate velocity profile than full waveform inversion. Results with the Society of Exploration Geophysicists/European Association of Geoscientists and Engineers overthrust model shows that multiscale phase inversion more clearly resolves meandering channels in the depth slices. However, the data and model misfit functions achieve about the same values after 50 iterations. The results with three-dimensional ocean-bottom cable data show that, compared to the full waveform inversion tomogram, the three-dimensional multiscale phase inversion tomogram provides a better match to the well log, and better flattens angle-domain common image gathers. The problem is that the tomograms at the well log provide an incomplete low-wavenumber estimate of the log's velocity profile. Therefore, a good low-wavenumber estimate of the velocity model is still needed for an accurate multiscale phase inversion tomogram.     

The effectiveness of a pseudo-inverse extended Born operator to handle lateral heterogeneity for imaging and velocity analysis applications

Wave equation–based migration velocity analysis techniques aim to construct a kinematically accurate velocity model for imaging or as an initial model for full waveform inversion applications. The most popular wave equation–based migration velocity analysis method is differential semblance optimization, where the velocity model is iteratively updated by minimizing the unfocused energy in an extended image volume. However, differential semblance optimization suffers from artefacts, courtesy of the adjoint operator used in imaging, leading to poor convergence. Recent findings show that true amplitude imaging plays a significant role in enhancing the differential semblance optimization's gradient and reducing the artefacts. Here, we focus on a pseudo-inverse operator to the horizontally extended Born as a true amplitude imaging operator. For laterally inhomogeneous models, the operator required a derivative with respect to a vertical shift. Extending the image vertically to evaluate such a derivative is costly and impractical. The inverse operator can be simplified in laterally homogeneous models. We derive an extension of the approach to apply the full inverse formula and evaluate the derivative efficiently. We simplified the implementation by applying the derivative to the imaging condition and utilize the relationship between the source and receiver wavefields and the vertical shift. Specifically, we verify the effectiveness of the approach using the Marmousi model and show that the term required for the lateral inhomogeneity treatment has a relatively small impact on the results for many cases. We then apply the operator in differential semblance optimization and invert for an accurate macro-velocity model, which can serve as an initial velocity model for full waveform inversion.     

Estimation of the 3D correlation structure of an alluvial aquifer from surface-based multi-frequency ground-penetrating radar reflection data

Knowledge about the stochastic nature of heterogeneity in subsurface hydraulic properties is critical for aquifer characterization and the corresponding prediction of groundwater flow and contaminant transport. Whereas the vertical correlation structure of the heterogeneity is often well constrained by borehole information, the lateral correlation structure is generally unknown because the spacing between boreholes is too large to allow for its meaningful inference. There is, however, evidence to suggest that information on the lateral correlation structure may be extracted from the correlation statistics of the subsurface reflectivity structure imaged by surface-based ground-penetrating radar measurements. To date, case studies involving this approach have been limited to 2D profiles acquired at a single antenna centre frequency in areas with limited complementary information. As a result, the practical reliability of this methodology has been difficult to assess. Here, we extend previous work to 3D and consider reflection ground-penetrating radar data acquired using two antenna centre frequencies at the extensively explored and well-constrained Boise Hydrogeophysical Research Site. We find that the results obtained using the two ground-penetrating radar frequencies are consistent with each other, as well as with information from a number of other studies at the Boise Hydrogeophysical Research Site. In addition, contrary to previous 2D work, our results indicate that the surface-based reflection ground-penetrating radar data are not only sensitive to the aspect ratio of the underlying heterogeneity, but also, albeit to a lesser extent, to the so-called Hurst number, which is a key parameter characterizing the local variability of the fine-scale structure.     

Approximate traveltime inversion in downhole microseismic monitoring

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.     

联合GNSS和InSAR观测位移反演2008年汶川大地震断层位错模型参数

利用现代空间大地测量技术,尤其是卫星合成孔径雷达干涉测量,能够获取高精度、高空间分辨率的同震和孕震形变,为地震断层形变和破裂机制研究提供了前所未有的机遇。本文介绍了利用大地测量观测数据反演地震断层位错模型参数的贝叶斯反演方法。联合运用2008汶川大地震前后GNSS和InSAR技术观测获得的同震位移,反演了地震断层的几何参数和滑动位错分布。研究结果表明,汶川地震的断层滑动主要集中在倾角较陡的浅部,同时包含逆冲和右旋走滑,其中最大逆冲6.1m,最大右旋6.5m。根据断层滑动分布正演计算得到的上盘同震位移明显小于下盘,预示该断层两侧孕震形变可能存在较大的不对称性。