一次典型高压型海雾过程中海上大气波导的数值模拟

2009年4月9—12日黄海海域发生了一次受高压系统影响的海雾过程。利用卫星观测与探空数据、WRF模式(Weather Research and Forecasting Model)对此次海雾过程及相伴的大气波导进行了观测分析与数值模拟。海雾与波导发展可分为3个阶段:(1)大气波导先于海雾存在于黄海海面;受高压下沉影响,黄海上空存在逆温层和较强的湿度梯度,表现为较强的贴海表面波导和非贴海表面波导。(2)海雾始于高压西部,并随高压系统逐渐东移减弱,向黄海北部扩展;辐射冷却虽然使雾顶附近逆温增强,但海雾的机械湍流使其顶部湿度梯度减小,雾顶附近对应弱悬空波导或波导消失。(3)高压系统影响使干空气下沉到雾区导致黄海海雾消散;雾顶附近逆温仍存在,同时湿度梯度增大,黄海上空逐渐变为非贴海表面波导。本研究结果表明:高压系统不仅极易为波导的发生提供有利条件,而且有利于海雾的生成,在海雾演变过程中主要是雾顶水汽梯度的变化导致了波导类型及强度的变化。     

Observations and modelling of the travel time delay and leading negative phase of the 16 September 2015 Illapel,Chile tsunami

The systematic discrepancies in both tsunami arrival time and leading negative phase (LNP) were identified for the recent transoceanic tsunami on 16 September 2015 in Illapel, Chile by examining the wave characteristics from the tsunami records at 21 Deep-ocean Assessment and Reporting of Tsunami (DART) sites and 29 coastal tide gauge stations. The results revealed systematic travel time delay of as much as 22 min (approximately 1.7％ of the total travel time) relative to the simulated long waves from the 2015 Chilean tsunami. The delay discrepancy was found to increase with travel time. It was difficult to identify the LNP from the near-shore observation system due to the strong background noise, but the initial negative phase feature became more obvious as the tsunami propagated away from the source area in the deep ocean. We determined that the LNP for the Chilean tsunami had an average duration of 33 min, which was close to the dominant period of the tsunami source. Most of the amplitude ratios to the first elevation phase were approximately 40％, with the largest equivalent to the first positive phase amplitude. We performed numerical analyses by applying the corrected long wave model, which accounted for the effects of seawater density stratification due to compressibility, self-attraction and loading (SAL) of the earth, and wave dispersion compared with observed tsunami waveforms. We attempted to accurately calculate the arrival time and LNP, and to understand how much of a role the physical mechanism played in the discrepancies for the moderate transoceanic tsunami event. The mainly focus of the study is to quantitatively evaluate the contribution of each secondary physical effect to the systematic discrepancies using the corrected shallow water model. Taking all of these effects into consideration, our results demonstrated good agreement between the observed and simulated waveforms. We can conclude that the corrected shallow water model can reduce the tsunami propagation speed and reproduce the LNP, which is observed for tsunamis that have propagated over long distances frequently. The travel time delay between the observed and corrected simulated waveforms is reduced to <8 min and the amplitude discrepancy between them was also markedly diminished. The incorporated effects amounted to approximately 78％ of the travel time delay correction, with seawater density stratification, SAL, and Boussinesq dispersion contributing approximately 39％, 21％, and 18％, respectively. The simulated results showed that the elastic loading and Boussinesq dispersion not only affected travel time but also changed the simulated waveforms for this event. In contrast, the seawater stratification only reduced the tsunami speed, whereas the earth's elasticity loading was responsible for LNP due to the depression of the seafloor surrounding additional tsunami loading at far-field stations. This study revealed that the traditional shallow water model has inherent defects in estimating tsunami arrival, and the leading negative phase of a tsunami is a typical recognizable feature of a moderately strong transoceanic tsunami. These results also support previous theory and can help to explain the observed discrepancies.     

海浪数值模式研究回顾与进展

