A FINITE ELEMENT SHALLOW WATER CIRCULATION MODEL FOR TOLO HARBOUR, HONG KONG

Numerical solution of two-dimensional shallow water equations has been achieved by the finite element method. In this work, the governing equations were first converted into different equations suitable for explicit, forward marching, integration in the time domain. A second order Taylor's series expansion was used for the continuity equation and an "upwind" differencing was incorporated in the nonlinear convective terms in the momentum equations.The model was applied to Tolo Harbour in Hong Kong. Velocity vector plots and surface elevations obtained from the calculations are reasonably consistent with the expected tidal response of the systems.     

Full 2D Hydrodynamic Modelling of Rainfall-induced Flash Floods

Mountain catchments are prone to flash flooding due to heavy rainfall. Enhanced understanding of the generation and evolution processes of flash floods is essential for effective flood risk management. However, traditional distributed hydrological models based on kinematic and diffusion wave approximations ignore certain physical mechanisms of flash floods and thus bear excessive uncertainty. Here a hydrodynamic model is presented for flash floods based on the full two-dimensional shallow water equations incorporating rainfall and infiltration. Laboratory experiments of overland flows were modelled to illustrate the capability of the model. Then the model was applied to resolve two observed flash floods of distinct magnitudes in the Lengkou catchment in Shanxi Province, China. The present model is shown to be able to reproduce the flood flows fairly well compared to the observed data. The spatial distribution of rainfall is shown to be crucial for the modelling of flash floods. Sensitivity analyses of the model parameters reveal that the stage and discharge hydrographs are more sensitive to the Manning roughness and initial water content in the catchment than to the Green-Ampt head. Most notably, as the flash flood augments due to heavier rainfall, the modelling results agree with observed data better, which clearly characterizes the paramount role of rainfall in dictating the floods. From practical perspectives, the proposed model is more appropriate for modelling large flash floods.     

Dynamic process simulation with a Savage-Hutter type model for the intrusion of landslide into river

Natural damming of rivers by mass movements is a very common and potentially dangerous phenomena which has been documented all over the world. In this paper, a two-layer model of Savage-Hutter type is presented to simulate the dynamic procedure for the intrusion of landslide into rivers. The two-layer shallow water system is derived by depth averaging the incompressible Navier-Stokes equations with the hydrostatic assumption. A high order accuracy scheme based on the finite volume method is proposed to solve the presented model equations. Several numerical tests are performed to verify the realiability and feasibility of the proposed model. The numerical results indicate that the proposed method can be competent for simulating the dynamic process of landslide intrusion into the river. The interaction effect between both layers has a significant impact on the landslide movement, water fluctuation and wave propagation.     

Two-dimensional numerical simulation of flow around three-stranded rope

Three-stranded rope is widely used in fishing gear and mooring system. Results of numerical simulation are presented for flow around a three-stranded rope in uniform flow. The simulation was carried out to study the hydrodynamic characteristics of pressure and velocity fields of steady incompressible laminar and turbulent wakes behind a three-stranded rope. A three-cylinder configuration and single circular cylinder configuration are used to model the three-stranded rope in the two-dimensional simulation. The governing equations, Navier-Stokes equations, are solved by using two-dimensional finite volume method. The turbulence flow is simulated using Standard κ-ε model and Shear-Stress Transport κ-ω(SST) model. The drag of the three-cylinder model and single cylinder model is calculated for different Reynolds numbers by using control volume analysis method. The pressure coefficient is also calculated for the turbulent model and laminar model based on the control surface method. From the comparison of the drag coefficient and the pressure of the single cylinder and three-cylinder models, it is found that the drag coefficients of the three-cylinder model are generally 1.3–1.5 times those of the single circular cylinder for different Reynolds numbers. Comparing the numerical results with water tank test data, the results of the three-cylinder model are closer to the experiment results than the single cylinder model results.     

GENERALIZED WAVE EQUATION FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL TIDAL WAVES

The study of tidal circulation has a long history . The numerical simulation of tidal flow has been developed greatly with the development of computer techniques in the past two decades. The generalized wave equation finite-element method is a relatively new numerical model for studying shallow water flow . This method was used to simulate tidal waves of the Gulf of St. Lawrence in Canada . The very good agreement of the numerical results with the field data indicated that the model is an effective and promising numerical method for solving two-dimensional tidal wave problems .     

