INVESTIGATIONS ON MOISTURE TRANSPORTS, BUDGETS AND SOURCES RESPONSIBLE FOR THE DECADAL VARIABILITY OF PRECIPITATION IN SOUTHERN CHINA

In the context of global warming, apparent decdal-interdecdal variabilities can be detected in summer precipitation in southern China. Especially around the 1990s, precipitation in South China experienced a phase transition from a deficiency regime to an abundance regime in the early 1990s, while the Yangtze River Valley witnessed a phase shift of summer precipitation from abundance to deficiency in the late 1990s. Pertinent analyses reveal a close relationship between such decadal precipitation shifts and moisture budgets, which is mainly modulated by the meridional component. This relationship can be attributed to large-scale moisture transport anomalies. Further, the HYSPLIT model is utilized to quantitatively evaluate relative moisture contributions from diverse sources during different regimes. It can be found that during the period with abundant precipitation in South China, the moisture contribution from the source of Indochina Peninsula-South China Sea increased significantly, while during the deficient precipitation regime in the Yangtze River Valley, moisture from local source, western Pacific and Indochina Peninsula-South China Sea contributed less to precipitation. It means some new features of relative moisture contributions from diverse sources to precipitation anomaly in southern China took shape after 1990s.     

GRAPES_Meso区域模式积云计算方案引进及预报效果检验

云量在模式中的准确模拟既是一个重点问题,也是一个难点问题。经评估发现,在GRAPES_Meso（Global/Regional Assimilation and PrEdiction System）模式中存在模拟云量偏少的问题,而ECMWF模式中的积云云量计算方案对云量模拟具有一定优势。鉴于此,本文借鉴ECMWF积云云量计算方案对GRAPES中尺度模式中的云量计算方案进行了适当改进,并通过影响试验模拟三维云量和地面温度来验证其改进效果。结果表明:较改进前,改进云量方案模拟的云量增多,单层云量更加准确,垂直方向云层结构改进效果明显;地面温度与观测值对比误差减小。综合分析模拟结果,认为改进云量计算方案在GRAPES模式中对云量和地面温度的模拟有较好的改进。     

输出道成像方式的共反射面元叠加方法I——理论

共反射面元（Common Reflection Surface）叠加是一种独立于宏观速度模型的零偏移距成像方法，该方法属于典型的克希霍夫型成像方法. 根据成像方式的不同，克希霍夫型成像方法可以分为两大类：输出道成像方式和输入道成像方式. 考察共反射面元叠加方法，它属于输入道成像方式. 本文基于理论模型数据，实现了输出道成缘方式的CRS叠加方法. 相比传统的输入道成像方式，它具有能够保证大偏移距反射信息的成像精度和计算效率较高的优点，而且更加容易推广到三维情形.     

Optimization and rationalization of karst groundwater resource in Feicheng Basin,China

Feicheng Basin is an important energy base in Shandong Province, China. The Ordovician limestone is the main aquifer for water supply in the basin. As a result of the current practice in which individual water users control the exploitation of the karst water, the regional groundwater level has dropped more than 30 m over the last 20 years. The groundwater is overdrawn in approximately 300 km² because of the poor management from 1981. A numerical model with distributed parameters is developed for the karst aquifer to help properly manage the water resource and sustain its utilization. The responding matrix method is used to couple the prediction model and the management model. The objective of the modeling exercises is to obtain the maximum amount of water while keeping the water level drawdown within an acceptable level. The optimum management scheme is selected among two proposed designs. The recommended measures from this study include (1) improvement of water use efficiency in industries, (2) comprehensive utilization of the water pumped from the mines, and (3) improvement of the irrigation practice. When the domestic water supply and corporates water demand take the top priorities, the numerical modeling provides a useful tool to help with the rational distribution of the karst water resources among different water users.     

某体育场边坡滑坡分析与治理方案

某校体育运动场边坡因场地平整引起多次滑坡。经勘查,分析了其中的滑坡特征及滑坡形成的机理,提出治理方案。实践证明该方案达到了预期的目的。     

基于五株采样提升算法的图像二叉分解与重构

提出了一种基于五株采样的提升算法，实现了一分为二的分解与重构。通过此算法可以构造非线性的形态小波变换，保持图像的几何信息。     

Seismic data compression based on integer wavelet transform

