Seismoelectric numerical simulation in 2D vertical transverse isotropic poroelastic medium

Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ.     

曲流河沉积亚相的地震识别方法

河道边缘的识别是河流相储层精细描述的重点。受河道叠置、交叉,砂体厚度薄,地震信噪比、分辨率低等因素影响,传统的切片解释、相干技术难以满足精细勘探的需求,新发展的基于算子处理的边缘检测仍存在应用误区。本文以分析河道边缘的几何特征为切入点,明确一阶导数、模值及二阶导数的物理含义;建立具有不同速度特征的三维河道模型,并通过模型与实际资料相干属性的提取,指出该项技术在河道精细描述中存在的问题。针对该问题以Sobel算子为例,图解说明使用此项技术处理后河道边缘的标志性特征,提出采用直方图均衡化与模糊集理论的河道边缘相干增强技术,并取得较好的应用效果。该方法技术对深化河道边缘特征的认识,提高河流相储层识别能力,有一定的借鉴作用。     

Three-dimensional numerical simulation of glacial trough forming process

The glacial trough is a common glacier erosion landscape, which plays an important role in the study of glacier erosion processes. In a sharp contrast with the developing river, which is generally meandering, the developing glacial trough is usually wide and straight. Is the straightness of the glacial trough just the special phenomenon of some areas or a universal feature? What controls the straightness of the glacial trough? Until now, these issues have not been studied yet. In this paper, we conduct systematic numerical models of the glacier erosion and simulate the erosion evolution process of the glacial trough. Numerical simulations show that:(1) while the meandering glacier is eroding deeper to form the U-shaped cross section, the glacier is eroding laterally. The erosion rate of the ice-facing slope is bigger than that of the back-slope.(2) The smaller(bigger) the slope is, the smaller(bigger) the glacier erosion intensity is.(3) The smaller(bigger) the ice discharge is, the smaller(bigger) the glacier erosion intensity is. In the glacier erosion process, the erosion rate of the ice-facing slope is always greater than that of the back-slope. Therefore, the glacial trough always develops into more straight form. This paper comes to the conclusion that the shape evolution of the glacial trough is controlled mainly by the erosion mechanism of the glacier. Thereby, the glacial trough prefers straight geometry.     

Formation of a cohesive floodplain in a dynamic experimental meandering river

Field studies suggest that a cohesive floodplain is a necessary condition for meandering in contrast to braided rivers. However, it is only partly understood how the balance between floodplain construction by overbank deposition and removal by bank erosion and chutes leads to meandering. This is needed because only then does a dynamic equilibrium exist and channels maintain meandering with low width–depth ratios. Our objective is to understand how different styles of floodplain formation such as overbank deposition and lateral accretion cause narrower channels and prevent chute cutoffs that lead to meandering. In this study we present two experiments with a self‐forming channel in identical conditions, but to one we added cohesive silt at the upstream boundary. The effect of cohesive silt on bank stability was tested in auxiliary bank erosion experiments and showed that an increase in silt reduced erosion rates by a factor of 2. The experiment without silt developed to a braided river by continuous and extensive shifting of multiple channels. In contrast, in the meandering river silt deposits increased bank stability of the cohesive floodplain and resulted in a reduction of chute cutoffs and increased sinuosity by continuous lateral migration of a single channel. Overbank flow led to deposition of the silt and two styles of cohesive floodplain were observed: first, overbank vertical‐accretion of silt, e.g. levee, overbank sedimentation or splays; and second, lateral point bar accretion with silt on the scrolls and in the swales. The first style led to a reduction in bank erosion, while the second style reduced excavation of chutes. We conclude that sedimentation of fine cohesive material on the floodplain by discharge exceeding bankfull is a necessary condition for meandering. Copyright © 2013 John Wiley & Sons, Ltd.     

起伏地形下各向异性的2D大地电磁无网格法数值模拟

作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明：无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法.     

轨道不平顺随机过程的数值模拟

轨道不平顺是列车-轨道系统的动力响应的主要影响因素,不平顺样本是列车-轨道系统动力分析模型中不可缺少的参数.根据随机理论原理推导了轨道不平顺谱的空间、时间域的相互转化关系,得到了时域功率谱的表达式.分析了目前广泛应用几种轨道不平顺随机过程的模拟方法的原理及作法,通过对功率谱快速数值算法的分析给出了一种新三角级数模拟方法,该方法与逆傅氏变换法具有等效性,同时也论证了逆傅氏变换方法样本是高斯过程.最后通过实际算例的分析论证了新三角级数法与傅氏逆变换法在简便、合理等方面的优点.     