Numerical simulation of mechanical compaction of deepwater shallow sediments

To study the compaction law and overpressure evolution in deepwater shallow sediments, a large-strain compaction model that considers material nonlinearity and moving boundary is formulated. The model considers the dependence of permeability and material properties on void ratio. The modified Cam-Clay model is selected as the constitutive relations of the sediments, and the deactivation/reactivation method is used to capture the moving top surface during the deposition process. A one-dimensional model is used to study the compaction law of the shallow sediments. Results show that the settlement of the shallow sediments is large under their own weight during compaction. The void ratio decreases strictly with burial depth and decreases more quickly near the seafloor than in the deeper layers. The generation of abnormal pressure in the shallow flow sands is closely related to the compaction law of shallow sediments. The two main factors that affect the generation of overpressure in the sands are deposition rate and permeability of overlying clay sediments. Overpressure increases with an increase in deposition rate and a decrease in the permeability of the overlying clay sediment. Moreover, an upper limit for the overpressure exists. A two-dimensional model is used to study the differential compaction of the shallow sediments. The pore pressure will still increase due to the inflow of the pore fluid from the neighboring clay sediment even though the deposition process is interrupted.     

西北内陆盆地浅层地下水温度场及动态类型: 以酒泉东盆地为例

温度是地下水的固有属性，地下水温度场和动态特征是地下水流系统的外在表现。为揭示地下水开采等人类活动影响下西北内陆盆地浅层地下水温度场特征与地下水流系统的关系，基于多点位、长序列、高精度的地下水温度监测数据，在酒泉东盆地开展了地下水温度场及动态特征研究。结果表明:酒泉东盆地浅层地下水温度9.33~20.77℃不等，平均水温为13.54℃，自地下水补给区至排泄区，沿地下水径流方向，浅层地下水温度逐渐升高；循环深度相近的不同地下水流系统对比表明，浅层地下水温度与地下水动力条件呈负相关，地表水入渗补给大、水动力条件强的水流系统地下水平均温度低，入渗补给小、水动力条件弱的地下水平均温度高；浅层地下水温度动态受自然地下水循环和地下水开采等人类活动共同影响，从山前地下水补给区到中游绿洲区再到下游排泄区，浅部地下水温度动态可划分为4种基本类型，依次分别为河流补给型、水温稳定型、开采相关型、正弦波动型。      

隔档式构造区岩溶地下水流系统多级次嵌套结构的水化学识别

川东隔档式构造区岩溶地下水流系统表现出多级次嵌套结构的特征, 之前主要是定性的水文地质条件分析, 缺乏来自水化学方面的量化依据。利用41个浅层和19个深层岩溶水样的水化学资料开展分析, 发现浅层岩溶水为HCO₃-Ca·Mg型,ρ(Mg2+)较低, 深层岩溶水为SO₄-Ca·Mg型且ρ(Mg2+)明显偏高, 说明地下水在岩溶含水层中的滞留时间与ρ(Mg2+)具有正变关系。对典型剖面地下水流系统的分析表明, 优先采用ρ(Mg2+)来评价地下水滞留时间, 将地下水排泄点ρ(Mg2+)小于20, 50 mg/L分别作为划分局部-中间、中间-区域水流系统的依据。浅层岩溶水受地貌作用控制明显, 其中浅切沟谷泉点为局部地下水流系统的排泄点, 深切沟谷泉点一般属于中间或区域流动系统的排泄点。ρ(Mg2+)反映了泉水循环的滞留时间, 也能够反映钻孔所揭露的深循环特征。这种水化学识别方法可为相似岩溶区识别地下水流系统的多级次嵌套结构提供参考。      

A quasi single-phase model for debris flows and its comparison with a two-phase model

A depth-averaged quasi single-phase mixture model is proposed for debris flows over inclined bed slopes based on the shallow water hydrosediment-morphodynamic theory with multi grain sizes. The stresses due to fluctuations are incorporated based on analogy to turbulent flows, as estimated using the depth-averaged k-? turbulence model and a modification component. A fully conservative numerical algorithm, using wellbalanced slope limited centred scheme, is deployed to solve the governing equations. The present quasi single-phase model using four closure relationships for the bed shear stresses is evaluated against USGS experimental debris flow and compared with traditional quasi single-phase models and a recent physically enhanced two-phase model. It is found that the present quasi single-phase model performs much better than the traditional models, and is attractive in terms of computational cost while the two-phase model performs even better appreciably.     

