Implementation of perfectly matched layers in an arbitrary geometrical boundary for elastic wave modelling   总被引：1，自引：0，他引：1  

Hongwei Gao  Jianfeng Zhang 《Geophysical Journal International》2008,174(3):1029-1036

The perfectly matched layer (PML) absorbing boundary condition is incorporated into an irregular-grid elastic-wave modelling scheme, thus resulting in an irregular-grid PML method. We develop the irregular-grid PML method using the local coordinate system based PML splitting equations and integral formulation of the PML equations. The irregular-grid PML method is implemented under a discretization of triangular grid cells, which has the ability to absorb incident waves in arbitrary directions. This allows the PML absorbing layer to be imposed along arbitrary geometrical boundaries. As a result, the computational domain can be constructed with smaller nodes, for instance, to represent the 2-D half-space by a semi-circle rather than a rectangle. By using a smooth artificial boundary, the irregular-grid PML method can also avoid the special treatments to the corners, which lead to complex computer implementations in the conventional PML method. We implement the irregular-grid PML method in both 2-D elastic isotropic and anisotropic media. The numerical simulations of a VTI lamb's problem, wave propagation in an isotropic elastic medium with curved surface and in a TTI medium demonstrate the good behaviour of the irregular-grid PML method.  相似文献   

Seismic radiation from dynamic coalescence, and the reconstruction of dynamic source parameters on a planar fault     

Nobuki Kame  Koji Uchida 《Geophysical Journal International》2008,174(2):696-706

The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω⁻²  -type displacement spectra.  相似文献   

Quantitative accuracy analysis of the discontinuous Galerkin method for seismic wave propagation   总被引：1，自引：0，他引：1  

Martin Käser  Verena Hermann  Josep de la Puente 《Geophysical Journal International》2008,173(3):990-999

  相似文献   

Development of a new semi-analytical model for cross-borehole flow experiments in fractured media     

《Advances in water resources》2015

Analysis of borehole flow logs is a valuable technique for identifying the presence of fractures in the subsurface and estimating properties such as fracture connectivity, transmissivity and storativity. However, such estimation requires the development of analytical and/or numerical modeling tools that are well adapted to the complexity of the problem. In this paper, we present a new semi-analytical formulation for cross-borehole flow in fractured media that links transient vertical-flow velocities measured in one or a series of observation wells during hydraulic forcing to the transmissivity and storativity of the fractures intersected by these wells. In comparison with existing models, our approach presents major improvements in terms of computational expense and potential adaptation to a variety of fracture and experimental configurations. After derivation of the formulation, we demonstrate its application in the context of sensitivity analysis for a relatively simple two-fracture synthetic problem, as well as for field-data analysis to investigate fracture connectivity and estimate fracture hydraulic properties. These applications provide important insights regarding (i) the strong sensitivity of fracture property estimates to the overall connectivity of the system; and (ii) the non-uniqueness of the corresponding inverse problem for realistic fracture configurations.  相似文献   

2013年洱源M_S5.5地震序列M_L≥3.4地震的震源机制解研究     

杨军  苏有锦  李孝宾  叶泵  陈佳 《地震研究》2015,38(2):196-202,333

利用云南省地震台网的宽频带数据,采用CAP方法计算了2013年3月3日洱源MS5.5地震前后共10次ML≥3.4地震的震源机制解,分析了洱源MS5.5地震前后震源机制解的一致性变化特征,结合精定位结果、震后现场灾害评估报告及震区地质构造情况,判定维西—巍山断裂以西的苍山西麓炼铁盆地东缘主边界断裂带为洱源MS5.5地震的发震构造,震源机制解节面Ⅱ(走向160°、倾角48°、滑动解-99°)为其发震断层面。  相似文献   

A new series solution method for two-dimensional elastic wave scattering along a canyon in half-space     

《Soil Dynamics and Earthquake Engineering》2016

This paper presents a semi-analytical method for studying the two-dimensional problem of elastic wave scattering by surface irregularities in a half-space. The new method makes use of the member of a c-completeness family of wave functions to construct the scattering fields, and then applies equal but opposite tractions to those of the foregoing constructed scattering fields on the horizontal surface of the half-space to produce additional scattering fields. These additional scattering fields are a series of Lamb's solutions. Thus the whole scattering field constructed in the series automatically satisfies the Navier equations, the condition of zero traction on the half-space surface, and the radiation boundary conditions at infinity. Using the traction-free conditions along the canyon surface, the coefficients of the series solutions are determined via a least-squares method. For incident P, SV, and Rayleigh waves, the numerical results are presented for the scattering displacements in the vicinity of a semi-circular canyon in the half-space.  相似文献   

Computer thermodynamic modeling of the transport and deposition of Sb and Au during the formation of Au-Sb deposits   总被引：1，自引：0，他引：1  

A.A. Obolensky  L.V. Gushchina  A.S. Borisenko  A.A. Borovikov  P.A. Nevol&#x;ko 《Russian Geology and Geophysics》2009,50(11):950-965

Using the Chiller computer program, we performed modeling of the mechanisms of the joint transport and deposition of Au and Sb from various ore-forming solutions during the formation of Au-Sb deposits. Three models are considered by the example of the Uderei Au-Sb deposit in the Yenisei Ridge: (1) simple cooling (cooling only), (2) iso-enthalpy boiling (P = f(T)), and (3) solution–rock interaction (rock titration model). The behavior of Sb(III) and Au(I) in the system Au–Sb–Fe–Cu–Pb–Zn–As–H₂O–Cl–H₂S–CO₂ under hydrothermal conditions was studied. It is shown that both weakly alkaline (near-neutral) and reduced acidic Fe_aq²⁺-enriched low-chloride high-CO₂ and high-chloride hydrothermal solutions play a crucial role in the formation of gold parageneses of Au-Sb ores.  相似文献   

On the frontal condition for finite amplitude α2Ω-dynamo wave trains in stellar shells     

Andrew P. Bassom  Andrew M. Soward 《地球物理与天体物理流体动力学》2013,107(3-4):285-328

Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θ_F. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θ_f; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α²Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θ_f—θ_F) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D _c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α²Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R _α, explicitly ε^½ ? R ² _α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θ_f—θ_F small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D _min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D _min| ≤ |D _c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D _min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R _{2 α} ? ε^½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.  相似文献   

混凝土膨胀剂使用中存在的问题及解决措施     

莫庆瑗 《中国勘察设计》2012,(7):80-82

混凝土膨胀剂具有补偿混凝土干燥收缩的效果,但在地下室防水混凝土工程中,使用效果不理想,存在掺加膨胀剂后地下室底板和外墙出现混凝土开裂渗漏的现象。本文从设计源头,力求全面地、准确有效地提出解决措施。  相似文献   

H020U井钻进所遇问题及处理方法          下载免费PDF全文

张新刚  隆东  向军文  林修阔 《探矿工程》2011,38(1):28-31

土耳其贝帕扎里天然碱矿H020U井在造斜钻进和水平钻进过程中遇到井漏、钻头钻遇硅化灰岩层磨损严重、未绕过第一靶点V020而打中套管、井涌、仪器掉入裸眼、钻遇未知溶腔导致顶角下降,方位偏离而导致未直接连通第二靶点等问题。针对这些问题进行了分析,并提出了相应的处理方法,最后成功地使H020U井连通第一靶点V020井和第二靶点V020U井。  相似文献