混合油藏中井响应的有效算法

本文绘出了一种混合油藏生产井压力为响应的可靠算法。指出：单层解可以较容易地组合从而得到混合岩石类型的混合油藏中的压力分布。对每一层来说完井方式和外边界条件可以不同。文中通过实例证明了这一点。

三维井间电磁响应的灵敏度分析和成像方法

三维井间电磁测量能给出井间储层的电阻率分布。由于三维井间电磁响应的复杂性 ,大多数研究都基于二维或2 5维模型 ,这些简化模型只能在某些特定条件下近似实际储集体 ,在实际生产应用中受到较大限制。利用灵敏度方程的方法 ,导出了三维介质中电场和磁场对电阻率的灵敏度计算公式 ,分析了磁场垂直分量和水平分量灵敏度的空间分布。给出了三维井间电磁数据的正则化最小二乘反演及反演参数约束的实现方法。由数值模拟结果知道 ,磁场垂直分量和水平分量纵横向上表现出不同的灵敏度分布 ,在靠近源和接收器的区域具有高灵敏度 ,而在两井的中间区域为低灵敏区 ,稳定的三维成像结果显示了三维反演成像方法的正确性和有效性。

短时超大吊装动力荷载作用下盾构工作井结构响应

本文依托济南市穿黄隧道南岸盾构超深工作井,针对盾构机出井吊装过程,构建短时荷载模型,通过数值分析结合现场实测数据研究超大、短时荷载作用对超深工作井结构体系的影响。结果表明:超大短时荷载存在一定的动力特性,在荷载作用结束后主体结构变形及内力响应均出现“回弹”现象;吊装荷载对变形主要影响深度在洞门以上范围,随深度增加,墙体横向、竖向应力变化增量差值逐渐减小,表明吊装荷载影响有限;环框梁结构的存在能够明显约束主体结构变形。

Characteristics of resistivity log response of oil layers under polymer flooding

Resistivity log responses are different when the oil layers is filled by a polymer solution or water and it seriously affects the log interpretation of the water flooded layer. In this paper, we first analyze the electrical properties of the polymer solution. Then, according to the different processes for polymer flooding the reservoir, we perform rock resistivity experiments and analyze the rock resistivity variations in different displacement phases. We also compare the resistivity log responses of oil layers filled with the polymer solution to those filled with water. The results show that when displacing the oil-bearing core by different polymers, the resistivity changes monotonously decrease and show either ＂W＂ type or ＂S＂ type. The resistivity log responses are weak when displacing water and there is no flooded show if the oil layer is injected by a fresh water polymer solution. If the oil layer is injected by a sewage polymer solution, the resistivity log response has similar characteristic as a polluted water flooded layer.

川东南地区茅口组地震相分析

针对勘探程度不高、钻井分布不均且以二维地震资料为主的川东南地区,通过对茅口组地震相分析,发现了中高频席状地震相、中低频席状地震相、前积地震相、丘状地震相、波状地震相等多种类型。通过单井相分析和井旁地震响应模式的建立,结合区域地质背景分析并以沉积模式为指导开展了沉积相展布特征研究。认为茅口组沉积时期研究区在开阔台地相的整体背景下,局部发育的点礁、滩体、风化壳喀斯特均是有利储层发育部位,是隐蔽油气藏潜在勘探领域。

Analysis and application of the response characteristics of DLL and LWD resistivity in horizontal well

There exist different response characteristics in the resistivity measurements of dual laterolog （DLL） and logging while drilling （LWD） electromagnetic wave propagation logging in highly deviated and horizontal wells due to the difference in their measuring principles. In this study, we first use the integral equation method simulated the response characteristics of LWD resistivity and use the three dimensional finite element method （3D-FEM） simulated the response characteristics of DLL resistivity in horizontal wells, and then analyzed the response differences between the DLL and LWD resistivity. The comparative analysis indicated that the response differences may be caused by different factors such as differences in the angle of instrument inclination, anisotropy, formation interface, and mud intrusion. In the interface, the curves of the LWD resistivity become sharp with increases in the deviation while those of the DLL resistivity gradually become smooth. Both curves are affected by the anisotropy although the effect on DLL resistivity is lower than the LWD resistivity. These differences aid in providing a reasonable explanation in the horizontal well. However, this can also simultaneously lead to false results. At the end of the study, we explain the effects of the differences in the interpretation of the horizontal well based on the results and actual data analysis.

A study on the discrete image method for calculation of transient electromagnetic fields in geological media

We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized： the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the ＂ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.

Numerical simulation of the dual laterolog for carbonate cave reservoirs and response characteristics

Cave carbonate formations are characterized by heterogeneity, which makes electrical log prediction difficult. It is currently important to know how to use the dual laterolog to accurately identify and quantitatively evaluate caves. Using numerical simulation to calculate electrical log responses can provide a theoretical basis for cave identification and evaluation. In this paper, based on the dual laterolog principles, we first study different size spherical cave models using the finite element method （FEM）, determine a relation between resistivity and cave filling after comprehensively studying the log responses of cave models with different filling material, and finally study the dual laterolog responses on caves filled with shale, limestone, conglomerate, and thin laminated formation of sand and shale. The numerical results provide a theoretical basis for identification and evaluation of carbonate cave reservoirs.

Log identification method in complex sand reservoirs

