基于BP人工神经网络的水体遥感测深方法研究

利用Landsat7ETM＋遥感图像反射率和实测水深值之间的相关性,建立了动量BP人工神经网络水深反演模型,并对长江口南港河段水深进行了反演.结果表明：具有较强非线性映射能力的动量BP神经网络模型能较好地反演出长江口南港河段的水深分布情况;由于受长江口水体高含沙量的影响,模型对小于5 m的水深值反演精度较高,而对大于10 m的水深值反演精度较低.     

Stability of sandy slopes under seepage conditions

Stability against shallow mass sliding in saturated sandy slopes under seepage depends on the flow direction and hydraulic gradient, particularly near the ground surface. Two modes of instability i.e., Coulomb sliding and liquefaction have been studied and the critical flow directions discussed. The utility of the numerical approach in solving complex flow problems with irregular boundaries and surface topography is demonstrated by means of two slope examples with different internal drainage conditions. The numerical results for the seepage gradients at different points are compared with those predicted by the simple expression derived in this study, and the corresponding effects on the stability are evaluated.     

长江口水域多光谱遥感水深反演模型研究

利用Landsat-7 ETM 遥感影像反射率和实测水深值之间的相关性可以探测水深。该文介绍单波段、双波段比值和多波段3种线性回归模型以及动量BP人工神经网络水深反演模型。选择长江口北港河道上段作为研究区,利用上述模型,分两种情况进行水深反演:一是以河道全部历史样本建模;二是将河道按自然水深划分为浅水区和深水区分别建模。结果表明:神经网络模型预测精度高于线性回归模型;水深分区后线性回归和神经网络模型预测误差均有所减小。     

固体介质中弹性波传播机制的细胞自动机有限深势阱模型

分析了采用细胞自动机研究波动问题的建模方法，针对一维、均匀、各向同性固体介质中弹性纵波的微观机制，借用一经典弹簧振子模型、细胞自动机格子气模型，以及量子力学中的无限深势阱模型，建立了一个细胞自动机有限深势阱模型，从量子力学角度出发，基于介观物理和纳米概念，以微观精子的德布罗意假设为基础，利用薛定谔方程，讨论了该模型中粒子（分子组）的振动速度与粒子物质波波速之间的联系，给出了模型中的波动方程，得出ζ＝Vp（ζ为粒子振动速度，Vp为物质波纵波波速）。同时还讨论了模型中粒子的大小和能量传递问题，引入引力场，得出了能量及引力势的量子化条件，另外，对声波速度、格子气粒子振动速度和本文模型中分子组振动速度进行了比较；还对本文模型中的粒子能量分布作了分析。     

A time-space based approach for mapping rainfall-induced shallow landslide hazard

A rainfall-induced shallow landslide is a major hazard in mountainous terrain, but a time-space based approach is still an unsettled issue for mapping rainfall-induced shallow landslide hazards. Rain induces a rise of the groundwater level and an increase in pore water pressure that results in slope failures. In this study, an integrated infinite slope analysis model has been developed to evaluate the influence of infiltration on surficial stability of slopes by the limit equilibrium method. Based on this new integrated infinite slope analysis model, a time-space based approach has been implemented to map the distributed landslide hazard in a GIS (Geographic Information Systems) and to evaluate the shallow slope failure induced by a particular rainfall event that accounts for the rainfall intensity and duration. The case study results in a comprehensive time-space landslide hazard map that illustrates the change of the safety factor and the depth of the wetting front over time.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

Dynamic soil–structure interaction analysis via coupled finite-element–boundary-element method

In this paper, a study on the transient response of an elastic structure embedded in a homogeneous, isotropic and linearly elastic half-plane is presented. Transient dynamic and seismic forces are considered in the analysis. The numerical method employed is the coupled Finite-Element–Boundary-Element technique (FE–BE). The finite element method (FEM) is used for discretization of the near field and the boundary element method (BEM) is employed to model the semi-infinite far field. These two methods are coupled through equilibrium and compatibility conditions at the soil–structure interface. Effects of non-zero initial conditions due to the pre-dynamic loads and/or self-weight of the structure are included in the transient boundary element formulation. Hence, it is possible to analyse practical cases (such as dam–foundation systems) involving initial conditions due to the pre-seismic loads such as water pressure and self-weight of the dam. As an application of the proposed formulation, a gravity dam has been analysed and the results for different foundation stiffness are presented. The results of the analysis indicate the importance of including the foundation stiffness and thus the dam–foundation interaction.     

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

Seismic analysis of underground tunnels by the 2.5D finite/infinite element approach

A procedure for the seismic analysis of underground tunnels using recorded free-field earthquakes based on the 2.5D finite/infinite element approach is presented. The near and far fields of the half space are modeled by finite and infinite elements, respectively. Using the 1D wave theory, the nodal force and displacement on the near-field boundary are computed for each spectral frequency of the earthquake. Then, equivalent seismic forces are computed for the near-field boundary for the earthquake spectrum. By assuming the soil-tunnel system to be uniform along the tunnel axis, the 2.5D approach can account for the wave transmission along the tunnel axis, which reduces to the 2D case for infinite transmission velocity. The horizontal and vertical components of the 1999 Chi-Chi Earthquake (TCU068) are adopted as the free-field motions in the numerical analysis. The maximal stresses and distribution patterns of the tunnel section under the P- and SV-waves are thoroughly studied by the 2.5D and 2D approaches, which should prove useful to the design of underground tunnels.     

无限水深聚焦波完全非线性数值模拟

基于势流理论提出一种新的高阶边界元方法对无限水深的聚焦波浪进行完全非线性数值模拟.自由水面满足完全非线性边界条件,模拟波浪的非线性效果可以达到更高阶.利用镜像原理,建立一种全新的格林函数应用到无限水深的数值波浪水槽中,以致于两无限深水槽侧壁的积分可以被排除.为了产生相应的入射波和吸收出流波浪,一个由点源组成的造波装置被布置于计算域内,同时人工阻尼层被用来吸引出流波浪,由波浪聚焦的方法得到极限波浪.通过开展线性和完全非线性聚焦波浪的数值实验及与理论解对比,验证本数值模型可以用来模拟无限深水域的极限波浪,且在出流边界没有反射.