A numerical model study on multi-species harmful algal blooms coupled with background ecological fields

Based on systematized physical, chemical, and biological modules, a multi-species harmful algal bloom （HAB） model coupled with background ecological fields was established. This model schematically embod-ied that HAB causative algal species and the background ecological system, quantified as total biomass, were significantly different in terms of the chemical and biological processes during a HAB while the inter-action between the two was present. The model also included a competition and interaction mechanism between the HAB algal species or populations. The Droop equation was optimized by considering tempera-ture, salinity, and suspended material impact factors in the parameterization of algal growth rate with the nutrient threshold. Two HAB processes in the springs of 2004 and 2005 were simulated using this model. Both simulation results showed consistent trends with corresponding HAB processes observed in the East China Sea, which indicated the rationality of the model. This study made certain progress in modeling HABs, which has great application potential for HAB diagnosis, prediction, and prevention.     

Comparison of evolutionary and static modeling of stresses around a salt diapir

We compare an evolutionary with a static approach for modeling stress and deformation around a salt diapir; we show that the two approaches predict different stress histories and very different strains within adjacent wall rocks. Near the base of a rising salt diapir, significantly higher shear stresses develop when the evolutionary analysis is used. In addition, the static approach is not able to capture the decrease in the hoop stress caused by the circumferential diapir expansion, nor the increase in the horizontal stress caused by the rise of the diapir. Hence, only the evolutionary approach is able to predict a sudden decrease in the fracture gradient and identify areas of borehole instability near salt. Furthermore, the evolutionary model predicts strains an order of magnitude higher than the strains within the static model. More importantly, the evolutionary model shows significant shearing in the horizontal plane as a result of radial shortening accompanied by an almost-equivalent hoop extension. The evolutionary analysis is performed with ELFEN, and the static analysis with ABAQUS. We model the sediments using a poro-elastoplastic model. Overall, our results highlight the ability of forward evolutionary modeling to capture the stress history of mudrocks close to salt diapirs, which is essential for estimating the present strength and anisotropic characteristics of these sediments.     

Time-lapse processing of 2D seismic profiles with testing of static correction methods at the CO2 injection site Ketzin (Germany)

The Ketzin project provides an experimental pilot test site for the geological storage of CO₂. Seismic monitoring of the Ketzin site comprises 2D and 3D time-lapse experiments with baseline experiments in 2005. The first repeat 2D survey was acquired in 2009 after 22 kt of CO₂ had been injected into the Stuttgart Formation at approximately 630 m depth. Main objectives of the 2D seismic surveys were the imaging of geological structures, detection of injected CO₂, and comparison with the 3D surveys. Time-lapse processing highlighted the importance of detailed static corrections to account for travel time delays, which are attributed to different near-surface velocities during the survey periods. Compensation for these delays has been performed using both pre-stack static corrections and post-stack static corrections. The pre-stack method decomposes the travel time delays of baseline and repeat datasets in a surface consistent manner, while the latter cross-aligns baseline and repeat stacked sections along a reference horizon.Application of the static corrections improves the S/N ratio of the time-lapse sections significantly. Based on our results, it is recommended to apply a combination of both corrections when time-lapse processing faces considerable near-surface velocity changes. Processing of the datasets demonstrates that the decomposed solution of the pre-stack static corrections can be used for interpretation of changes in near-surface velocities. In particular, the long-wavelength part of the solution indicates an increase in soil moisture or a shallower groundwater table in the repeat survey.Comparison with the processing results of 2D and 3D surveys shows that both image the subsurface, but with local variations which are mainly associated to differences in the acquisition geometry and source types used. Interpretation of baseline and repeat stacks shows that no CO₂ related time-lapse signature is observable where the 2D lines allow monitoring of the reservoir. This finding is consistent with the time-lapse results of the 3D surveys, which show an increase in reflection amplitude centered around the injection well. To further investigate any potential CO₂ signature, an amplitude versus offset (AVO) analysis was performed. The time-lapse analysis of the AVO does not indicate the presence of CO₂, as expected, but shows signs of a pressure response in the repeat data.     

Use of non-linear filtering for the regional–residual separation of potential field data

Regional–residual separation is essential in gravity and magnetic data interpretation and a variety of techniques have been proposed. Graphical determination of the regional allows geological information to be taken into account. Upward continuation can be used to obtain the regional field either empirically or using some hypothesis about the geology. In some cases, a matched filter can be designed and used to separate deep and shallow sources. Simple low pass filtering has also been used but without much success. Here we propose to use a non-linear filter approach to remove gravity and magnetic anomalies smaller than a given width. This technique attempts to mimic the graphical separation method. The results from synthetic models are presented as well as the results from a case study in eastern Canada and compared to regional gravity and magnetic anomalies obtained by other techniques. Contrary to the regional fields obtained by upward continuation, non-linear filtering does not have any physical meaning. However, its main advantage is that it gives a regional component of the gravity or magnetic field similar to the one obtained from a graphical separation.     

汶川8.0级地震对鄂尔多斯块体西南缘断裂体系影响的数值模拟

