Structure of residual oil as a function of wettability using pore-network modelling

In the water flooding of mixed-wet porous media, oil may drain down to relatively low residual oil saturations (S_or). Various studies have indicated that such low saturations can only be reached when oil layers in pore corners are included in the pore-scale modelling. These processes within a macroscopic porous medium can be modelled at the pore-scale by incorporating the fundamental physics of capillary dominated displacement within idealised pore network models. Recently, the authors have developed thermodynamic criteria for oil layer existence in pores with non-uniform wettability which takes as input geometrically and topologically representative networks, to calculate realistic S_or values for mixed-wet and oil-wet sandstones [16, 21]. This previous work is developed in this paper to include (i) the visualisation of the 3D structure of this residual oil, and (ii) a statistical analysis of this “residual/remaining” oil. Both the visualisation and the statistical analysis are done under a wide range of wettability conditions, which is reported for the first time in this paper.The structure of residual oil for strongly water wet systems is well known (where residual = remaining oil) and our model agrees with this but this structure changes radically for mixed wet systems (where residual ≠ remaining) and this has not yet been visualised experimentally. We find that for more water-wet systems high final residual oil saturations are reached at relatively small amounts of water injected and this oil is present in the pores as bulk oil. On the other hand, for more oil-wet systems we find a slow decrease of the amount of remaining oil with increasing amounts of injected water. During the process, the remaining connectivity of the oil phase is increasingly provided by oil layers only, hence the slow drainage. The final residual oil saturation, only reached in the theoretical limit of an infinite amount of injected water, is almost entirely contained in large number of (relatively low volume) oil layers, which are present in pores of most radius sizes.     

Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals

Possible configurations of the magnetic field in the outer magnetosphere during geomagnetic polarity reversals are investigated by considering the idealized problem of a magnetic multipole of order m and degree n located at the centre of a spherical cavity surrounded by a boundless perfect diamagnetic medium. In this illustrative idealization, the fixed spherical (magnetopause) boundary layer behaves as a perfectly conducting surface that shields the external diamagnetic medium from the compressed multipole magnetic field, which is therefore confined within the spherical cavity. For a general magnetic multipole of degree n, the non-radial components of magnetic induction just inside the magnetopause are increased by the factor 1 + [(n + 1)/n] relative to their corresponding values in the absence of the perfectly conducting spherical magnetopause. An exact equation is derived for the magnetic field lines of an individual zonal (m = 0), or axisymmetric, magnetic multipole of arbitrary degree n located at the centre of the magnetospheric cavity. For such a zonal magnetic multipole, there are always two neutral points and n – 1 neutral rings on the spherical magnetopause surface. The two neutral points are located at the poles of the spherical magnetopause. If n is even, one of the neutral rings is coincident with the equator; otherwise, the neutral rings are located symmetrically with respect to the equator. The actual existence of idealized higher-degree (n > 1) axisymmetric magnetospheres would necessarily imply multiple (n + 1) magnetospheric cusps and multiple (n) ring currents. Exact equations are also derived for the magnetic field lines of an individual non-axisymmetric magnetic multipole, confined by a perfectly conducting spherical magnetopause, in two special cases; namely, a symmetric sectorial multipole (m = n) and an antisymmetric sectorial multipole (m = n – 1). For both these non-axisymmetric magnetic multipoles, there exists on the spherical magnetopause surface a set of neutral points linked by a network of magnetic field lines. Novel magnetospheric processes are likely to arise from the existence of magnetic neutral lines that extend from the magnetopause to the surface of the Earth. Finally, magnetic field lines that are confined to, or perpendicular to, either special meridional planes or the equatorial plane, when the multipole is in free space, continue to be confined to, or perpendicular to, these same planes when the perfectly conducting magnetopause is present.Also Honorary Research Associate, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK and Visiting Reader in Physics. University of Sussex, Falmer, Brighton BN1 9QH, UK     

Analysis of thick, non-planar boundaries using the discontinuity analyser

The advent of missions comprised of phased arrays of spacecraft, with separation distances ranging down to at least mesoscales, provides the scientific community with an opportunity to accurately analyse the spatial and temporal dependencies of structures in space plasmas. Exploitation of the multi-point data sets, giving vastly more information than in previous missions, thereby allows unique study of their small-scale physics. It remains an outstanding problem, however, to understand in what way comparative information across spacecraft is best built into any analysis of the combined data. Different investigations appear to demand different methods of data co-ordination. Of the various multi-spacecraft data analysis techniques developed to affect this exploitation, the discontinuity analyser has been designed to investigate the macroscopic properties (topology and motion) of boundaries, revealed by multi-spacecraft magnetometer data, where the possibility of at least mesoscale structure is considered. It has been found that the analysis of planar structures is more straightforward than the analysis of non-planar boundaries, where the effects of topology and motion become interwoven in the data, and we argue here that it becomes necessary to customise the analysis for non-planar events to the type of structure at hand. One issue central to the discontinuity analyser, for instance, is the calculation of normal vectors to the structure. In the case of planar and ‘thin’ non-planar structures, the method of normal determination is well-defined, although subject to uncertainties arising from unwanted signatures. In the case of ‘thick’, non-planar structures, however, the method of determination becomes particularly sensitive to the type of physical sampling that is present. It is the purpose of this article to firstly review the discontinuity analyser technique and secondly, to discuss the analysis of the normals to thick non-planar structures detected in magnetometer data.     

High-Altitude CUSP: The Extremely Dynamic Region in Geospace

