Response of circular tunnel with imperfect interface bonding in layered ground subjected to obliquely incident plane P or SV wave

A semianalytical procedure is proposed for evaluating the internal forces of circular tunnel with imperfect interface bonding in layered ground subjected to an obliquely incident plane P or SV wave. In this study, the hoop bending moment and hoop axial force are related to the free‐field responses of the ground. A one‐dimensional numerical approach is firstly presented to obtain the free‐field responses of a layered half‐space with an obliquely incident plane P or SV wave propagation. Then, the free‐field stress state is transformed and decomposed in the polar coordinate system. The internal forces of tunnel caused by the isotropic stress state and the pure shear state are calculated and then summed up to obtain the overall analytical solutions. Finally, the validity of the proposed semianalytical procedure is demonstrated by numerical examples.     

Coupled and splitting bedload sediment transport models based on a modified flux-wave approach

Numerical modeling of free-surface flow over a mobile bed with predominantly bedload sediment transport can be done by solving the shallow water and Exner equations using coupled and splitting approaches.The coupled method uses a coupling of the governing equations at the same time step leading to a non-conservative solution.The splitting method solves the Exner and the shallow water equations in a separate manner,and is only capable of modeling weak free-surface and bedload interactions.In the current study,an extended version of a Godunov-type wave propagation algorithm is presented for modeling of morphodynamic systems using both coupled and splitting approaches.In the introduced coupled method the entire morphodynamic system is solved in the form of a conservation law.For the splitting technique,a new wave Riemann decomposition is defined which enables the scheme to be utilized for mild and strong interactions.To consider the bedload sediment discharge within the Exner equation,the Smart and Meyer-Peter&Müller formulae are used.It was found that the coupled solution gives accurate predictions for all investigated flow regimes including propagation over a dry-state using a Courant-Friedrichs-Lewy(CFL)number equal to 0.6.Furthermore,the splitting method was able to model all flow regimes with a lower CFL number of 0.3.     

基于全球海洋模型GOCTM的海啸传播过程数值模拟——以2004印尼大海啸为例

文中利用在有限体积近海模型FVCOM基础上拓展开发的全球海洋模型GOCTM(Global Ocean Circulation and Tide Model)进行了海啸波传播过程模拟,GOCTM采用全球无结构三角形网格,避免了开边界条件引入带来的误差,利用德国AWI研究所提供的海啸源作为初始水位场,模拟了2004年12月26日苏门答腊-安达曼Mw 9.2地震引发的海啸传播过程。通过模拟结果与印度沿岸潮位站数据以及海啸发生过程观测到的卫星高度计数据进行了对比,发现模拟结果与观测值相近,相关系数最高达0.82,相关性较好。模拟的海啸波到达苏门答腊岛北部的时间与日本的TUNAMI 模型和德国AWI研究所的TsunAWI模型的模拟结果相符,时间相差不到30 min,证明GOCTM全球模型可以较好地对海啸传播过程进行模拟,结果令人满意,希望本工作可以为我国海啸预报和预警提供参考。     

基于近邻传播聚类元胞自动机模型的武汉城市扩散和聚合过程同步模拟

元胞自动机（CA）作为城市时空动态模拟应用最广泛的模型,可以有效模拟填充式和边缘式城市扩张过程,但是在飞地式扩张模拟方面稍显不足。本文提出一种改进CA模型—APCA,在传统CA基础上利用近邻传播聚类（AP）搜寻城市扩散增长的“种子点”,实现城市增长扩散过程和聚合过程的同步模拟。以武汉市为研究区域,使用APCA模拟其在2005—2025年间城市扩张的时空过程。结果显示：① APCA在设置“种子点”数量为1~8个时模拟总体精度均高于Logistics-CA,当“种子点”数量为6时,模拟新增部分精度最高,达到0.5217;② 2015—2025年武汉市飞地型增长面积约为8.67 km²,占新增城市用地总面积比例为6.30%;③ 武汉市1995—2025年间“先扩散后聚合”的城市扩张过程符合城市增长相位理论。APCA在一定程度上了完善了传统二维平面CA框架,将城市扩张模拟维度由面维扩展到点维,为准确展现城市用地空间扩展规律提供参考。     

Finite amplitude instability thresholds in penetrative convection

Abstract

A nonlinear energy stability analysis is presented for the penetrative convection model of Veronis (1963). For top temperatures between 4°C and 8°C the nonlinear stability boundary obtained is very close to the linear one of Veronis and enables a region of possible sub-critical instabilities to be determined.     

Positioning precision analysis of GNSS multi-frequency carrier phase combinations

GPS positioning precision is affected by various error sources, and traditional combinations of GPS carrier phase observations have their own limitations such as the wide-lane, the narrow-lane and the ionospheric-free combinations. To obtain the optimal positioning precision, a new linear combination method is addressed through the variance-covariance (VCV) of the GPS multi-frequency carrier phase combination equations, and the impact of the positioning precision is analyzed with the changing of the observation errors deduced by the law of error propagation. For the high precision positioning with only one carrier phase combination, the optimal combination method is deduced and further validated by an example of a baseline resolution with 60 km length. The result indicates that this method is the simplest, and the positioning precision is the best. Therefore, it is useful for long baseline quick positioning for different precision requirements in various distances.     