利用有限元方法及静态位移解理论,采用前人反演的汶川8.0级地震破裂过程时空参数,定量计算分析汶川地震对鄂尔多斯块体西南缘断裂体系的影响.数值计算表明,具逆冲右旋性质的汶川地震,在汶川-宝鸡-线以东产生顺时针的位移扰动场,以西产生逆时针的位移扰动场,这与该分界两侧附近的活断层整体走滑方向相一致.水平剪切应力场的局部极值区...     

折射延迟时与时深曲线反演模型法静校正技术

地震波在沙丘中的传播速度与沙丘厚度密切相关.利用延迟时反演高速顶界时,采用常速或插值方法求取表层速度,不能正确反应这种变化关系,使反演模型底界存在误差,影响静校正精度.利用微测井成果数据拟和得到时深曲线,再利用延迟时从时深曲线求取表层速度和厚度可有效解决此问题.方法为:解释和分析地震微测井数据→应用地震微测井或有代表性表层结构线段拟合时深曲线→利用野外单炮初至时间获取表层延迟时→应用延迟时与时深曲线反演近地表结构模型→静校正量计算.应用结果表明,该方法快速有效,在M10JQ3D中获得较好效果.     

Calculation of CO2 column heights in depleted gas fields from known pre-production gas column heights

Depleted gas fields have been identified as potential targets for CO₂ storage. In order to maximise storage capacity, a target field must be deep enough to ensure that the CO₂ is in a dense phase (either liquid or supercritical). Accurate assessment of the storage capacity also requires an estimation of the amount of CO₂ that can be safely stored beneath the reservoir seal which can be estimated in several ways. In this paper we develop a methodology to convert known pre-production gas column heights into CO₂ column heights in order to estimate storage capacities. Several authors have correctly identified that the differences in interfacial tension and wettability act to reduce the threshold capillary entry pressure for CO₂ compared to natural gas, consequently reducing column height estimates. However, under reservoir conditions the density of CO₂ is substantially higher than natural gas so the buoyancy force on the seal for a fixed column height is much lower for dense phase CO₂ increasing column height estimates; we investigate the effects of this trade off on storage estimates and apply it to a typical dataset, in this case the UK North Sea.     

原子模型对散射偏振理论轮廓的影响

太阳磁场的诊断是太阳物理研究中非常重要的一部分.斯托克斯参量I、Q、U和V可以用来完整描述偏振光的信息,对观测得到的斯托克斯光谱进行反演可以诊断太阳磁场的信息.近几十年来,塞曼效应成为磁场诊断的主要手段,利用塞曼效应使谱线产生的分裂可以诊断强度达到几百高斯的强磁场,但是在太阳宁静区中存在大量强度小于100 Gs的弱磁场,对于弱磁场可以利用汉勒效应来诊断,应用汉勒效应诊断弱磁场一直是磁场诊断的主要内容之一,需要对偏振的产生机制有更完整的理解.文章的主要内容是研究在采用不同原子模型的假设下,由散射产生的谱线轮廓的区别,以及在存在磁场时,不同原子模型谱线的汉勒效应的区别.以中性镁为研究对象,选取了7能级、4能级和2能级原子模型,分别研究了这几个原子模型对散射偏振和汉勒效应的影响.发现2能级原子模型和多能级(7能级和4能级)原子模型产生的谱线的偏振度有很大的区别,并且在2能级原子模型假设下b_4线不存在下能级的汉勒效应,但是在多能级的原子模型下b_4线仍存在下能级的汉勒效应.对比7能级原子模型和4能级原子模型,谱线的偏振度只有很小的变化,汉勒效应的表现也基本相同.主要结论如下:在利用中性镁b 3线的线偏振轮廓Q/I反演磁场时,至少需要采用4个能级的原子模型.这主要依赖于谱线形成区域的原子的能级分布,原子能级占有数多的能级在研究中是必须考虑的.     

Numerical simulations of three-dimensional magnetic swirls in a solar flux-tube

We aim to numerically study evolution of Alfv′en waves that accompany short-lasting swirl events in a solar magnetic flux-tube that can be a simple model of a magnetic pore or a sunspot. With the use of the FLASH code we numerically solve three-dimensional ideal magnetohydrodynamic equations to simulate twists which are implemented at the top of the photosphere in magnetic field lines of the flux-tube. Our numerical results exhibit swirl events and Alfv′en waves with associated clockwise and counterclockwise rotation of magnetic lines, with the largest values of vorticity at the bottom of the chromosphere, and a certain amount of energy flux.     

Using spherical wavelets to search for magnetically-induced alignment in the arrival directions of ultra-high energy cosmic rays

Due to the action of the intervening cosmic magnetic fields, ultra-high energy cosmic rays (UHECRs) can be deflected in such a way as to create clustered energy-ordered filamentary structures in the arrival direction of these particles, the so-called multiplets. In this work we propose a new method based on the spherical wavelet transform to identify multiplets in sky maps containing arrival directions of UHECRs. The method is illustrated in simulations with a multiplet embedded in isotropic backgrounds with different numbers of events. The efficiency of the algorithm is assessed through the calculation of Type I and II errors.