土水特征曲线(SWCC)的滞回特性模拟研究

为了模拟非饱和土的水力、力学特性在多次降雨、蒸发过程中的变化特性,本文以传统域模型的基本原理为基础,推到得到了一个能够模拟SWCC滞后性的计算模型,该模型计算方法简单,易于通过程序实现,且通过与试验数据和其他模型的计算结果对比,验证了该模型的合理性,同时,这个模型可为研究在复杂吸力变化状态下非饱和土的渗流特性、强度特性以及本构模型的研究打下一定的基础。     

The Coarse-Grained Lacustrine Slope Apron Deposits in the Moliqing Area, Yitong Basin, Northeast China

Overlapping gravity accumulation bodies were formed on the northwestern steep slope of the Shuangyang Formation in the Moliqing fault depression of northeast China. This study analyzed in detail the spatial distribution of the lithofacies and lithofacies associations of these accumulation bodies based on more than 600 m of core sections, and summarized 12 major types of lithofacies and three types of lithofacies associations: (1) the proximal zone consists of gravelly debris flows dominated by alluvial channel conglomerates; (2) the middle zone is dominated by various gravity flow deposits and traction flow deposits; and (3) the distal zone is dominated by mudstones with intercalations of sandy debris and turbidites. Combining with the grain size cumulative probability curves analysis, we determined the transformation of debris flows to sandy debris flows and to turbidity currents in the slope zone of the basin margin, and further proposed a lacustrine slope apron model that is characterized by (1) an inconstant multiple source (line source), (2) an alternation of gravity flow deposits and traction flow deposits dominated by periodical changes in a source flood flow system, and (3) the transformation of sandy debris flow deposits into distal turbidity current deposits. This sedimentary model may be applicable to other fault depressions for predicting reservoir distribution.     

以钻代槽勘查方法研究与应用

为解决槽、井探等传统地表揭露手段自身局限性问题及减少对环境的破坏,提出以钻代槽勘查方法研究的技术路线,总结归纳出浅钻技术方法与钻孔布孔方案,同时初步探讨了基于钻孔数据的三维地质建模方法。最后将以钻代槽方法成果在西藏等地进行了推广应用,验证了可行性和有效性。     

Non-symmetrical levee breaching processes in a channel bend due to overtopping

A series of experiments were done to reveal the overtopping breaching process of non-cohesive and cohesive levees in a U-bend flume. The flood hydrograph and breaching geometry were measured and analyzed in detail. The results show that the levee breaching processes can be briefly divided into four stages: slope erosion, longitudinal headward gully-cutting, lateral erosion, and relative stabilization. For non-cohesive levees, non-symmetrical lateral development of the breach occurs throughout the four stages, and the final non-symmetrical coefficient is approximately 2.2–2.6. Larger flow discharge or higher water level can accelerate the breaching process, while coarser sands tend to accelerate the process initially but depress the process at the end. The fluvial erosion rate of a non-cohesive breach shows a power-function relation with the excess wall shear stress. For cohesive levees, a plateau forms in the breach partially blocking the flow in the first two stages. The breach flow is approximately perpendicular to the levee body, and, thus, the erosion rates of the two breach sides are almost the same. Non-symmetrical lateral development mainly occurs in the third stage when the deep gully forms. The final non-symmetrical coefficient is approximately 2.7–3.3. It is expected that these findings can provide a valuable experimental dataset and a theoretical basis for breach closure and flood alleviation.     

桶形基础沉贯室内模型试验研究

桶形基础采用负压沉贯法施工,其贯入阻力与压桩,打桩等施工方法显著不同。本文介绍了在粉土和粉质粘土地层中的一组模型试验情况,认为负压法施工可显著降低粉土的抗贯阻力,但在粉质粘土中减阻效果不明显。     

IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

Tests and applications of an approach to absorbing reflected waves towards incident boundary

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.     

长江口附近海域台风浪的数值模拟--以鹿沙台风和森拉克台风为例

利用SWAN波浪模型计算长江口附近海域的台风浪,鉴于长江河口岸界和地形复杂,拟采用曲线网格.为证实曲线网格下的SWAN模型对于复杂地形的有效性,首先选用美国特拉华大学波浪水池实验资料对SWAN模型进行检验,结果表明利用曲线网格能不过多增加计算量而提高关键区域的计算精度.以0215号鹿沙台风和0216号森拉克台风为例,将SWAN模型应用到长江口附近海域,进行台风浪的数值模拟.通过浮标测站实测资料验证,表明有效波高计算值与实测值符合良好.通过综合分析模型计算的波浪场,说明SWAN模型能合理地反映长江口附近海域台风浪的分布.     

资源潜力评价中典型矿床与区域矿产编图思路讨论

根据新一轮全国矿产资源潜力评价工作需要,首先按矿产预测类型选择典型矿床编制典型矿床成矿要素图—典型矿床预测要素图—典型矿床成矿模式图和预测模型图,然后再按预测研究工作区范围编制区域地质建造构造矿产图—区域成矿要素图—区域预测要素图—区域成矿模式图和预测模型图.讨论了有关图件的编制思路,并提出了图中突出显示成矿要素和预测要素的方法.     

GIS工作流系统的长事务处理研究

分析了Sagas事务模型及空间数据版本树技术,设计了一种灵活的长事务补偿控制方案,为集成GIS,MIS及OA的工作流系统提供了有力的事务处理技术支持。