In order to obtain effective parameters for complex sand reservoirs,a log evaluation method for relevant reservoir parameters is established based on an analysis in the gas-bearing sandstone with high porosity and low permeability,low porosity and permeability and on various characteristics of log responses to reservoir lithologies and physical properties in the Neopleozoic sand reservoir of the Ordos basin.This log evaluation method covers the Cook method that is used to evaluate the porosity and oiliness in high porosity and low permeability reservoirs and another method in which the mineral content,derived from geochemical logs,is used to identify formation lithologies.Some areas have high calcium and low silt content,not uniformly distributed,the results of which show up in the complex formation lithologies and conventional log responses with great deviation.The reliability of the method is verified by comparison with conventional log data and core analyses.The calculation results coincide with the core analytical data and gas tests,which indicate that this log evaluation method is available,provides novel ideas for study of similar complex reservoir lithologies and has some reference value.

The motion performance and mooring system of deepwater semi-submersible drilling unit

This paper analyzes the motion performance and mooring system of deepwater semi-submersible drilling unit in the district of the South China Sea using the MOSES procedure system. After the 3-D panel model of the unit was built, the 3-D diffraction-radiation theory was used to obtain the hydrodynamic loads on the wet surfaces and the response amplitude operators （RAO） of the unit. According to the environmental data, the short-term motion response to motion performance of the unit is predicted by the spectral method. Then a time-domain calculation was done to analyze the motion of the unit with its mooring system. The research results can be a reference for the model test of unit.

A case study of multi-seam coal mine entry stability analysis with strength reduction method

In this paper,the advantage of using numerical models with the strength reduction method(SRM) to evaluate entry stability in complex multiple-seam conditions is demonstrated.A coal mine under variable topography from the Central Appalachian region is used as a case study.At this mine,unexpected roof conditions were encountered during development below previously mined panels.Stress mapping and observation of ground conditions were used to quantify the success of entry support systems in three room-and-pillar panels.Numerical model analyses were initially conducted to estimate the stresses induced by the multiple-seam mining at the locations of the affected entries.The SRM was used to quantify the stability factor of the supported roof of the entries at selected locations.The SRM-calculated stability factors were compared with observations made during the site visits,and the results demonstrate that the SRM adequately identifies the unexpected roof conditions in this complex case.It is concluded that the SRM can be used to effectively evaluate the likely success of roof supports and the stability condition of entries in coal mines.

Numerical simulation of the factors influencing dust in drilling tunnels:Its application

Gas-solid two-phase flow theory was used to predict dust distribution and movement at the working face of a mine.The software package FLUENT was used to numerically simulate dust motion and the results were compared to observed data.The simulation agrees with the data taken from an actual working face,which confirms the choice of mathematical model and numerical simulation method. Using the model we predict a set of conditions optimum for reducing dust concentrations at the mine working face.

Computer simulation of evacuation in underground coal mines

In order to study the evacuation behavior of miners during accidents,we analyzed the rational layout of a safe chamber,and improved evacuation efficiency.A visual model of evacuation was developed and used to simulate the evacuation process of mines,given their special underground environments.In our simulation,the working faces of a coal mine and nearby tunnels were divided into 0.5 m×0.5 m grids to reflect the location of occupants and obstacles.Based on a "flow potential of the occupant",we determined the direction of movement and routes of occupants.In the model,evacuation speed changed as a function of crowd density,with an average speed of 1.30 m/s.The model also takes into account detection time of the disaster and the personnel response time.Evacuation time,exit flow rate and evacuation path were obtained by simulation.The results indicate that using simulation technology can present a more visual evacuation process and show the behavior of occupants.Our establishment of a mine safety evacuation system has merit as reference value.

Analysis on the Characteristics of the Co-seismic Response of Water Levels in the Jiaji Well

In this study, we analyze the co-seismic response of water levels in the Jiaji well to strong earthquakes （Ms ≥7. 8） from 2001 to 2010 at an epicentrai distance less than 8000km. We investigated the co-seismic variation form of water levels, and analyzed the relationship between the amplitude of water level variation and the magnitude and the epicentral distance. We then checked the seismic wave phases when the changes of water level occurred. It was shown that： （1） the water level＇s co-seismic response is mainly characterized by escalation with no oscillation; （2） the amplitude of water level change has a certain connection with epicentral distance and magnitude; （3） co-seismic response of water levels in the Jiaji well shows a certain directivity; （4） most of the co-seismic responses were caused by surface waves, and some by long-period S waves.

Tidal Response Characteristics of the Well Water Level in the Sichuan-Yunnan Region

In this paper,the long time series data of the well water-level data of 12 wells in the Sichuan and Yunnan area is analyzed by the Baytap-G tidal analysis software,and well water level tidal response characteristic parameters( amplitude ratio and phase change)are extracted. We analyzed the features of the shape and stage change,and characteristic parameters of the tidal response of well water level before and after the earthquakes,which can provide a new method and approach to analyzing the response relationships between well water level and earth tide and barometric pressure. The results show that Luguhu Well and 9 other wells are affected by earth tides,and their well water level amplitude ratios and phases are relatively stable; the Nanxi Well and Dayao Well water level changes are affected by the barometric pressure combined with tide force,and their well water level amplitude ratios and phases are more discrete. The water level amplitude ratios and phases of Jiangyou Well,Luguhu Well and Dongchuan Well are significant to large earthquakes,and the relationship between seismic energy density and water level amplitude ratios and phases of M_2 wave of the three wells are presented.