The high-altitude dayside cusps (both northern and southern) are extremely dynamic regions in geospace. Large diamagnetic cavities with significant fluctuations of the local magnetic field strength have been observed there. These cusp diamagnetic cavities are always there day after day and are as large as 6 R_E Associated with these cavities are charged particles with energies from 20 keV up to 10 MeV. The intensities of the cusp energetic ions have been observed to increase by as much as four orders of the magnitude when compared with regions adjacent to the cusp which includes the magnetosheath. Their seed populations are a mixture of ionospheric and solar wind particles. The measured energetic ion fluxes in the high-altitude cusp are higher than that in both the regions upstream and downstream from the bow shock. Turbulent electric fields with an amplitude of about 10 mV/m are also present in the cusp, and a cusp resonant acceleration mechanism is suggested. The observations indicate that the dayside high-altitude cusp is a key region for transferring the solar wind mass, momentum, and energy into the Earth’s magnetosphere.     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

Are our ideas about Dst correct?

The idea of two separate storm time ring currents, a symmetric and an asymmetric one has accepted since the 1960s. The existence of a symmetric equatorial ring current was concluded from Dst. However, the asymmetric development of the low-latitude geomagnetic disturbance field during storms lead to the assumption of the real existence of an asymmetric ring current. I think it is time to inquire whether this conception is correct. Thus, I have investigated the development of the low-latitude geomagnetic field during all the magnetic local times under disturbed and quiet conditions. The storm on February 6–9, 1986 and a statistical analysis of many storms has shown that the asymmetry does not vanish during the storm recovery phase. The ratio between the recovery phase asymmetry and the main phase asymmetry is low only for powerful storms. Storms of moderate intensity show the opposite. The global picture of the field evolution of the February storm shows clear differences at different local times. For instance the main phase and recovery phase start time does not coincide with Dst. Also the ring current decay is not the same at different local times. Therefore, Dst gives an incorrect picture of the field development. Moreover, asymmetry does not disappear during international quiet days as the investigation of the low-latitude geomagnetic field shows. Considering all these observations, I think we must revise our ideas about the ring current. In my opinion only one ring current exists and this is an asymmetric one. This asymmetry increases during storms and develops rather fast to more or less symmetric conditions. However, in no case is itjustified to conclude from Dst that a symmetric ring current exists.     

Spatial distribution of upstream magnetospheric ≥50 keV ions

We present for the first time a statistical study of 50 keV ion events of a magnetospheric origin upstream from Earths bow shock. The statistical analysis of the 50–220 keV ion events observed by the IMP-8 spacecraft shows: (1) a dawn-dusk asymmetry in ion distributions, with most events and lower intensities upstream from the quasi-parallel pre-dawn side (4 LT-6 LT) of the bow shock, (2) highest ion fluxes upstream from the nose/dusk side of the bow shock under an almost radial interplanetary magnetic field (IMF) configuration, and (3) a positive correlation of the ion intensities with the solar wind speed and the index of geomagnetic index Kp, with an average solar wind speed as high as 620 km s^–1 and values of the index Kp 2. The statistical results are consistent with (1) preferential leakage of 50 keV magnetospheric ions from the dusk magnetopause, (2) nearly scatter free motion of 50 keV ions within the magnetosheath, and (3) final escape of magnetospheric ions from the quasi-parallel dawn side of the bow shock. An additional statistical analysis of higher energy (290–500 keV) upstream ion events also shows a dawn-dusk asymmetry in the occurrence frequency of these events, with the occurrence frequency ranging between 16%-34% in the upstream region.     

Pc3 pulsations during variable IMF conditions

Pc3 geomagnetic field fluctuations detected at low latitude (L’Aquila, Italy) during the passage of a high velocity solar wind stream, characterized by variable interplanetary magnetic field conditions, are analyzed. Higher frequency resonant fluctuations and lower frequency phenomena are simultaneously observed; the intermittent appearance and the variable frequency of the longer period modes can be well interpreted in terms of the variable IMF elements; moreover their polarization characteristics are consistent with an origin related to external waves propagating in antisunward direction. A comparison with simultaneous observations performed at Terra Nova Bay (Antarctica) provides additional evidence for a clear relationship between the IMF and Pc3 pulsations also at very high latitudes.     

Signatures of storm sudden commencements in geomagnetic H, Y and Z fields at Indian observatories during 1958–1992

The work describes an intensive study of storm sudden commencement (SSC) impulses in horizontal (H), eastward (Y) and vertical (Z) fields at four Indian geomagnetic observatories between 1958–1992. The midday maximum of δH has been shown to exist even at the low-latitude station Alibag which is outside the equatorial electrojet belt, suggesting that SSC is associated with an eastward electric field at equatorial and low latitudes. The impulses in Y field are shown to be linearly and inversely related to δH at Annamalainagar and Alibag. The average SC disturbance vector is shown to be about 10–20°W of the geomagnetic meridian. The local time variation of the angle is more westerly during dusk hours in summer and around dawn in the winter months. This clearly suggests an effect of the orientation of shock front plane of the solar plasma with respect to the geomagnetic meridian. The δZ at δC have a positive impulse as in δH. The ratio of δZ/δH are abnormally large exceeding 1.0 in most of the cases at Trivandrum. The latitudinal variation of δZ shows a tendency towards a minimum over the equator during the nighttime hours. These effects are explained as (1) resulting from the electromagnetic induction effects due to the equatorial electrojet current in the subsurface conducting layers between India and Sri Lanka, due to channelling of ocean currents through the Palk Strait and (2) due to the concentration of induced currents over extended latitude zones towards the conducting graben between India and Sri Lanka just south of Trivandrum.