Estimation of Land Subsidence Based on Groundwater Flow Model

This article presents an approach for estimating land subsidence due to withdrawal of groundwater. The proposed method calculates the groundwater seepage in 3-D-condition and calculates the land subsidence one-dimensionally. The governing equation on groundwater seepage is based on the three-dimensional mass conservation law and the principle of effective stress. The land subsidence calculation method is derived based on the following assumptions: (1) displacements occur only in the vertical direction, and (2) in vertical direction the total stresses do not change. The governing equation is solved by numerical method, i.e., finite element method (FEM) in spatial discretization and finite difference method (FDM) in time series discretization. In FEM Galerkin method is adopted and in FDM, lumped matrix method is employed. The proposed method is calibrated via analyzing 1-D consolidation problem and the results are compared with those from Terzaghi's one-dimensional consolidation theory and oedometer test. The proposed method is employed to analyze the consolidation of a soft layer due to withdrawal of groundwater from an aquifer under it. Moreover, this method is also applied to a field case of land subsidence due to groundwater pumping in a gas production field in Japan. The analytical results are compared with the field observed data. The results show that this approach simulates the field case well.     

Sediment movement under unidirectional flows: an assessment of empirical threshold curves

A selection of threshold data sets from previous investigations is used to re-examine, selectively, some of the empirical curves that are most commonly used for the prediction of sediment threshold. Simple analytical formulae are derived describing single line curves for mean threshold values and envelopes for the initial movement of discrete particles and the commencement of mass sediment transport phases of the critical condition.The complexity of factors governing the prevailing conditions at threshold, together with the randomness of events, limits deterministic solutions; however, the curves provide an approach that can be used readily in practical applications. In contrast, the envelope presentation is more of a stochastic approach to the definition of threshold. It recognises that there is no definitive threshold condition under which a particular particle size can be displaced, but rather that a range of threshold values exist.The empirical threshold relationships were derived on the basis of a relatively large amount of data. However, the sources utilised have interpreted the ‘critical condition’ in various ways, implying that the formulae must be treated with caution in applications.     

Does subterranean flow initiate mud banks off the southwest coast of India?

Coastal waters off the southwest coast of India draw special attention because of the occurrence of mud banks at certain locations during southwest monsoon period. The present study puts forward a hypothesis of a subterranean flow, which could be a plausible mechanism to initiate the mud banks. The subterranean flow is believed to be coupled with activated trending faults and originate from the adjacent watershed (Vembanad Lake) separated from the sea by a narrow strip of land where submerged porous lime shell beds are present. When the lake water injection occurs through the fault, the mud/clay gets excited by its thixotrophic properties (as the overlying water looses its electrolyte) and transforms into a flowing fluid. The lowering of salinity due to the introduction of fresh water keeps the mud suspension in the water column for longer duration, leading to the formation of mud banks. The idea of subterranean flow through lime shell beds initiating formation of mud banks may apply globally to any coastal regions hugged by wetlands and of similar geological conditions.     

A Numerical Model for Natural Backfill of Pipeline Trenches Subjected to Unidirectional／Oscillatory Flows

A numerical model for the self-burial of a pipeline trench is developed. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The oscillatory flow in this study represents the action of waves. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-e turbulence closure, as well as the sediment transport equations, are solved with the finite difference method in a eurvilirrear coordinate system. Both bed and suspended loads of sediment transport are included in the morphological model. Because of the lack of experimental data on the backfilling of pipeline trenches, the numericalmodel is firstly verified against three closely-relevant experiments available in literature. A detailed measurement of the channel migration phenomenon under steady currents is employed for the assessment of the integral performance of themodel. The two experimental results from U-tube tests are used to validate the model‘s ability in predicting oscillatory flows. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.     

Estimating the geostrophic velocity obtained by HF radar observations in the upstream area of the Kuroshio