A COMBINED REFRACTION-DIFFRACTION-DISSIPATION MODEL OF WAVE PROPAGATION

A numerical model based on the mild-slope equation of water wave propagation over complicated bathymetry,taking into account the combined effects of refraction,diffraction and dissipation due to wavebreaking is presented.Wave breaking is simulated by modifying the wave height probability density func-tion and the wave energy dissipation mechanism is parameterized according to that of the hydraulic jumpformulation.Solutions of the wave height,phase function,and the wave direction at every grid point areobtained by finite difference approximation of the governing equations,using Gauss-Seidel Iterative Method(GSIM)row by row.Its computational convenience allows it to be applied to large coast regions tostudy the wave transformation problem.Several case studies have been made and the results compare verywell with the experiment data and other model solutions.The capability and utility of the model forreal coast areas are illustrated by application to a shallow bay of northeast Australia.     

A numerical model for tidal currents in the Changjiang River estuary

A generalized finite difference scheme is presented for the solution of two-dimensioned tidal current equations (averaged over the depth of water). This method has been used to simulate the tidal currents in the Changjiang River estuary, with computed results being in good agreement with the field observations. The two-dimensional domain consists of triangular elements with acute angles. First- and second-order partial derivatives are defined at nodal points. In addition, the basic tidal current equations contain an extra term to account for horizontal eddy viscosity. The method has the great advantage of allowing a wide flexibility in configuration simulation, thereby exhibiting the capability to arbitrarily subdivide the region of interest as well as to provide a more accurate boundary shape, and what is more, it offers savings in computing costs. This paper was published in Chinese in the Ocean. et Limn. Sinica.16(1):18–27, 1985     

Model of wind-generated ambient noise in stratified shallow water

In order to build a rapid ocean ambient noise model adapted for a stratified shallow water, a hybrid model of normal mode method (for far field) and ray method (for near field) is suggested which combines the advantages of both methods. Since the near field of wind-generated noise is not sensitive to the sound speed profile, the sound speed profile is regarded as a constant; which makes the model rapid and accurate. The simulation results are in agreement with those of the wave model.     

使用L0范数代价函数的二维相位快速解缠方法

针对最小范数解缠方法的代价函数模型以及迭代效率问题,提出一种符合L0范数准则的高效的二维相位全局最优解缠方法。在对最小范数代价函数模型特征进行分析的基础上,给出一种符合L0范数准则的代价函数,满足在相位不连续边界的切方向上相对法方向具有更强的约束作用,保证解缠迭代处理时连续相位按照代价函数的约束条件进行趋近的同时,不连续边界不会被破坏。为了解决线性方程中低频误差收敛慢的问题,根据最小范数方法独立性的特点,采用数据分块处理方法,注重高频信息以提高迭代处理效率,将低频信息处理转移到数据块间偏移量计算中。实验对比分析表明,新解缠方法在可靠性和效率方面都能够达到较好的水平。     

Algal genetic engineering in China: Progress and cooperation