Due to the particularity of the seismic data, they must be treated by lossless compression algorithm in some cases. In the paper, based on the integer wavelet transform, the lossless compression algorithm is studied. Comparing with the traditional algorithm, it can better improve the compression rate. CDF (2, n) biorthogonal wavelet family can lead to better compression ratio than other CDF family, SWE and CRF, which is owe to its capability in canceling data redundancies and focusing data characteristics. CDF (2, n) family is suitable as the wavelet function of the lossless compression seismic data. Contribution No.04FE1019, Institute of Geophysics, China Earthquake Administration.     

渝怀铁路圆梁山隧道桐麻岭背斜东翼岩溶涌水突泥灾害与整治方案比选

圆梁山隧道是在建铁路重庆至怀化线的关键性控制工程,隧道全长11.068km,最大埋深约780m。隧道施工穿越桐麻岭背斜东翼时,遇到了水平循环带岩溶强烈发育地段,该地段岩溶十分发育,且连通性好;加之数次强降雨,DK361 764、DK360 873等段多次突发大规模的涌水突泥砂、块石灾害,造成部分施工机具掩埋损毁,施工掘进严重受阻,蒙受巨大的经济损失,初步分析认为涌水突泥与地下暗河或岩溶泉群的主管道网络发生较为密切的联系,可以直接接受大气降水补给。为了确保隧道施工和运营安全,针对背斜东翼段岩溶发育情况、涌水突泥特征及其降雨影响程度,进行了泄水洞排水、平导排水、封堵和堵排结合4个整治方案比选,最终采用了泄水洞方案;2003年5～9月雨季,背斜地区多次普降(特)大暴雨,DK361 764、DK360 873等地段又再次突发大规模的涌突水(泥、砂、块石),都通过兴修的泄水洞排泄掉了,隧道及其村砌结构安然无恙,表明采用泄水洞方案整治背斜东翼段的岩溶涌水突泥灾害是正确的选择。     

Response Surface Designs for Scenario Management and Uncertainty Quantification in Reservoir Production

Response surface experimental designs provide a framework for evaluating sensitivities and assessing uncertainties in reservoir-production forecasts for continuous parameters (i.e. permeability, flow rate, etc.). In this paper, the method is extended in order to integrate both continuous and discrete parameters (i.e. fault status: open/close, injection scheme: SWAG/WAG, etc.). This paper presents an appropriate experimental designs approach, notably the regression models associated with, and the statistical interpretation (sensitivity study, Monte Carlo simulations, etc.). The method has been successfully applied to a reservoir oil-production simulation problem. The objective was to define the best production scheme by optimizing the well-completion level. This application has highlighted the advantages of this new approach, both in terms of decreasing simulation cost and improving the interpretation quality.     

Runge-Kutta算法与Li差分法不同阶数配合对计算精度影响研究

为了充分发挥高阶Li空间微分方案（Li, 2005）的优点,实现了时间积分为2~6阶Runge-Kutta（简称RK）格式的偏微分方程求解算法（简称RKL算法）。然后通过多组数值试验,研究了时间积分阶数对计算误差的影响。线性平流方程的试验结果表明对于方波函数型初值,2、4、5和6阶RK算法能获得和3阶精度差不多的结果,而对于高斯函数型的初值,高阶RKL算法可以取得较好的计算效果。RK为5（6）阶时,对应的Li微分阶数可达9（10）阶,总误差控制在10-7（10-8）以内。随RK阶数增加Li微分有效阶数有增加的趋势,而总误差在逐渐减小。计算非线性无粘Burgers方程时,RKL算法能否获得好的计算结果,除了受初始场形式的影响,还与计算的目标时刻有关。当目标时刻解的各阶导数连续（且未出现无穷大数值时）,高阶（RK为4~6阶）算法是有效的;若出现了导数间断、或导数为无穷大,就会碰到冲击波解类型的问题,此时高阶RK算法也无法获得很高精度的数值解。此非线性的算例中,Li微分阶数仍然随RK阶数增加而增加,但增加的趋势不是线性的,具体变化关系可以通过实验结果拟合而获得。研究发现时间积分方案阶数大于3之后,对应的最优空间差分精度阶数可以比6阶提高很多,这再次证明了以前研究中6阶以上空间差分格式对结果无改进的现象,是由于没有使用足够高精度的时间积分方案引起的。相比于Taylor-Li（Wang,2017）算法,5~6阶的RK方法编程和实现简单,计算结果的精度比3阶算法要提高很多,因此,它是一种能够对复杂方程适用的简易高阶算法方案,具有一定的实用价值。