A method to extract geostrophic current in the daily mean HF radar data in the Kuroshio upstream region is established by comparison with geostrophic velocity determined from the along-track altimetry data. The estimated Ekman current in the HF velocity is 1.2% (1.5%) and 48° (38°)-clockwise rotated with respect to the daily mean wind in (outside) the Kuroshio. Furthermore, additional temporal smoothing is found necessary to remove residual ageostrophic currents such as the inertial oscillation. After removal of the ageostrophic components, the HF geostrophic velocity agrees well with that from the altimetry data with rms difference 0.14 (0.12) m/s in (outside) the Kuroshio.     

1998年春夏南海温盐结构及其变化特征

利用1998年5～8月“南海季风试验”期间“科学1”号和“实验3”号科学考察船两个航次CTD资料,分析了1998年南海夏季风暴发前后南海主要断面的温盐结构及其变化特征.观测发现,南海腹地基本被典型的南海水团所控制,但在南海东北部尤其是吕宋海峡附近,表层和次表层水明显受到西太平洋水的影响.季风暴发以后,南海北部表面温度有显著升高,升幅由西向东递减,而南海中部和南部表面温度基本没变,这使得南海北部东西向温度梯度和整个海盆南北向温度梯度均减小.北部断面表层盐度普遍由34以上降低到34以下,混合层均有所发展,是季风暴发后降水和风力加剧的结果.观测期间黑潮水跨越吕宋海峡的迹象明显但变化剧烈.4～5月,黑潮次表层水除在吕宋海峡中北部出现外,在吕宋岛以西亦有发现,表明有部分黑潮水从吕宋海峡南端沿岸向西进而向南进入南海.6～7月,次表层高盐核在吕宋海峡中北部有极大发展,但在吕宋岛以西却明显萎缩;虽然看上去黑潮水以更强的流速进、出南海,但对南海腹地动力热力结构的影响未必更大.一个超过34.55的表层高盐水体于巴拉望附近被发现,似与通过巴拉望两侧水道入侵南海的西太平洋水有关.     

电磁波信号在对流层中路径弯曲改正的级数展开式

大气折射延迟是空间大地测量中的一个主要误差源。本文严格地给出了中性大气路径弯曲改正的原理性公式。通过对各种因素量级的分析,在亚毫米量级的精度上给出了一个以视天顶距正切和正割乘积为参数的级数展开模型。新模型在球对称大气模型下,级数展开模型的精度达到了毫米量级,理论上可以展开到任意阶次,且可允许观测高度角小到5°。     

大港油田官142断块巨厚砂岩的储层流动单元

以黄骅坳陷官142断块中生界油藏为例,探讨巨厚砂岩储层流动单元的研究方法。巨厚砂岩储层流动单元研究包括两个层次,一为确定渗流屏障和连通体的分布,二为连通体内部储层渗流差异分析。研究表明,官142断块渗流屏障主要有泥质屏障、钙质砂岩胶结屏障和钙质砂砾岩胶结屏障3种类型。通过对连通体内部砂体渗流差异性分析,将连通体划分为3类流动单元。其中A类流动单元以粗孔、粗喉类型为主,渗流性能好,吸水强度大;B类流动单元储层为中孔、中喉型,渗流性能中等,吸水强度次之;C类流动单元储层多为粉细砂岩、砂砾岩或钙质胶结稍差的储层,吸水强度较差。通过流动单元的划分与研究,对预测该区的剩余油分布规律和优化调整方案提供依据。     

低渗透储层水平井—直井联合布井极限注采井距研究

低渗透储层水平井—直井联合布井技术具有广泛的应用前景,其布井方式至关重要。首先利用劈分流场的方法直观地给出常规储层水平井—直井联合布井的产能评价方法;然后针对低渗透储层的渗流特点,在考虑启动压力梯度条件下建立新的井网渗流物理模型并求解得到产能评价公式;最后经分析计算给出启动压力梯度影响的面积井网极限注采井距关系曲线。研究结果为低渗透储层部署水平井—直井联合井网提供理论依据。