Seismic channel edge detection using 3D shearlets—a study on synthetic and real channelised 3D seismic data

Automatic feature detection from seismic data is a demanding task in today's interpretation workstations. Channels are among important stratigraphic features in seismic data both due to their reservoir capability or drilling hazard potential. Shearlet transform as a multi‐scale and multi‐directional transformation is capable of detecting anisotropic singularities in two and higher dimensional data. Channels occur as edges in seismic data, which can be detected based on maximizing the shearlet coefficients through all sub‐volumes at the finest scale of decomposition. The detected edges may require further refinement through the application of a thinning methodology. In this study, a three‐dimensional, pyramid‐adapted, compactly supported shearlet transform was applied to synthetic and real channelised, three‐dimensional post‐stack seismic data in order to decompose the data into different scales and directions for the purpose of channel boundary detection. In order to be able to compare the edge detection results based on three‐dimensional shearlet transform with some famous gradient‐based edge detectors, such as Sobel and Canny, a thresholding scheme is necessary. In both synthetic and real data examples, the three‐dimensional shearlet edge detection algorithm outperformed Sobel and Canny operators even in the presence of Gaussian random noise.     

基于无人机倾斜影像的三维建筑物震害精细信息提取

无人机倾斜摄影技术建模生成的三维影像较好地展现了建筑物侧面和顶面的震害细节信息,然而影像的高维度特性难以直接基于三维影像提取震害信息,经过降低维度转换的二维纹理影像往往会导致建筑物震害信息的不完整性和破碎性。针对这些问题,本文以2017年九寨沟M_S7.0地震为例,提出了一种直接从九寨沟震后三维影像获取侧面纹理信息的方法,即将三维模型打散,实现纹理与不规则三角网分离,从而获取完整的纹理影像,然后利用金字塔模型的瓦片坐标范围、瓦片命名规则和建筑物单体的空间位置选取最优纹理影像,再使用加权均值方差法确定纹理影像中建筑物的外墙最佳分割尺度后,采用面向对象方法提取建筑物外墙和墙皮脱落信息,最后通过对这些建筑物震害特征的分析,判定单体建筑物的破坏等级。结果显示,该方法成功获取了建筑物完整的侧面震害纹理影像,并基于纹理影像提取了外墙、裂缝和墙皮脱落区域信息判定建筑物单体为中等、严重两个破坏等级。      

Effective solutions to a global 3D visual system in networking environments

Today, with the proposal of digital earth, digital city and digital community, and rapid advances in informa-tion technologies and collection of spatial data, to a great degree, the traditional 2D GIS is promoted to Web and 3D GIS, which greatly challenges the capaci-ties of the system in transmitting and handling massive datasets. A great number of 3D GIS prototype systems have been developed and applied in such areas as mine resource management, digital cities. However, appli-cation sys…     

Validation of 3D synthetic seismograms based on the ray-Born approximation

The first-order Born approximation is a weak scattering perturbation method which is a powerful tool. The combination of the Born approximation and the ray theory enables to extend the applicability of the ray theory in terms of the required smoothness of the model and ensures faster computations than with, e.g., the finite difference method. We are motivated to describe and explain the effects of the numerical discretization of the Born integral on the resulting seismograms. We focus on forward modelling and study the cases in which perturbation from the background model contains the interface. We restrict ourselves to isotropic models that contain two homogeneous layers. We compare the 2D and 3D ray-based Bornapproximation seismograms with the ray-theory seismograms. The Born seismograms are computed using a grid of finite extent. We anticipate that the computational grid should contain an appropriate number of gridpoints, otherwise the seismogram would be inaccurate. We also anticipate that the limited size of the computational grid can cause problems. We demonstrate numerically that an incorrect grid can produce significant errors in the amplitude of the wave, or it can shift the seismogram in time. Moreover, the grid boundaries work as interfaces, where spurious waves can be generated. We also attempt to explain these phenomena theoretically. We give and test the options of removing the spurious waves. We show that it is possible to compute the Born approximation in a sparser grid, if we use elastic parameters averaged from some dense grid.     

地震信息二、三维一体化系统设计

运用Echarts可视化插件和全景技术,将GIS技术与地震信息服务进行有效结合,构建基于WebGIS的地震信息可视化地图显示、标注、查询等功能,实现历史大规模地震目录数据和台站信息在可视化系统中的二、三维一体化展示效果,通过发布震源深度参数的分布图、三维动态折纸模型和全景地震台站,加强公众结合区域特征观看和理解地震活动的直观性和感官效果。     

Lake ice simulation using a 3D unstructured grid model

Ocean Dynamics - We develop a single-class ice and snow model embedded inside a 3D hydrodynamic model on unstructured grids and apply it to lake studies using highly variable mesh resolution. The...     

A probabilistic framework to make a decision on the post-earthquake functionality of bridges considering the damage,residual displacement,and aftershock

Bulletin of Earthquake Engineering - Making a decision on the post-earthquake functionality of structures has always been one of the most challenging issues in earthquake engineering. There is an...     

