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

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

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.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

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

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

Distribution characteristics of delta reservoirs reshaped by bottom currents: A case study from the second member of the Yinggehai Formation in the DF1-1 gas field,Yinggehai Basin,South China Sea

The Dongfang1-1 gas field (DF1-1) in the Yinggehai Basin is currently the largest offshore self-developed gas field in China and is rich in oil and gas resources. The second member of the Pliocene Yinggehai Formation (YGHF) is the main gas-producing formation and is composed of various sedimentary types; however, a clear understanding of the sedimentary types and development patterns is lacking. Here, typical lithofacies, logging facies and seismic facies types and characteristics of the YGHF are identified based on high-precision 3D seismic data combined with drilling, logging, analysis and testing data. Based on 3D seismic interpretation and attribute analysis, the origin of high-amplitude reflections is clarified, and the main types and evolution characteristics of sedimentary facies are identified. Taking gas formation upper II (IIU) as an example, the plane distribution of the delta front and bottom current channel is determined; finally, a comprehensive sedimentary model of the YGHF second member is established. This second member is a shallowly buried “bright spot” gas reservoir with weak compaction. The velocity of sandstone is slightly lower than that of mudstone, and the reflection has medium amplitude when there is no gas. The velocity of sandstone decreases considerably after gas accumulation, resulting in an increase in the wave impedance difference and high-amplitude (bright spot) reflection between sandstone and mudstone; the range of high amplitudes is consistent with that of gas-bearing traps. The distribution of gas reservoirs is obviously controlled by dome-shaped diapir structural traps, and diapir faults are channels through which natural gas from underlying Miocene source rocks can enter traps. The study area is a delta front deposit developed on a shallow sea shelf. The lithologies of the reservoir are mainly composed of very fine sand and coarse silt, and a variety of sedimentary structural types reflect a shallow sea delta environment; upward thickening funnel type, strong toothed bell type and toothed funnel type logging facies are developed. In total, 4 stages of delta front sand bodies (corresponding to progradational reflection seismic facies) derived from the Red River and Blue River in Vietnam have developed in the second member of the YGHF; these sand bodies are dated to 1.5 Ma and correspond to four gas formations. During sedimentation, many bottom current channels (corresponding to channel fill seismic facies) formed, which interacted with the superposed progradational reflections. When the provenance supply was strong in the northwest, the area was dominated by a large set of delta front deposits. In the period of relative sea level rise, surface bottom currents parallel to the coastline were dominant, and undercutting erosion was obvious, forming multistage superimposed erosion troughs. Three large bottom current channels that developed in the late sedimentary period of gas formation IIU are the most typical.     

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模型能合理地反映长江口附近海域台风浪的分布.     

非平稳条件下北京市最大月降水量频率特征分析

为探究气候变化下极端降水的频率变化特征,基于北京市22个雨量站实测月降水量数据,以时间为协变量构建平稳和非平稳GEV模型,对北京市最大月降水量序列(极值降水序列)进行模拟和频率分析,并采用Bootstrap方法对频率分析结果的不确定性进行评价.结果表明:所有极值降水序列的最优概率分布模型均为非平稳GEV模型,该模型能够...     

应用重力场模型GL1500E和LLR天平动参数估算月核大小及其组成

针对月核大小及密度估计问题,利用高分辨率重力场模型GL1500E的二阶位系数和LLR(lunar laser ranging)天平动参数,考虑月核分层为外核和内核的情况,结合非线性粒子群优化算法,对月核大小和组成进行了估计。大批量的统计结果表明,大概率分布的外核半径rc1约为469 km,内核半径rc2约为303 km,外核密度ρc1约为4613 kg/m3,内核密度ρc2约为7004 kg/m3,月幔密度ρm约为3340 kg/m3。月幔密度非常接近地质研究结果3360 kg/m3,内外核半径与近期其他研究结果相近,估算的月核大小与密度组成具有一定的参考价值。若月核由纯铁和硫化亚铁组成,则研究表明,内核大部分由纯铁组成,外核大部分由硫化亚铁组成。