Modulation of the flow structure by progressive bedforms in the Kinoshita meandering channel

An in‐house fully three‐dimensional general‐purpose finite element model is applied to solve the hydrodynamic structure in a periodic Kinoshita‐generated meandering channel. The numerical model solves the incompressible Reynolds‐averaged Navier–Stokes equations for mass and momentum, while solving the k ? ε equations for turbulence. The free surface is described by the rigid‐lid approximation (using measured water surface data) for flat (smooth‐bed) and self‐formed (rough‐bed) conditions. The model results are compared against experimental measurements in the ‘Kinoshita channel’, where three‐dimensional flow velocities and turbulence parameters were measured. This validation was carried out for the upstream‐valley meander bend orientation under smooth (flat bed) conditions. After validation, several simulations were carried out to predict the hydrodynamics in conditions where either it was not possible to perform measurements (e.g. applicability of the laboratory acoustic instruments) and to extrapolate the model to other planform configurations. For the flat smooth‐bed case, a symmetric (no skewness) planform configuration was modeled and compared to the upstream‐skewed case. For the self‐formed rough‐bed case, prediction of the hydrodynamics during the progression of bedforms was performed. It appears that the presence of bedforms on a bend has the following effects: (i) the natural secondary flow of the bend is disrupted by the presence of the bedforms, thus depending on the location of the dune, secondary flows might differ completely from the traditional orientation; (ii) an increment on both the bed and bank shear stresses is induced, having as much as 50% more fluvial erosion, and thus a potential increment on the migration rate of the bend. Implications on sediment transport and bend morphodynamics are also discussed in the paper. Copyright © 2013 John Wiley & Sons, Ltd.     

3D SPH numerical simulation of the wave generated by the Vajont rockslide

A 3D numerical modeling of the wave generated by the Vajont slide, one of the most destructive ever occurred, is presented in this paper. A meshless Lagrangian Smoothed Particle Hydrodynamics (SPH) technique was adopted to simulate the highly fragmented violent flow generated by the falling slide in the artificial reservoir. The speed-up achievable via General Purpose Graphic Processing Units (GP-GPU) allowed to adopt the adequate resolution to describe the phenomenon. The comparison with the data available in literature showed that the results of the numerical simulation reproduce satisfactorily the maximum run-up, also the water surface elevation in the residual lake after the event.Moreover, the 3D velocity field of the flow during the event and the discharge hydrograph which overtopped the dam, were obtained.     

Spatial random field simulation by a numerical series representation

An approach to the simulation of spatial random fields is proposed. The target random field is specified by its covariance function which need not be homogeneous or Gaussian. The technique provided is based on an approximate Karhunen–Loève expansion of spatial random fields which can be readily realized. Such an approximate representation is obtained from a correction to the Rayleigh–Ritz method based on the dual Riesz basis theory. The resulting numerical projection procedure improves Rayleigh–Ritz algorithm in the approximation of second-order random fields. Simulations are developed to illustrate the convergence and accuracy of the method presented.

Calculation of primary and secondary flow and boundary shear stresses in a meandering channel

Turbulent flow in a meandering channel is computed with two Computational Fluid Dynamics (CFD) codes solving the Navier–Stokes equations by employing different turbulence closure approaches. The first CFD code solves the steady Reynolds-Averaged Navier–Stokes equations (RANS) using an isotropic turbulence closure. The second code is based on the concept of Large Eddy Simulation (LES). LES resolves the large-scale turbulence structures in the flow and is known to outperform RANS models in flows in which large-scale structures dominate the statistics. The results obtained from the two codes are compared with experimental data from a physical model study. Both, LES and RANS simulation, predict the primary helical flow pattern in the meander as well as the occurrence of an outer-bank secondary cell. Computed primary as well as secondary flow velocities are in reasonably good agreement with experimental data. Evidence is given that the outer-bank secondary cell in a meander bend is the residual of the main secondary cell of the previous bend. However, the RANS code, regardless of the turbulence model employed, overpredicts the size and strength of the outer-bank secondary cell. Furthermore, only LES is able to uphold the outer-bank second secondary cell beyond the bend apex until the exit of the bend as turbulence anisotropy contributes to its persistence. The presence of multiple secondary cells has important consequences for the distribution of shear stresses along the wetted perimeter of the channel, and thereby the sediment transport in meandering channels. Consequently, even though LES is expected to compute the bed-shear stresses along the wetted perimeter of the channel with a higher degree of accuracy than the RANS model, comparisons between LES and RANS computed wall shear stresses agree well. These findings are useful for practitioners who need to rely on RANS model predictions of the flow in meandering channels at field scale.     

3D numerical simulation of flow field around twin piles