Analytical study of performance evaluation for seismic retrofitting of reinforced concrete building using 3D dynamic nonlinear finite element analysis

This paper presents three-dimensional fi nite element (FE) analyses of an all-frame model of a three-story reinforced concrete (RC) building damaged in the 1999 Taiwan Chi-Chi Earthquake. Non-structural brick walls of the building acted as a seismic resistant element although their contributions were neglected in the design. Hence, the entire structure of a typical frame was modeled and static and dynamic nonlinear analyses were conducted to evaluate the contributions of the brick walls. However, the results of the analyses were considerably overestimated due to coarse mesh discretizations, which were unavoidable due to limited computer resources. This study corrects the overestimations by modifying (1) the tensile strengths and (2) shear stiffness reduction factors of concrete and brick. The results indicate that brick walls improve frame strength although shear failures are caused in columns shortened by spandrel walls. Then, the effectiveness of three types of seismic retrofi ts is evaluated. The maximum drift of the firstoor is reduced by 89.3%, 94.8%, and 27.5% by Steel-confi ned, Full-RC, and Full-brick models, respectively. Finally, feasibility analyses of models with soils were conducted. The analyses indicated that the soils elongate the natural period of building models although no signifi cant differences were observed.     

3D resistivity inversion using 2D measurements of the electric field

Field and 'noisy' synthetic measurements of electric-field components have been inverted into 3D resistivities by smoothness-constrained inversion. Values of electrical field can incorporate changes in polarity of the measured potential differences seen when 2D electrode arrays are used with heterogeneous 'geology', without utilizing negative apparent resistivities or singular geometrical factors. Using both the X - and Y -components of the electric field as measurements resulted in faster convergence of the smoothness-constrained inversion compared with using one component alone. Geological structure and resistivity were reconstructed as well as, or better than, comparable published examples based on traditional measurement types. A 2D electrode grid (20 × 10), incorporating 12 current-source electrodes, was used for both the practical and numerical experiments; this resulted in 366 measurements being made for each current-electrode configuration. Consequently, when using this array for practical field surveys, 366 measurements could be acquired simultaneously, making the upper limit on the speed of acquisition an order of magnitude faster than a comparable conventional pole–dipole survey. Other practical advantages accrue from the closely spaced potential dipoles being insensitive to common-mode noise (e.g. telluric) and only 7% of the electrodes (i.e. those used as current sources) being susceptible to recently reported electrode charge-up effects.     

利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性

大地电磁野外实测数据目前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解释二维剖面数据.对大地电磁二维剖面的张量数据进行三维反演,不仅可以沿剖面获得较好的二维断面结果,还能够得到二维反演所不能获得的剖面附近的三维电阻率结构信息.合成数据的反演算例表明:对二维剖面数据进行三维反演时,对角元素对于圈定剖面附近三维异常体的空间分布具有独特作用,应尽量反演所有的张量元素.     

Calculations of hydrodynamic time parameters in a semi-opened coastal zone using a 3D hydrodynamic model

Hydrodynamic time parameters (HTs) in a semi-opened aquatic ecosystem are synthetic indicators offering the opportunity to bring out the links between its physical functioning and its biology. The generic term “residence time” is frequently used through literature to mention HTs resulting in various calculation methods. This article presents different computing methods relying on the use of a 3D numerical hydrodynamic model and the HTs to which they give access. Several large-scale (water exchange time, average water export time, e-flushing time) and local time parameters (export time, flushing lag, local e-flushing time) are defined. The applications presented are carried out within the south-west lagoon of New Caledonia (SLNC), on three embedded control volumes. The definition of the control volume is more important for the values of local HTs than for their comparative distribution. The comparison of the global hydrodynamic time scales applied to a control volume provides information on the mixing processes inside the control volume.     

3D multi-fractal analysis of porous media using 3D digital images: considerations for heterogeneity evaluation

Pore structure heterogeneity is a critical parameter controlling mechanical, electrical and flow transport behaviour of rock. Multi-fractal analysis method was used for a heterogeneity comparison of three-dimensional rock samples with different lithology. Six real digital samples, containing three sandstones and three carbonates, were used. Based on the mercury injection capillary pressure test on these samples, we found that the carbonate samples are more heterogeneous than sandstones, but primary results of multi-fractal behaviours for all samples were similar. We show that if multi-fractal is used to evaluate and compare heterogeneity of different samples, one needs to follow some considerations such as (1) all samples must have the same size in pixel, (2) samples volume must be bigger than representative volume element, (3) multi-fractal dimensions should be firstly normalized to a determined porosity value and (4) multi-fractal results should be interpreted based on resolution of the imaging tool (effects of fine scale sub-resolution pores are missed). Results revealed that using normalized fractal dimensions, the real samples were divided to less and high heterogeneous groups. Moreover, the study of scale effect also showed that porous structures of these samples are scale invariant in a wide range of scales (from one to eight times bigger).     

基于局部一维模型与走时标定的区域 三维速度模型构建技术研究