基于神经网络的水力压裂监测数据校正方法研究

油气田井况地表复杂,为准确推断水力压裂裂缝方位,提出一种基于BP神经网络的压裂裂缝监测数据校正方法.本文将复杂地表条件下测点布置不均匀时的电位梯度响应,分解为正常场响应和畸变场响应.利用神经网络估算因测点布置不均匀引起的畸变场响应,通过剔除电位梯度响应中的畸变场量,获得测点布置均匀条件下的电位梯度响应,即正常场响应.采用井地电阻率法三维正演模拟裂缝推断过程,结果表明神经网络可有效估计畸变场响应,提高了利用电位梯度响应推断裂缝方位的准确度,具有良好的应用前景.     

Circulation in the South China Sea is in a state of forced oscillation: Results from a simple reduced gravity model with a closed boundary

The South China Sea (SCS) is a narrow semi-enclosed basin, ranging from 4°–6°N to 21°–22°N meridionally. It is forced by a strong annual cycle of monsoon-related wind stress. The Coriolis parameter f increases at least three times from the southern basin to the northern basin. As a result, the basin-cross time for the first baroclinic Rossby wave in the southern part of the basin is about 10-times faster than that in the northern part, which plays the most vitally important role in setting the circulation. At the northernmost edge of SCS, the first baroclinic Rossby wave takes slightly less than 1 year to move across the basin, however, it takes only 1–2 months in the southernmost part. Therefore, circulation properties for a station in the model ocean are not solely determined by the forcing at that time instance only; instead, they depend on the information over the past months. The combination of a strong annual cycle of wind forcing and large difference of basin-cross time for the first baroclinic Rossby wave leads to a strong seasonal cycle of the circulation in the SCS, hence, the circulation is dominated by the forced oscillations, rather than the quasi-steady state discussed in many textbooks.The circulation in the SCS is explored in detail by using a simple reduced gravity model forced by seasonally varying zonal wind stress. In particular, for a given time snap the western boundary current in the SCS cannot play the role of balancing mass transport across each latitude nor balancing mechanical energy and vorticity in the whole basin. In a departure from the steady wind-driven circulation discussed in many existing textbooks, the circulation in the SCS is characterized by the imbalance of mechanical energy and vorticity for the whole basin at any part of the seasonal cycle. In particular, the western boundary current in the SCS cannot balance the mass, mechanical energy, and vorticity in the seasonal cycle of the basin. Consequently, the circulation near the western boundary cannot be interpreted in terms of the wind stress and thermohaline forcing at the same time. Instead, circulation properties near the western boundary should be interpreted in terms of the contributions due to the delayed wind stress and the eastern boundary layer thickness. In fact, there is a clear annual cycle of net imbalance of mechanical energy and vorticity source/sink. Results from such a simple model may have important implications for our understanding of the complicated phenomena in the SCS, either from in-situ observations or numerical simulations.     

Modeling hydraulic fracture propagation in permeable media with an embedded strong discontinuity approach

The paper presents an embedded strong discontinuity approach to simulate single hydraulic fracture propagation in the poroelastic medium under plane-strain conditions. The method enriches the strain field with the discontinuous deformation mode and allows the fracture to be modeled inside elements. The Mode-I fracture initiation and propagation are described by the trilinear cohesive law, which is implemented by the penalty method. The enhanced permeability inside the fractured elements is dependent on the fracture aperture. Hydraulic fracture propagation is driven by the high pressure gradient near the fracture. Fluid transfer between the fracture and bulk rock is automatically captured within the poroelastic framework. The numerical framework is verified by the comparisons with the asymptotic analytical solutions for single hydraulic fracture propagation.     

Rocking vibration of a rigid circular disc in a transversely isotropic full‐space

In this paper, forced rocking vibration of a rigid circular disc placed in a transversely isotropic full‐space, where the axis of material symmetry of the full‐space is normal to the surface of the plate, is analytically investigated. Because of using the Fourier series and Hankel integral transforms, the mixed boundary‐value problem is transformed into two separate pairs of integral equations called dual integral equations. The dual integral equations involved in this paper are reduced to Fredholm integral equations of the second kind. With the aid of contour integration, the governing integral equation is numerically evaluated in the general dynamic case. The reduced static case of the dual integral equations is solved analytically and the vertical displacement, the contact pressure and the static impedance/compliance function are explicitly determined, and it is shown that the pressure in between the plate and the full‐space and the compliance function reduced for isotropic half‐space are identical to the previously published solutions. The dynamic contact pressure in between the disc and the space and also the related impedance function are numerically evaluated in general dynamic case and illustrated. It is shown that the singularity exists in the contact pressure at the edge of the disc is the same as the static case. To show the effect of material anisotropy, the numerical evaluations are given for some different transversely isotropic materials and compared. Copyright © 2010 John Wiley & Sons, Ltd.