Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

Analytical study of ground motion caused by seismic wave propagation across faulted rock masses

Studying seismic wave propagation across rock masses and the induced ground motion is an important topic, which receives considerable attention in design and construction of underground cavern/tunnel constructions and mining activities. The current study investigates wave propagation across a rock mass with one fault and the induced ground motion using a recursive approach. The rocks beside the fault are assumed as viscoelastic media with seismic quality factors, Q_p and Q_s. Two kinds of interactions between stress waves and a discontinuity and between stress waves and a free surface are analyzed, respectively. As the result of the wave superposition, the mathematical expressions for induced ground vibration are deduced. The proposed approach is then compared with the existing analysis for special cases. Finally, parametric studies are carried out, which includes the influences of fault stiffness, incident angle, and frequency of incident waves on the peak particle velocities of the ground motions.     

长江上游洪涝暴雨环流形势中期演变趋势分析

本文通过对1980年以来长江上游10次洪涝过程的500hpa环流形势分析,并用车比雪夫多项式对500hpa高度场进行展开,研究了各类洪涝过程的中期演变特征,发现代表不同洪涝过程环流形势变化的不同的车氏系数值都有明显的特征变化,结合欧洲中心中期数值预告96小时500hpa高度场的车氏系数分析,可为洪涝暴雨的中期预报提供一种新的思路。     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Precipitation patterns in flood season over China associated with the El Niño/Southern Osciliation

The precipitation patterns in flood season over China associated with the El Niño/Southern Oscillation (ENSO) are investigated, especially in the eastern China, using the rather long period rainfall data in this century. The results show that there were remarkable differences between the precipitation patterns in flood seasons of ENSO warm phase (El Niño year) and cold phase (La Niña year), as well as between the patterns in El Niño years and their following years. The most parts of China received below normal rainfall in flood season of the onset years of El Niño events, but the coastal area of Southeast China received above normal amounts. Comparatively, the most parts of China received above normal rainfall in flood season of the following years of El Niño events, but the eastern part of the reaches among the Huanghe (Yellow) River, the Huaihe River and the Haihe River, and the Northeast China received less. During ENSO cold phase, the reaches of the Changjiang (Yangtze) River and the North China received more amounts than normal rainfall in flood season of the onset years of La Niña events, and the other regions of China received less. In the following years of La Niña events, the coastal area of the Southeast China, the most part of the Northeast China and the regions between the Huanghe River and the Huaihe River received more precipitation during flood seasons, but the other parts received below normal precipitation.     

The response of an elastic half‐space under a momentary shear line impulse

The response of an ideal elastic half‐space to a line‐concentrated impulsive vector shear force applied momentarily is obtained by an analytical–numerical computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. The shear force is concentrated along an infinite line, drawn on the surface of the half‐space, while being normal to that line as well as to the axis of symmetry of the half‐space. An exact loading model is introduced and built into the computational method for this shear force. With this model, a compatibility exists among the prescribed applied force, the geometric decay of the shear stress component at the precursor shear wave, and the boundary conditions of the half‐space; in this sense, the source configuration is exact. For the transient boundary‐value problem described above, a wave characteristics formulation is presented, where its differential equations are extended to allow for strong discontinuities which occur in the material motion of the half‐space. A numerical integration of these extended differential equations is then carried out in a three‐dimensional spatiotemporal wavegrid formed by the Cartesian bicharacteristic curves of the wave characteristics formulation. This work is devoted to the construction of the computational method and to the concepts involved therein, whereas the interpretation of the resultant transient deformation of the half‐space is presented in a subsequent paper. Copyright © 2003 John Wiley & Sons, Ltd.     

The efifect of wave breaking on surface wave imaging by Synthetic Aperture Radar

The effect of ocean wave breaking as a non-Bragg mechanism on backscattering cross-section and modulation transfer functions (MTF) of radar was investigated based on Bragg resonance theory and parametric method. The result showed that the additional effect of wave breaking on backscattering cross-section is not more than 20% except for the small incident angle of VV polarized electromagnetic (e.m.) wave but is significant for HH polarized e.m. wave. Breaking waves lead to increase in the modulus of tilt modulation MTF and the larger the wind speed, the faster the increase. For large incident angle, the modulus of tilt modulation MTF with wave breaking decreases quickly with incident angle for HH polarization and approach to that without wave breaking for VV polarization. The hydrodynamic MTF increases 30%-60% when considering wave breaking and the increase is larger for HH polarization than for VV polarization.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

富氏谱分析法模拟地壳垂直运动的地学依据及其策应力模型

应用地壳波浪与镶嵌构造学说对富氏谱分析法提取地壳垂直形变信息的科学性做了地质学意义上的阐释 ,并提出了根据多期形变资料提取特定波段上构造策应力的数学模型     

Rift-related Jurassic basaltic phreatomagmatic volcanism in the central Transantarctic Mountains: precursory stage to flood-basalt effusion

 The Middle Jurassic Kirkpatrick flood basalts and comagmatic Ferrar intrusions in the Transantarctic Mountains represent a major pulse of tholeiitic magmatism related to early stages in the breakup of Gondwana. A record of the volcano-tectonic events leading to formation of this continental flood-basalt province is provided by strata underlying and only slightly predating the Kirkpatrick lavas. In the central Transantarctic Mountains, the lavas rest on widespread (≥7500 km²) tholeiitic pyroclastic deposits of the Prebble Formation. The Prebble Formation is dominated by lahar deposits and is an unusual example of a regionally developed basaltic lahar field. Related, partly fault-controlled pyroclastic intrusions cut underlying strata, and vents are represented by the preserved flanks of two small tephra cones associated with a volcanic neck. Lahar and air-fall deposits typically contain 50–60% accidental lithic fragments and sand grains derived from underlying Triassic – Lower Jurassic strata in the upper part of the Beacon Supergroup. Juvenile basaltic ash and fine lapilli consist of nonvesicular to scoriaceous tachylite, sideromelane, and palagonite, and have characteristics indicating derivation from hydrovolcanic eruptions. The abundance of accidental debris from underlying Beacon strata points to explosive phreatomagmatic interaction of basaltic magma with wet sediment and groundwater, which appears to have occurred in particular where rising magma intersected upper Beacon sand aquifers. Composite clasts in the lahar deposits exhibit complex peperitic textures formed during fine-scale intermixing of basaltic magma with wet sand and record steps in subsurface fuel-coolant interactions leading to explosive eruption. The widespread, sustained phreatomagmatic activity is inferred to have occurred in a groundwater-rich topographic basin linked to an evolving Jurassic rift zone in the Transantarctic Mountains. Coeval basaltic phreatomagmatic deposits of the Mawson and Exposure Hill Formations, which underlie exposures of the Kirkpatrick Basalt up to 1500 km to the north along strike in Victoria Land, appear to represent other parts of a regional, extension-related Middle Jurassic phreatomagmatic province which developed immediately prior to rapid outpouring of the flood basalts. This is consistent with models which assign an important role to lithospheric stretching in the generation of flood-basalt provinces. Received: 28 August 1995 / Accepted: 18 April 1996