构建区域介质三维速度模型并以之获得准确的区域震相走时, 是提高区域地震定位精度的重要手段之一. 为充分利用已有的一维模型、 GT事件、 地质资料等实现三维模型构建, 尝试基于目标区域内已有的部分局部一维模型, 通过克里金空间插值建立初始三维模型, 然后利用GT事件走时数据并参考其它地震地质资料对其不断进行修正, 使得其走时偏差图与GT事件走时偏差图一致, 进而获得能够提高区域地震定位精度的三维模型. 使用不同模型进行的地震定位实验表明, 以此方法建立的三维模型的定位偏差较初始模型减少约20%, 较好地起到了减小区域震相走时残差, 提高区域地震定位精度的作用.     

基于局部一维模型与走时标定的区域三维速度模型构建技术研究

构建区域介质三维速度模型并以之获得准确的区域震相走时,是提高区域地震定位精度的重要手段之一.为充分利用已有的一维模型、GT事件、地质资料等实现三维模型构建,尝试基于目标区域内已有的部分局部一维模型,通过克里金空间插值建立初始三维模型,然后利用GT事件走时数据并参考其它地震地质资料对其不断进行修正,使得其走时偏差图与GT事件走时偏差图一致,进而获得能够提高区域地震定位精度的三维模型.使用不同模型进行的地震定位实验表明,以此方法建立的三维模型的定位偏差较初始模型减少约20％,较好地起到了减小区域震相走时残差,提高区域地震定位精度的作用.     

迭代法在三维网络模型最短路径射线追踪中的应用

目前有许多基于二维网格模型的射线追踪算法,但应用与三维网格模型的射线追踪算法较少,主要是因为三维网格射线追踪算法在算法实现、编程以及图形显示方面都有一定难度.本文提出了一种适用于三维网格射线追踪的迭代算法,并应用于几种理论模型,同时实现了三维绘图显示,通过研究证明迭代算法能在一定程度上提高三维网格射线追踪的精度和速度.     

3D inversion of DC data using artificial neural networks

In this paper, we investigate the applicability of artificial neural networks in inverting three-dimensional DC resistivity imaging data. The model used to produce synthetic data for training the artificial neural network (ANN) system was a homogeneous medium of resistivity 100 Ωm with an embedded anomalous body of resistivity 1000 Ωm. The different sizes for anomalous body were selected and their location was changed to different positions within the homogeneous model mesh elements. The 3D data set was generated using a finite element forward modeling code through standard 3D modeling software. We investigated different learning paradigms in the training process of the neural network. Resilient propagation was more efficient than any other paradigm. We studied the effect of the data type used on neural network inversion and found that the use of location and the apparent resistivity of data points as the input and corresponding true resistivity as the output of networks produces satisfactory results. We also investigated the effect of the training data pool volume on the inversion properties. We created several synthetic data sets to study the interpolation and extrapolation properties of the ANN. The range of 100–1000 Ωm was divided into six resistivity values as the background resistivity and different resistivity values were also used for the anomalous body. Results from numerous neural network tests indicate that the neural network possesses sufficient interpolation and extrapolation abilities with the selected volume of training data. The trained network was also applied on a real field dataset, collected by a pole-pole array using a square grid (8 ×8) with a 2-m electrode spacing. The inversion results demonstrate that the trained network was able to invert three-dimensional electrical resistivity imaging data. The interpreted results of neural network also agree with the known information about the investigation area.     

3D seismic response of a limited valley via BEM using 2.5D analytical Green's functions for an infinite free-rigid layer

This paper presents analytical solutions for computing the 3D displacements in a flat solid elastic stratum bounded by a rigid base, when it is subjected to spatially sinusoidal harmonic line loads. These functions are also used as Greens functions in a boundary element method code that simulates the seismic wave propagation in a confined or semi-confined 2D valley, avoiding the discretization of the free and rigid horizontal boundaries.The models developed are then used to simulate wave propagation within a rigid stratum and valleys with different dimensions and geometries, when struck by a spatially sinusoidal harmonic vertical line load. Simulations are performed in the frequency domain, for varying spatial wave numbers in the axial direction of the valley. Time results are obtained by means of inverse Fourier transforms, to help understand how the geometry of the valley may affect the variation of the displacement field.     

Numerical simulation of organic chemicals in a marine environment using a coupled 3D hydrodynamic and ecotoxicological model

For ecotoxicological risk assessment in a marine ecosystem, we constructed a coupled three-dimensional hydrodynamic and ecotoxicological model (EMT-3D), and applied it to Tokyo Bay. The model was calibrated with field data obtained in 2002. The results of sensitivity analysis for dissolved Bisphenol A showed that biodegradation rate was the most important factor for concentration change. Bioconcentration coefficient was the most important factor for Bisphenol A in phytoplankton. Therefore, the parameters must be carefully considered in the modeling. The mass balance results showed that standing stocks of Bisphenol A in water, in particulate organic carbon and in phytoplankton are 7.85 x 10(4), 1.78 x 10(2) and 3.44 x 10(-1) g, respectively. With respect to flux, biodegradation in the water column had the highest value of 1.06 x 10(3) g/day, and next were effluent to the open sea, partition to particulate organic carbon, and bioconcentration in phytoplankton.