In this study to identify the flow pattern and local scour mechanism around pile groups, the flow field was simulated using FLOW-3D software. A pair of pile on a flat-bed channel with side by side and tandem arrangements was investigated. To establish Navier–Stokes equations, the RNGk-ε turbulence model was used and the results were verified using experimental data. In case of FLOW-3D capability, it was found that the software was able to properly simulate the expected interaction between the pile groups. The results of flow field simulation showed that Reynolds number and the pile spacing are the most influential variables in forming vortices. The flow around tandem pile and the downward flow around wake vortices were more intense and complicate in comparison with side by side arrangements and single pile.     

隐伏正断层错动引发上覆土体破裂过程的三维数值模拟

通过建立三维数值模型,对隐伏正断层在均匀错动和倾斜错动方式下土体的破裂过程进行研究。利用应力罗德参数和等效塑性应变分别对断层错动过程中上覆土体的应力状态和破坏形式进行分析,并提出土体破裂的判别方法。通过对数值模拟结果的分析得到以下结论:① 在断层错动过程中,下盘一侧受断层错动影响的上覆土体的应力状态经压剪→纯剪→拉剪逐渐变化,而上盘一侧上覆土体的应力状态变化较为复杂,经压剪→纯剪→拉剪→纯剪→压剪重复变化;② 在断层均匀错动过程中,断层下盘一侧土体的破裂率先出现在地表拉剪区内,随错动量的增大,破裂带向两侧、向深部扩展;同时,下盘一侧土体的底部产生破坏,并斜向上扩展,逐渐与顶部破裂相连;③ 在断层倾斜错动过程中,地表破裂出现的位置和上覆土体的厚度有关。对于厚度较大的土体,正断层倾斜错动能够在地表形成与断层走向有一定夹角、且与断层长度相比长度很短的地表破裂或地裂缝,而数值模拟可对正断层错动导致的地表破裂的模式加以补充,为研究地裂缝的形成机理和分布形式提供依据。      

震源动力学破裂过程数值模拟研究

首先阐述了震源动力学过程研究的重要意义, 在此基础上, 研究了复杂的断层几何形态及介质模型对动力学破裂过程的影响, 并对常用的有限元方法、 离散元方法、 有限差分方法和边界积分方程方法等进行了相应介绍. 讨论了这些数值模拟方法各自的优缺点, 建议在方法的选择上应视具体问题及计算的精度而定. 最后对动力学数值模拟的关键部分, 滑动摩擦准则进行了论述. 常用的滑动摩擦准则有滑动弱化准则、 速率弱化准则和速率-状态依赖摩擦准则. 在单纯考虑某个地震的动力学破裂传播过程时, 滑动弱化准则较为常用, 其中滑动弱化距离的选取至关重要. 但若考虑整个地震循环, 速率-状态依赖摩擦准则更为合适.      

基于BISQ模型的三维双相各向异性介质数值模拟

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.     

东北冷涡引起的对流输送过程个例分析和数值模拟

平流层-对流层物质交换是影响全球大气成分收支的重要过程.过去的研究认为大尺度的交换过程在平流层-对流层物质交换中最为重要,但是近些年的研究表明,中小尺度过程对平流层-对流层物质交换也有重要贡献.本文利用OMI和MLS数据、ERA-Interim再分析资料,结合中尺度WRF模式综合分析了东北地区发生在冷涡前部和冷涡后部的两次强对流天气过程.结果表明：发生在冷涡前部暖锋云系中的强对流持续时间长,对流垂直尺度小,下平流层静力稳定度高;发生在冷涡后部的孤立强对流持续时间短,水平尺度较小,且在对流层顶附近,静力稳定度小,对流可穿出热力学对流层顶.从示踪物分布情况来看,两次强对流都可将示踪物输送到对流层顶附近,但是冷涡前部对流可将示踪物从边界层输送到整个对流层,而孤立对流是把示踪物输送到对流层顶,而不与自由对流层空气发生混合.

     

Resolution of 2D Wenner resistivity imaging as assessed by numerical modelling

Modelling of 2D resistivity imaging was done in order to understand the principle resolution of the technique in different geological situations, and for assessing the behaviour of the interpretation methods under controlled circumstances. The Wenner array was used throughout. The results show that the 1D approximation only provides reasonable results in environments with very gradual lateral resistivity changes, otherwise the result may be strongly misleading. Inversion using the 2D quasi-Newton technique results in adequate resolution of the structures in moderately complex environments, but the Gauss–Newton method holds a significant advantage in some complicated cases. The data density can also be of crucial importance for the resolution capability, notably of narrow structures.     