本文简单介绍了海浪的生成、发展理论,以及海浪数值模式的分类;回顾总结了国内外在海浪数值模式研究方面取得的成果。     

Identifying Rocky Seabed Using GIS-Modeled Predictor Variables

Mapping the seabed along the Norwegian coast is costly and time consuming. Hence, finding a modeling method to separate rocky seabed from other substrate types will provide digital maps that may be used to develop cost-effective sampling designs to predict species and habitat distribution. Our approach was to use geophysical data that were quantitative and objectively defined, generalized additive models (GAMs), and Akaike information criterion (AIC) to develop statistical models and select among them. We found that slope, terrain curvature, wave exposure, and depth predicted rocky seabed occurrence with a high degree of certainty.     

Modeling an atypical petroleum system: A case study of hydrocarbon generation,migration and accumulation related to igneous intrusions in the Neuquen Basin,Argentina

In the Altiplanicie del Payún area (Neuquen Basin, Argentina), immature source rock sections intruded by up to 600 m thick Tertiary laccoliths show full spectrum maturity aureoles over hundreds of meters from the contacts. Commercial oil accumulations (20–33°API) and oil shows are located along the entire column, both in sandstone/carbonate and fractured igneous reservoirs. A challenging numerical model that included the emplacement of the intrusive bodies, with extreme temperature ranges and unusually short calculation time steps, has been done with the aim to better understand hydrocarbon generation and migration processes related to these thermal anomalies.     

A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

北部湾盆地涠西南凹陷C洼烃源岩热史及成熟史模拟

北部湾盆地涠西南凹陷C洼流沙港组烃源岩热史及成熟史研究,对C洼深水油气勘探具有指导意义。在恢复涠西南凹陷C洼地史和热史的基础上,利用EASY%Ro模型计算了流沙港组烃源岩的成熟度史。研究结果表明,在涠西南凹陷发展的裂陷阶段初始期热流值较高,最大值约为77mW/m2,其后热流值逐渐减小,现今热流值约为54mW/m2;涠西南凹陷C洼流沙港组烃源岩开始生烃(Ro=0.5%)时间为51MaBP,达到生烃高峰(Ro=1%)时间为42MaBP,达到高成熟演化阶段(Ro=1.3%)时间为17MaBP;对比涠1井流沙港组烃源岩演化特征,处于C洼深水勘探区的流沙港组烃源岩成熟度较高,生烃能力较强,拥有广阔的油气勘探前景。     

Ocean current estimation using a Multi-Model Ensemble Kalman Filter during the Grand Lagrangian Deployment experiment (GLAD)

In the summer and fall of 2012, during the GLAD experiment in the Gulf of Mexico, the Consortium for Advanced Research on Transport of Hydrocarbon in the Environment (CARTHE) used several ocean models to assist the deployment of more than 300 surface drifters. The Navy Coastal Ocean Model (NCOM) at 1 km and 3 km resolutions, the US Navy operational NCOM at 3 km resolution (AMSEAS), and two versions of the Hybrid Coordinates Ocean Model (HYCOM) set at 4 km were running daily and delivering 72-h range forecasts. They all assimilated remote sensing and local profile data but they were not assimilating the drifter’s observations. This work presents a non-intrusive methodology named Multi-Model Ensemble Kalman Filter that allows assimilating the local drifter data into such a set of models, to produce improved ocean currents forecasts. The filter is to be used when several modeling systems or ensembles are available and/or observations are not entirely handled by the operational data assimilation process. It allows using generic in situ measurements over short time windows to improve the predictability of local ocean dynamics and associated high-resolution parameters of interest for which a forward model exists (e.g. oil spill plumes). Results can be used for operational applications or to derive enhanced background fields for other data assimilation systems, thus providing an expedite method to non-intrusively assimilate local observations of variables with complex operators. Results for the GLAD experiment show the method can improve water velocity predictions along the observed drifter trajectories, hence enhancing the skills of the models to predict individual trajectories.     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.