INTRODUCTIONTheWellknownthooriesofPritChard(l952)andHansenandRattray(l965)destribeestUarine~ationasfashenedseawatCrflowingoutinanearsurfacelayerandanwaterflowinlandwardbeow.WhenthelateralvariationofbathyInetryistakenintoarmunt,theverticalshatificationrnaybereplatal,orpadiallyreplaed,byalateralstIatfficatfonalongwithalateralvariatinnofveodtyfiekl,whichrnayshowalandwardflowofsaltwaterinthedeepchannelandanWhrdflowoffasherWhterontheshoaIs.fumodeofbarotaedationWhsfustsuggestalbyFischer(l9…     

利用层次分析法评价东营市浅层地热能适宜性

利用定性、半定量相结合的层次分析法(AHP),建立地埋管换热系统适宜性递阶层次结构模型,在野外试验、室内测试等基础数据的基础上,选取地下水埋藏深度、砂层厚度、成孔难易程度、地面沉降、热扩散率、比热容、综合热导率等7个评价指标,将东营市浅层地热能开发利用划分为适宜区、较适宜区。这是东营地区第一次利用实测数据开展浅层地热能适宜性评价,为今后的合理开发利用提供了科学依据。     

Underwater topography detection of Taiwan Shoal with SAR images

Under suitable conditions of tidal current and wind, underwater topography can be detected by synthetic aperture radar (SAR) indirectly. Underwater topography SAR imaging includes three physical processes: radar ocean surface backscattering, the modulation of sea surface short wave spectrum by the variations in sea surface currents, and the modulation of sea surface currents by the underwater topography. The first process is described usually by Bragg scattering theory because the incident angle of SAR is always between 20°–70°. The second process is described by the action balance equation. The third process is described by an ocean hydrodynamic model. Based on the SAR imaging mechanism for underwater topography, an underwater topography SAR detection model and a simplified method for its calculation are introduced. In the detection model, a two-dimensional hydrodynamic model — the shallow water model is used to describe the motion of tidal current. Due to the difficulty of determining the expression of SAR backscattering cross section in which some terms can not be determined, the backscattering cross section of SAR image used in the underwater topography SAR detection is pro-processed by the simulated SAR image of the coarse-grid water depth to simplify the calculation. Taiwan Shoal, located at the southwest outlet of Taiwan Strait, is selected as an evaluation area for this technique due to the occurrence of hundreds of sand waves. The underwater topography of Taiwan Shoal was detected by two scenes of ERS-2 SAR images which were acquired on 9 January 2000 and 6 June 2004. The detection results are compared with in situ measured water depths for three profiles. The average absolute and relative errors of the best detection result are 2.23 m and 7.5 %, respectively. These show that the detection model and the simplified method introduced in the paper is feasible.     

Environmental capacity of chemical oxygen demand in the Bohai Sea: modeling and calculation

A three-dimensional advection-diffusion model coupled with the degradation process is established for describing the transport of chemical oxygen demand (COD). Comparison of the simulated distribution of COD at the surface in the Bohai Sea in August, 2001 with field observations, shows that the model simulates the dataset reasonably well. The Laizhou Bay, Bohai Bay, and Liaodong Bay were contaminated heavily near shore. Based on the optimal discharge flux method, the Environmental Capacity (EC) and allocate...     

干热岩勘查评价指标与形成条件

干热岩具有利用率高、无污染、储量巨大、分布较广、持续稳定、安全性好等特点,被全球公认为在21世纪能够取代化石能源的一种最优质的、可再生的新型清洁能源。目前对于干热岩的成因机理还没有定论,而对于干热岩的勘探寻找、远景区的圈定以及资源评价也有不同的观点,为了有利于干热岩勘探开发,本文在综述前人研究工作基础上,总结归纳了干热岩勘查评价指标,几个重要的指标及其特征如下。第一个指标是岩石圈厚度和莫霍面埋深,岩石圈厚度较小并且莫霍面埋深比较浅,是评价干热岩远景区的一个重要指标。莫霍面埋深较浅指示深部热源（幔源热）更加接近地表,故埋深较浅（较薄的厚度）并具有上隆的特征有利于深部的热量向上传导,为干热岩的孕热环境提供良好的条件。第二个指标是居里等温面,埋深较浅的居里等温面是干热岩远景区一个重要的评价指标。居里等温面是地球内部一个非常特殊的温度（热物质）界面,它不仅能指示地下温度场的分布特征,还可指示地壳深部热能分布特征,对干热岩及地震的成因研究具有十分重要的意义;如果居里面埋深较浅,则热量传导到地表的距离比较短,深部的热流活动更容易向地表传送,不仅是有利的高温干热岩孕热环境,也有利于储存的热能快速向上传导。第三个指标为地温梯度,地温梯度较大是寻找干热岩远景区的一个重要指标;如果一个地区具有较高的地温梯度（≥ 35℃/km）,则随着深度的增加深部的地温增加较快,在相对较浅的地方就可以获得温度较高的岩石体。第四个指标是大地热流,大地热流值较大（≥ 75 mW/m2）就指示地球深部有存在高温岩石体的可能;大地热流是地温场的综合性热参数,能够准确地反映区域内的地温场特征。第五个指标是新构造运动,这也是人们寻找干热岩时容易忽略的一个重要指标。新构造运动包括火山、地震及活动断裂构造等。地震和火山是极具破坏性的自然灾害,两者的发生都表明了地球内部的热能汇聚到一定程度,从而打破了地球内部平衡状态而以地震或火山的形式把热量进行释放的地球系统行为。该指标中,如地震震级大（>3级）、震源深度浅（10~15 km）、频度大,火山活动时间新（活火山、休眠火山、中新世以来的死火山）规模大都说明地球深部存在不稳定的高热状态,易形成干热岩;如果能提前找到该区域的干热岩,可以先取出其中的热,那么地震和火山就有可能不会发生,这样可以达到取热减灾减排的作用。活动的深大断裂即能产生一部分热,也能将深部的热传输到浅部,尤其是活动性强的走滑拉张断裂,其深部具有韧性剪切特征,直接指示了深部的高温体的存在。第六个指标是高温温泉与气田等。温泉、气田的形成通常与深部的热储关系密切,一般认为地下水沿某个通道向下渗透接触到深部高温热储被加热后再沿某一通道流出地表而形成温泉;所以,温泉的出露指示了深部存在高温的岩石体（干热岩）;区域地温场异常明显,地表热泉等高温水热型地热田较密集的区域有望在深部寻找到干热岩,这也是一种就热（水热）找热（干热）的常规方法。作为固体矿产资源的干热岩,其形成具备四个必要条件:源、通、储、盖。第一个条件是要有丰富的动态热源如来源于深部地幔（幔源热）,来源于晚新生代活动的控热构造系统-活动的韧性剪切带,来源于地壳内的低速低阻体（中下地壳热）以及来源于高放射性中新生代花岗岩体（壳源热）。第二个条件是要有优良的导热通道,如壳内15~25 km低速层不仅是热源,同时具有将深部地幔热能向上传导的作用;软流圈地幔上隆时具活动性的深大断裂（深部具韧性特征、浅部具脆性特征）常常具有很好的导热功能;地壳浅表层次的脆性断裂系统往往不是干热岩的热通道,而是水热型地热能的导水、释热构造。第三个条件是要有巨大的储热岩石体,除埋深要适中（3~6 km）并具有较高温度（≥150℃）外,其规模要大（蕴含丰富的热能）,热导率大（>2 W/mK）,裂隙少（不含水或含少量不流动的水）;当然热储层可以是变质岩、岩浆岩,也可以为沉积岩。第四个条件是要有良好的保热盖层,盖层（被子）导热率低（< 2 W/mK）、厚度适中（>1 km）（具有良好的保温效果）,地温梯度高（≥40℃/km）、大地热流值高（≥70 mW/m2）（指示深部存在高温特征）是深部赋存有高温地热资源的必要条件。      

On the Lagrangian residual current and residual transport in a multiple time scale system of shallow seas

A multiple time scale perturbation method is used to discuss the Lagrangian residual current and residual transport on the basis of a weakly nonlinear dynamic model of shallow seas. The governing equations for the long-term variation of zero order “apparent concentration” (which is a linear combination of salinity, temperature of seawater and the concentration of any tracer which is conservative and passive) and its mean value over tidal cycles are obtained for the system with single tidal constituent, and for the one with multi-constituents, winds and thermohaline. The equations for the two cases are in the same form and show this long-term variation resulted from the cumulative effect of residual convection and turbulent diffusion. The multiple time scale variation of current is caused by tides, winds, and the thermohaline and the nonlinear effects of the system. The derived set of governing field equations of the Lagrangian current for this multiple time scale system is also in the same form as that for a single time scale system.     

THE EFFECTS OF A SLOPING BOTTOM ON THE EQUATORIAL WAVE-INDUCED DEEP RESIDUAL CURRENTS

ImODUcnONThedeepequatorialoceanhasobvioussignilicantflowapparentlycarryinghacelsa1ongandacrosstheequator(WissCtal.,l985).RmtfloatmsurementSshoWedthattheflowishigh1yvariable(Richardsonetal.,l993).ThomPsonandKawase(l993)pro-posedthatthelargeinstantaneousandfloatvelocitiesasWellasthevariabilitysuggestthatthetracersignaIsreflCCtreCtificationoftimedependentmotionsandniinginsteadofrneanEulerianflow.TheresultSoftheirstudyonthegenerationofmeancurrentSbyperi-odicfordngintheequatorialoceaninasir…