基于非结构化网格的边界积分方程方法的断层自发破裂模拟

地震自发破裂模拟是震源动力学研究的重要内容,了解复杂的断层动力学破裂过程对深入认识震源特征和解释运动学反演结果具有重要意义.基于边界积分方程方法的破裂模拟已经被广泛使用,大多采用的是平面断层模型的结构化网格划分.由于实际的断层往往具有较为复杂的几何特征,为了更为灵活地刻画断层几何复杂性,我们建立断层模型的三角形网格离散方案,通过精确的解析解形式来计算断层各个单元之间的应力格林函数,联立滑动弱化摩擦准则和非奇异边界积分方程,对断层的自发破裂过程进行了模拟.在简单的平面断层模型下,将计算结果与前人的结果进行了对比,验证了方法的正确性与有效性.对于几种常见的复杂断层模型,例如弯折、阶跃、含障碍体断层等,我们模拟了其破裂过程并对计算结果进行了比较与分析.模拟结果表明,非结构化网格划分的边界积分方程方法能够很好地模拟平面矩形断层或由其组成的规则断层,同时也能成功地模拟具有复杂几何形状的不规则断层上的动力学破裂过程.本研究的结果显示了边界积分方程方法在模拟复杂断层系统的动力学破裂问题上具有较广阔的应用前景.

     

Adjustment of self-formed bankfull channel geometry of meandering rivers: modelling study

We consider the evolution of the hydraulic geometry of sand-bed meandering rivers. We study the difference between the timescale of longitudinal river profile adjustment and that of channel width and depth adjustment. We also study the effect of hydrological regime alteration on the evolution of bankfull channel geometry. To achieve this, a previously developed model for the spatiotemporal co-evolution of bankfull channel characteristics, including bankfull discharge, bankfull width, bankfull depth and down-channel bed slope, is used. In our modelling framework, flow variability is considered in terms of a specified flow duration curve. Taking advantage of this unique feature, we identify the flow range responsible for long-term bankfull channel change within the specified flow duration curve. That is, the relative importance of extremely high short-duration flows compared to moderately high longer duration flows is examined. The Minnesota River, MN, USA, an actively meandering sand-bed stream, is selected for a case study. The longitudinal profile of the study reach has been in adjustment toward equilibrium since the end of the last glaciation, while its bankfull cross-section is rapidly widening due to hydrological regime change in the last several decades. We use the model to demonstrate that the timescale for longitudinal channel profile adjustment is much greater than the timescale for cross-sectional profile adjustment due to a lateral channel shift. We also show that hydrological regime shift is responsible for the recent rapid widening of the Minnesota River. Our analysis suggests that increases in the 5–25% exceedance flows play a more significant role in recent bankfull channel enlargement of the Minnesota River than increase in either the 0.1% exceedance flow or the 90% exceedance flow. © 2020 John Wiley & Sons, Ltd.     

3D DEMAP 数值模拟与观测实验

为推动电磁剖面(EMAP)技术向找矿勘探领域发展, 本文针对密集阵列剖面(DEMAP)测量方式, 利用三维积分方程法对层状介质中赋存三维地质异常体的电性结构进行了数值模拟, 并在青海野马泉铁锌矿集区西部进行了大地电磁(MT)与DEMAP的野外对比观测实验.数值模拟结果显示, DEMAP观测方式获得的结果受偏移距影响, 造成的视电阻率和相位误差与偏移距大小和地下电性结构的复杂性密切相关, 但在整体上, 统计误差≤10%.通过对野马泉矿集区的实验结果分析可得, DEMAP与MT的视电阻率断面和相位断面具有较好的一致性.本文从理论和实验上都表明DEMAP观测技术是一种有潜力的、经济、快速、有效的找矿手段.     

矿井瞬变电磁法三维时域有限差分数值模拟

给出了矿井全空间瞬变电磁场的时域有限差分法(FDTD)算法,并推导了Mur吸收边界条件.利用Mur吸收边界条件,选用均匀全空间电偶极源作为初始激发源,模拟了均匀介质中巷道底板岩层内部和层状介质中三维低阻异常体的全空间响应特性,分析了瞬变电磁场在均匀介质中1.9μs和27.6μs两时刻的传播规律及在层状介质中1.2μs和0.023 ms两时刻的传播规律.结果表明:巷道对电场的影响在初期并不明显,在20μs后才表现出来;瞬变场能较好地分辨低阻层,对高阻层的穿透能力强,并且对低阻体反映灵敏;吸收边界条件只有在场域较大、网格节点较多时才能产生明显的效果.所研究成果为矿井瞬变电磁法资料的解释提供了理论依据.