About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Helmholtz算子逆的多方向分裂

Fourier有限差分(FFD)方法是一个重要的波场深度延拓方法，它需要用到包括Taylor展开、连分式晨开、迭代展开等一系列数学展开，将FFD方法从二维推广到三维时，需要解决二维Helmholtz算子的求逆问题。本文利用多方向分裂方法导出了二维Helmholtz算子逆的近似公式，可以应用于三维波场深度延拓问题的研究。     

时域瞬变场电磁场有限差分法

模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。     

西安市黑河水库左坝肩渗漏三维数值模拟及坝肩坡体稳定性分析

直接针对西安市黑河水库左坝肩渗漏这一重大工程实际问题展开研究与计算。应用三维有限差分方法对左坝肩单薄山梁的初始渗流场、蓄水后的渗流场进行了数值模拟 ,分析渗漏方式、计算相应的渗漏量 ,并通过计算说明了断层是集中渗漏通道。而后利用三维数值计算提供的水位值 ,应用极限平衡分析法 ,对左坝肩岩质斜坡稳定性进行计算与分析 ,计算结果对实际工程的决策与治理提供重要的依据。     

网格剖分及其精度和计算量分析

网格剖分程度直接影响着地震波正演数值模拟的计算精度及其计算量.以均匀介质模型为例,分析不同网格大小对波场模拟精度和计算量的影响,得出精细化网格剖分是实现高精度地震波正演模拟的有效方法,然而其计算量较大.以均匀倾斜介质模型为例,探讨倾斜地层网格剖分问题,数值实例分析不同震源频率对不同网格剖分方案引起的波场传播精度的影响,...     

Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

Seismic,petrologic, and geodetic analyses of the 1999 dome-forming eruption of Guagua Pichincha volcano,Ecuador

Guagua Pichincha, located 14 km west of Quito, Ecuador, is a stratovolcano bisected by a horseshoe-shaped caldera. In 1999, after some months of phreatic activity, Guagua Pichincha entered into an eruptive period characterized by the extrusion of several dacitic domes, vulcanian eruptions, and pyroclastic flows. We estimated the three-dimensional (3-D) P-wave velocity structure beneath Guagua Pichincha using a tomographic inversion method based on finite-difference calculations of first-arrival times. Hypocenters of volcano-tectonic (VT) earthquakes and long-period (LP) events were relocated using the 3-D P-wave velocity model. A low-velocity anomaly exists beneath the caldera and may represent an active volcanic conduit. Petrologic analysis of eruptive products indicates a magma storage region beneath the caldera, having a vertical extent of 7–8 km with the upper boundary at about sea level. This zone coincides with the source region of deeper VT earthquakes, indicating that a primary magma body exists in this region. LP swarms occurred in a cyclic pattern synchronous with ground deformation during magma extrusions. The correlation between seismicity and ground deformation suggests that both respond to pressure changes caused by the cyclic eruptive behavior of lava domes.