New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation. A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.     

利用自洽方程计算法分析天山中段分量钻孔应变资料可靠性

根据天山中段8套四分量钻孔应变仪观测资料,利用自洽方程计算应变观测的信度,用信度指标检验数据的可靠性,从而进一步用自洽方程计算实际各元件灵敏度的相对大小,并用相对标定的结果,对观测数据进行校正。结果表明:①石场、巴仑台、巴音沟、雀儿沟、榆树沟和库米什台经过校正后,数据观测信度有了明显提升,其中石场和雀儿沟台的数据信度最高;②库米什台校正前数据观测信度较高,因个别“坏点”影响,校正后,观测信度反而下降;③小泉沟和新源台数据信度较低,存在从自洽到失洽的情况,分析认为这可能与系统自身存在问题有关,也有可能是前兆异常,有待进一步研究。     

基于偏微分方程的模型光滑处理在偏移成像中的应用

本文针对射线类偏移成像当中的速度模型光滑处理问题,借鉴数字图像处理当中的偏微分方程法,基于能量泛函,应用变分方法导出基于速度模型的偏微分方程实现射线类偏移成像当中的速度模型的光滑处理.由于偏微分方程法具有线性叠加特性、模型解的唯一性和局部特征保持性,因此,应用该算法可以实现基于原始速度模型空间结构的模型光滑处理.通过在原始速度模型以及光滑处理后的速度模型上计算速度的空间分布以及地震波走时、射线路径可以得出,偏微分方程法对速度模型的光滑处理能够很好地保持原始模型的空间结构,偏移成像结果也证明了该方法的实用性.

     

The mean and turbulence structure simulation of the monsoon trough boundary layer using a one-dimensional model withe-l ande-ε closures

An attempt has been made here to study the sensitivity of the mean and the turbulence structure of the monsoon trough boundary layer to the choice of the constants in the dissipation equation for two stations Delhi and Calcutta, using one-dimensional atmospheric boundary layer model withe-ε turbulence closure. An analytical discussion of the problems associated with the constants of the dissipation equation is presented. It is shown here that the choice of the constants in the dissipation equation is quite crucial and the turbulence structure is very sensitive to these constants. The modification of the dissipation equation adopted by earlier studies, that is, approximating the Tke generation (due to shear and buoyancy production) in theε-equation by max (shear production, shear + buoyancy production), can be avoided by a suitable choice of the constants suggested here. The observed turbulence structure is better simulated with these constants. The turbulence structure simulation with the constants recommended by Aupoixet al (1989) (which are interactive in time) for the monsoon region is shown to be qualitatively similar to the simulation obtained with the constants suggested here, thus implying that no universal constants exist to regulate dissipation rate. Simulations of the mean structure show little sensitivity to the type of the closure parameterization betweene-l ande-ε closures. However the turbulence structure simulation withe-ε. closure is far better compared to thee-l model simulations. The model simulations of temperature profiles compare quite well with the observations whenever the boundary layer is well mixed (neutral) or unstable. However the models are not able to simulate the nocturnal boundary layer (stable) temperature profiles. Moisture profiles are simulated reasonably better. With one-dimensional models, capturing observed wind variations is not up to the mark.     

Pore characterization of organic-rich Lower Cambrian shales in Qiannan Depression of Guizhou Province,Southwestern China

Nine organic-rich shale samples of Lower Cambrian black shales were collected from a recently drilled well in the Qiannan Depression, Guizhou Province where they are widely distributed with shallower burial depth than in Sichuan Basin, and their geochemistry and pore characterization were investigated. The results show that the Lower Cambrian shales in Qiannan Depression are organic rich with TOC content ranging from 2.81% to 12.9%, thermally overmature with equivalent vitrinite reflectance values in the range of 2.92–3.25%, and clay contents are high and range from 32.4% to 53.2%. The samples have a total helium porosity ranging from 2.46% to 4.13% and total surface area in the range of 9.08–37.19 m²/g. The estimated porosity in organic matters (defined as the ratio of organic pores to the volume of total organic matters) based on the plot of TOC vs helium porosity is about 10% for the Lower Cambrian shales in Qiannan Depression and is far lower than that of the Lower Silurian shales (36%) in and around Sichan Basin. This indicates that either the organic pores in the Lower Cambrian shale samples have been more severely compacted than or they did not develop organic pores as abundantly as the Lower Silurian shales. Our studies also reveal that the micropore volumes determined by Dubinin–Radushkevich (DR) equation is usually overestimated and this overestimation is closely related to the non-micropore surface area of shales (i.e. the surface area of meso- and macro-pores). However, the modified BET equation can remove this overestimation and be conveniently used to evaluate the micropore volumes/surface area and the non-micropore surface areas of micropore-rich shales.     

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.     

Static bay beach concept for scientists and engineers: A review

Headland-bay beach (HBB) is one of the most prominent physiographic features on the oceanic margin of many countries in the world. Under the influence of a predominant swell, its curved periphery in natural environment may reach static equilibrium and remains stable without sediment supply from updrift and/or a riverain source within its own embayment. Coastal scientists and engineers have attempted to develop mathematical expressions to quantify this ideal bay shape since the 1940s. As the scenario with depleting sediment supply has become a common reality on many parts of the world coastline in more recent time, some coastal engineers have advocated a rational approach to mimicking the static bay shape found in nature in order to mitigate beach erosion as well as for coastal management. Nowadays, many useful applications have emerged since the publication of the parabolic bay shape equation (PBSE) developed for static equilibrium planform (SEP) in late 1980s. The advance in modern computer technologies and international collaboration has further facilitated the exchange of knowledge and applications of this static bay beach concept (SBBC).     

On the equation of state for the Λ field

The recently detected accelerated expansion of the Universe is related to the existence of a new type of matter called the Λ field or quintessence. Constraints were obtained on its equation of state from the absence of clustering of this matter on scales much smaller than the cosmological horizon. The question of how these constraints affect the possibility of fitting the accelerated expansion by such cosmological models as the Chaplygin gas model is discussed.     

Harmonic wave propagation through viscoelastic heterogeneous media exhibiting mild stochasticity - I. Fundamental solutions

This work examines the propagation of time harmonic, horizontally polarized shear waves through a naturally occurring heterogeneous medium that exhibits viscous behaviour as well as random fluctuations of its elastic modulus about a mean value. As a first step, the governing equation, which is a heterogeneous Helmholtz equation, is solved using algebraic transformations and the relevant Green's function is obtained for two sets of boundary conditions, one corresponding to a finite depth layer and the other to an infinite layer. Viscous material behaviour is introduced by considering the depth-dependent elastic modulus to be a complex quantity. Subsequently, material stochasticity in the medium is handled through the perturbation approach by assuming that the elastic modulus has a small random fluctuation about its mean value. The final results are closed-form expressions for the mean value and covariance matrix of both the wave speed profile in the medium and the corresponding Green's function. In Part II, (Soil Dynam. Earth. Engng, 1996,15, 129-39), two examples concerning seismic wave propagation in soft topsoil and in sandstone serve to illustrate the methodology and comparisons are made with Monte Carlo simulations.