首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 78 毫秒
1.
由于泥质所造成的附加导电现象,泥质含量及其分布形式对电阻率增大系数I和含水饱和度Sw关系具有重要影响,由于岩石物理实验中岩心孔隙结构及其组分构成、分布的微观不可调性,因而泥质分布形式所造成的影响很难通过岩心实验来单独研究。基于数字岩心的格子气自动机方法是一种有效的微观数值模拟方法,本研究利用储层岩心薄片的骨架颗粒尺寸信息资料建立数字岩心模型,结合格子气自动机技术对数字岩心不同饱和流体情况下电的传输特性进行数值模拟研究,揭示了不同泥质含量和泥质分布形式对孔隙介质导电特性非阿尔奇现象产生的影响,建立饱和度指数和泥质含量之间的关系模型,其良好的吻合性表明该方法在岩石物理研究中是一种十分有效的研究方法,而新模型适于在非阿尔奇储层进行准确的饱和度评价。  相似文献   

2.
We use the technique of direct minimization of the Gibbs free energy of the 8-component (K2O-Na2O-Fe2O3-FeO-CaO-MgO-Al2O3-SiO2) multiphase system in order to determine the equilibrium mineral assemblages of rocks of different bulk chemical compositions equilibrated at various P-T conditions. The calculated modal compositions of rocks and experimental data on elastic moduli of single crystals are then used to calculate densities and isotropic elastic wave velocities of rocks together with their pressure and temperature derivatives. Sufficient accuracy of the calculations is confirmed by comparison with experimental data on the gabbro-eclogite transformation and precise ultrasonic measurements of elastic wave velocities in a number of magmatic and metamorphic rocks.We present calculated phase diagrams with isolines of density, elastic wave velocities, and their pressure and temperature derivatives for several anhydrous magmatic rocks, from granite to lherzolite. Density and elastic properties of rocks are controlled by their chemical compositions, especially the SiO2 content, and by P-T of equilibration, and they increase with pressure due to mineral reactions changing mineral assemblages from plagioclase-bearing and garnet-free to garnetbearing and plagioclase-free. TheV p -density correlation is high, and shows two clear trends: one for iron-poor ultramafic rocks and another for all the other rocks considered. Mineral reactions, which occur at high pressures, changeV p and density of anhydrous magmatic rocks following the well-known Birch (or a similar) law.Felsic, intermediate and mafic rocks can be well distinguished in theV p -V p /V s - diagram, although their values ofV p can be close to one another. TheV p -V p /V s -density diagrams together with calculated phase diagrams can serve as efficient instruments for petrologic interpretation of seismic velocities.  相似文献   

3.
The pressure dependence of P- and S-wave velocities, velocity anisotropy, shear wave splitting and crack-porosity has been investigated in a number of samples from different crustal rock types for dry and wet (water saturated) conditions. At atmospheric pressure, P-wave velocities of the saturated, low-porosity rocks (< 1%) are significantly higher than in dry rocks, whereas the differences for S-wave velocities are less pronounced. The effect of intercrystalline fluids on seismic properties at increased pressure conditions is particularly reflected by the variation of the Poisson's ratio because P-wave velocities are more sensitive to fluids than S-wave velocities in the low-porosity rocks. Based on the experimental data, the respective crack-density parameter (), which is a measure of the number of flat cracks per volume unit contained within the background medium (crack-free matrix), has been calculated for dry and saturated conditions. There is a good correlation between the calculated crack-densities and crack-porosities derived from the experimentally determined volumetric strain curves. The shear wave velocity data, along with the shear wave polarisation referred to a orthogonal reference system, have been used to derive the spatial orientation of effective oriented cracks within a foliated biotite gneiss. The experimental data are in reasonable agreement with the self consistent model of O'Connell and Budiansky (1974). Taking the various lithologies into account, it is clear from the present study, that combined seismic measurements ofV p andV s , using theV p V s -ratio, may give evidence for fluids on grain boundaries and, in addition, may provide an estimate on the in-situ crack-densities.  相似文献   

4.
An inclusion model, based on the Kuster–Toksöz effective medium theory along with Gassmann theory, is tested to forward model velocities for fluid-saturated rocks. A simulated annealing algorithm, along with the inclusion model, effectively inverts measured compressional velocity (VP) to achieve an effective pore aspect ratio at each depth in a depth variant manner, continuously along with depth. Early Cretaceous syn-rift clastic sediments at two different depth intervals from two wells [well A (2160–2274 m) and well B (5222–5303 m)], in the Krishna–Godavari basin, India, are used for this study. Shear velocity (VS) estimated using modelled pore aspect ratio offers a high correlation coefficient (>0.95 for both the wells) with measured data. The modelled pore aspect ratio distribution suggests the decrease in pore aspect ratio for the deeper interval, mainly due to increased effective vertical stress. The pore aspect ratio analysis in relation to total porosity and volume of clay reveals that the clay volume has insignificant influence in shaping the pore geometry in the studied intervals. An approach based on multiple linear regression method effectively predicts velocity as a linear function of total porosity, the volume of clay and the modelled pore-space aspect ratio of the rock. We achieved a significant match between measured and predicted velocities. The correlation coefficients between measured and modelled velocities are considerably high (approximately 0.85 and 0.8, for VP and VS, respectively). This process indicates the possible influence of pore geometry along with total porosity and volume of clay on velocity.  相似文献   

5.
Summary In a series of triaxial experiments we have measuredV p ,V s and volumetric strain simultaneously in dilating dry and saturated rocks. For the first time these data permit quantitative comparison of seismic velocities or their ratio and dilatant volumetric strain. In air-dry samplesV p /V s decreases by a few per cent at strains of 10–3; in saturated materials with high pore pressure,V p /V s increases by a comparable amount. Decreases in seismic velocity ratio are difficult to generate in initially saturated rocks even with low pore pressures and at strain rates of 10–4/sec. A liquid-vapor transition will not produce a significant drop inV p /V s . If dilatancy and fluid flow are responsible for seismic travel time anomalies prior to earthquakes, our results suggest that such anomalies will occur only in regions where pore fluid source to sink dimensions are of the order of 10 km or more, or in regions where the rocks are not saturated to begin with.  相似文献   

6.
Converted-wave imaging in anisotropic media: theory and case studies   总被引:1,自引:0,他引:1  
Common‐conversion‐point binning associated with converted‐wave (C‐wave) processing complicates the task of parameter estimation, especially in anisotropic media. To overcome this problem, we derive new expressions for converted‐wave prestack time migration (PSTM) in anisotropic media and illustrate their applications using both 2D and 3D data examples. The converted‐wave kinematic response in inhomogeneous media with vertical transverse isotropy is separated into two parts: the response in horizontally layered vertical transverse isotrophy media and the response from a point‐scatterer. The former controls the stacking process and the latter controls the process of PSTM. The C‐wave traveltime in horizontally layered vertical transverse isotrophy media is determined by four parameters: the C‐wave stacking velocity VC2, the vertical and effective velocity ratios γ0 and γeff, and the C‐wave anisotropic parameter χeff. These four parameters are referred to as the C‐wave stacking velocity model. In contrast, the C‐wave diffraction time from a point‐scatterer is determined by five parameters: γ0, VP2, VS2, ηeff and ζeff, where ηeff and ζeff are, respectively, the P‐ and S‐wave anisotropic parameters, and VP2 and VS2 are the corresponding stacking velocities. VP2, VS2, ηeff and ζeff are referred to as the C‐wave PSTM velocity model. There is a one‐to‐one analytical link between the stacking velocity model and the PSTM velocity model. There is also a simple analytical link between the C‐wave stacking velocities VC2 and the migration velocity VCmig, which is in turn linked to VP2 and VS2. Based on the above, we have developed an interactive processing scheme to build the stacking and PSTM velocity models and to perform 2D and 3D C‐wave anisotropic PSTM. Real data applications show that the PSTM scheme substantially improves the quality of C‐wave imaging compared with the dip‐moveout scheme, and these improvements have been confirmed by drilling.  相似文献   

7.
The relationship between Vp and Vs may be used to predict Vs where only Vp is known. Vp/Vs is also used to identify pore fluids from seismic data and amplitude variation with offset analysis. Theoretical, physical, as well as statistical empirical Vp‐Vs relationships have been proposed for reservoir characterization when shear‐wave data are not available. In published work, the focus is primarily on the Vp‐Vs relationship of quartzitic sandstone. In order to broaden the picture we present Vp‐Vs relationships of greensand composed of quartz and glauconite by using data from the Paleocene greensand Nini oil field in the North Sea. A Vp‐Vs relationship derived from modelling is compared with empirical Vp‐Vs regressions from laboratory data as well as from log data. The accuracy of Vs prediction is quantified in terms of root‐mean‐square error. We find that the Vp‐Vs relationship derived from modelling works well for greensand shear‐wave velocity prediction. We model the seismic response of glauconitic greensand by using laboratory data from the Nini field. Our studies here reveal that brine‐saturated glauconitic greensand can have a similar seismic response to that from oil‐saturated quartzitic sandstone and that oil‐saturated strongly cemented greensand can have a similar amplitude variation with offset response to that from brine‐saturated weakly cemented greensand.  相似文献   

8.
A practical method is presented for determining three‐dimensional S‐wave velocity (VS) profile from microtremor measurements. Frequency–wave number (fk) spectral analyses of microtremor array records are combined, for this purpose, with microtremor horizontal‐to‐vertical (H/V) spectral ratio techniques. To demonstrate the effectiveness of the proposed method, microtremor measurements using arrays of sensors were conducted at six sites in the city of Kushiro, Japan. The spectral analyses of the array records yield dispersion characteristics of Rayleigh waves and H/V spectra of surface waves, and joint inversion of these data results in VS profiles down to bedrock at the sites. Conventional microtremor measurements were performed at 230 stations within Kushiro city, resulting in the H/V spectra within the city. Three‐dimensional VS structure is then estimated from inversion of the H/V spectra with the VS values determined from the microtremor array data. This reveals three‐dimensional VS profile of Kushiro city, together with an unknown hidden valley that crosses the central part of the city. The estimated VS profile is consistent with available velocity logs and results of subsequent borings, indicating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

9.
Three component recordings from an array of five ocean bottom seismographs in the northwestern part of the Vøring basin have been used to obtain a 2-D shear-wave (S-wave) velocity-depth model. The shear waves are identified by means of travel-time differences compared to the compressional (P) waves, and by analyzing their particle motions. The model has been obtained by kinematic (travel-time) ray-tracing modelling of the OBS horizontal components.The shear-wave modelling indicates that mode conversions occur at several high velocity interfaces (sills) in the 4–10 km depth range, previously defined by a compressional-wave velocity-depth model using the same data set.An averageV p /V s ratio of 2.1 is inferred for the layers above the uppermost sill, indicative of both poorly consolidated sediments and a low sand/shale ratio. A significant decrease in theV p /V s ratio (1.7) below the first sill may in part be atributed to well consolidated sediments, and to a change in lithology to more sandy sediments. This layer is interpreted to lie within the lower Cretaceous sequence. At 5–10 km depthV p /V s ratios of 1.85 indicate a lower sand/shale ratio consistent with the expected lithologies. The averageV p /V s ratio inferred for the crust is 1.75, which is consistent with values obtained north of Vøring, in the Lofoten area. An eastward thinning of the crystalline basement is supported by the shear-wave modelling.  相似文献   

10.
In crystalline rocks seismic velocityV p and densityp increase, whereas radioactive heat productionA decreases from acidic to basic compositions. From the velocity-density systematics for crustal rocks at different pressures an empiricalA(V p) relationship has been derived for the range 5.0–8.0 km/sec which follows the exponential law:A(V p )=a exp (-bV p ), where the numerical factorsa andb depend onin situ pressure. A graph is given by means of which the heat production distributionA(z) can be obtained for any givenV p (z) structure.Contribution No. 207, Institute of Geophysics, ETH Zurich.  相似文献   

11.
In the first part of the paper [Molodenskii, 2011], we considered the problem of ambiguity in the solution of the inverse problem of retrieval of density distribution in the Earth’s core and mantle and determination of the Q factors in the mantle from the entire set of modern data on seismic velocities (V p and V S ), the frequencies f i and quality factors Q i of free oscillations of the Earth, and the amplitudes and phases of its forced nutations. We have constructed the model distributions of these parameters, in which the root-meansquared (rms) deviations of all observed values from the predicted ones are much smaller than in the PREM model. Below, we compare the observed amplitudes of the forced nutation with the values predicted by our model. In order to understand how rigid are the constraints imposed by the amplitudes of forced nutation, we not only calculate the deviations of the observed amplitudes of nutation from the predictions by our model but also estimate the changes in these deviations caused by small variations in several parameters of the model. To the parameters to be varied we refer those which have no or barely any effect on the periods and damping constants of free oscillations but have a pronounced effect on the amplitudes of forced nutation. These parameters include (1) the rheological properties of the mantle in the interval of periods from an hour to a day; (2) the dynamical flattening of the liquid core; (3) the dynamic flattening of the solid inner core; (4) the viscosity of the liquid core; and (5) the moment of inertia of the solid inner core. In addition, we estimate the effects of variations in the moment of inertia of the liquid core to be small (±0.2%) and not to affect, within the observation error, the periods of free oscillations. We show that the uncertainty in the model depth distributions of density considerably decreases when the new data on the amplitudes and phases of the forced nutation of the Earth are taken into account. With these data, it is possible to estimate the creep function for the lower mantle in a wide range of periods from a second to a day.  相似文献   

12.
The dispersion of exhaust products of rocket fuel in the direction perpendicular to the motion of a rocket is investigated in this work. A comparison of the results of numerical calculations with a self-similar approximation of a strong cylindrically symmetric explosion is fulfilled. It is shown that at sufficiently high rocket velocity V , which exceeds the sum of gas exhaust velocity V e from the nozzle and sound speed V s (V > V e +V s ), a gasdynamic hole can arise around the rocket trajectory in the upper atmosphere, inside which the total concentration of gas becomes less than the equilibrium concentration of gas at a given altitude. The dynamics of the profiles of density and temperature of the exhaust products inside a rocket plume is calculated.  相似文献   

13.
The Khalkhab–Neshveh (KN) pluton is a part of Urumieh–Dokhtar Magmatic Arc and was intruded into a covering of basalt and andesite of Eocene to early Miocene age. It is a medium to high‐K, metaluminous and I‐type pluton ranging in composition from quartz monzogabbro, through quartz monzodiorite, granodiorite, and granite. The KN rocks show subtle differentiation trends strongly controlled by clinopyroxene, plagioclase, hornblende, apatite, and titanite, where most major elements (except K2O) are negatively correlated with SiO2; and Al2O3, Na2O, Sr, Eu, and Y define curvilinear trends. Considering three processes of magmatic differentiation including mixing and/or mingling between basaltic and dacitic magmas, gravitational fractional crystallization and in situ crystallization revealed that the latter is the most likely process for the evolution of KN magma. This is supported by the occurrence of all rock types at the same level, the lack of mafic enclaves in the granitoid rocks, the curvilinear trends of Na2O, Sr, and Eu, and the constant ratios of (87Sr/86Sr)i from quartz monzodiorite to granite (0.70475 and 0.70471, respectively). In situ crystallization took place via accumulation of plagioclase and clinopyroxene phenocrysts and concentration of these phases in the quartz monzogabbro and quartz monzodiorite at the margins of the intrusion at T ≥ 1050°C, and by filter pressing and fractionation of hornblende, plagioclase, and later biotite in the granitoids at T = ~880°C.  相似文献   

14.
Seismic modelling of the shallow subsurface (within the first few metres) is often challenging when the data are dominated by ground-roll and devoid of reflection. We showed that, even when transmission is the only available phase for analysis, fine-scale and interpretable P-wave velocity (VP) and attenuation (QP−1) models can still be prepared using full-waveform inversion, with data being preconditioned for ground-roll. To prove this idea, we suppressed the ground-roll in two different ways before full-waveform inversion modelling: first, through a bottom mute; second, through a novel wavelet transform-based method known as the redundant-lifting scheme. The applicability of full-waveform inversion is tested through imaging two buried targets. These include a pair of utility water pipes with known diameters of 0.8 m and burial depths of 1.5 m, respectively. The second target is the poorly documented backfill, which was the former location of the pipe(s). The data for full-waveform inversion are acquired along a 2D profile using a static array of 24, 40 Hz vertical component geophones and a buried point source. The results show that (a) the redundant-lifting scheme better suppresses the ground roll, which in turn provides better images of the targets in full-waveform inversion; and (b) the VP and QP−1 models from full-waveform inversion of redundant-lifting scheme data could detect the two targets adequately.  相似文献   

15.
Summary The travel-time curves of characteristic moments of stress wave patterns were investigated. The arrival time of the onset of stress wavet 0, of the maximum amplitudet 1 and of the moment terminating the pressure part (the first half-wave) of this wavet 2 were taken to be representative. These travel-time curves were used to determine the velocities of propagation of these moments (V 0,V 1 andV 2) as functions of the distance from the source. According to their variations it is possible to appoint the size of the cavity created by the explosion and to determine the distance at which the elastic source surface of stress waves is to be found. The radius of the cavity is given by 1.5 times the distance at whichV 1 separates fromV 0. The elastic source surface of stress waves is defined by the distance from whichV 0 is constant and at whichV 1 andV 2 have minimum values. These two distances determine the points at which the character of the stress wave pattern changes: from shock wave in gases, through shock wave in a solid medium to a seismic wave.Paper presented at the General Assembly of the European Seismological Commission, Luxembourg, Sept. 21–29, 1970.  相似文献   

16.
We describe an approach to the construction of an engineering geological expert system for identification of sub-bottom soils in accordance with some predefined nomenclature. The following principles of integrated interpretation of engineering geophysical and geotechnical data are presented: Firstly, the transformation of physical data (compressional- and shear-wave velocities, compressional-wave attenuation coefficients, electrical conductivity, etc.) for each of the medium points into subjective probabilities for the soil belonging to each type listed in the nomenclature, and secondly, the extrapolation of local geotechnical data (primarily drilling data) to the surrounding space by means of diffusion of the initial membership function distribution, resulting in the same set of probabilities for soil types at each point in the medium under consideration. Aggregation of the fuzzy information obtained, sufficient for reaching a conclusion for most points in the medium, is carried out by means of Bayesian summation. An example is given of integrated interpretation of real data obtained from four different sources (compressional- and shear-wave velocity sections Vp(x, z) and Vs(x, z), and two boreholes) related to the same profile.  相似文献   

17.
Velocities of compressional and shear waves in limestones   总被引:2,自引:1,他引:2  
Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (Vp and Vs, respectively) decrease with increasing porosity and that Vp decreases approximately twice as fast as Vs. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model.  相似文献   

18.
Improvements in the joint inversion of seismic and marine controlled source electromagnetic data sets will require better constrained models of the joint elastic‐electrical properties of reservoir rocks. Various effective medium models were compared to a novel laboratory data set of elastic velocity and electrical resistivity (obtained on 67 reservoir sandstone samples saturated with 35 g/l brine at a differential pressure of 8 MPa) with mixed results. Hence, we developed a new three‐phase effective medium model for sandstones with pore‐filling clay minerals based on the combined self‐consistent approximation and differential effective medium model. We found that using a critical porosity of 0.5 and an aspect ratio of 1 for all three components, the proposed model gave accurate model predictions of the observed magnitudes of P‐wave velocity and electrical resistivity and of the divergent trends of clean and clay‐rich sandstones at higher porosities. Using only a few well‐constrained input parameters, the new model offers a practical way to predict in situ porosity and clay content in brine saturated sandstones from co‐located P‐wave velocity and electrical resistivity data sets.  相似文献   

19.
Controls on sonic velocity in carbonates   总被引:2,自引:0,他引:2  
Compressional and shear-wave velocities (V p andV s) of 210 minicores of carbonates from different areas and ages were measured under variable confining and pore-fluid pressures. The lithologies of the samples range from unconsolidated carbonate mud to completely lithified limestones. The velocity measurements enable us to relate velocity variations in carbonates to factors such as mineralogy, porosity, pore types and density and to quantify the velocity effects of compaction and other diagenetic alterations.Pure carbonate rocks show, unlike siliciclastic or shaly sediments, little direct correlation between acoustic properties (V p andV s) with age or burial depth of the sediments so that velocity inversions with increasing depth are common. Rather, sonic velocity in carbonates is controlled by the combined effect of depositional lithology and several post-depositional processes, such as cementation or dissolution, which results in fabrics specific to carbonates. These diagenetic fabrics can be directly correlated to the sonic velocity of the rocks.At 8 MPa effective pressureV p ranges from 1700 to 6500 m/s, andV s ranges from 800 to 3400 m/s. This range is mainly caused by variations in the amount and type of porosity and not by variations in mineralogy. In general, the measured velocities show a positive correlation with density and an inverse correlation with porosity, but departures from the general trends of correlation can be as high as 2500 m/s. These deviations can be explained by the occurrence of different pore types that form during specific diagenetic phases. Our data set further suggests that commonly used correlations like Gardner's Law (V p-density) or the time-average-equation (V p-porosity) should be significantly modified towards higher velocities before being applied to carbonates.The velocity measurements of unconsolidated carbonate mud at different stages of experimental compaction show that the velocity increase due to compaction is lower than the observed velocity increase at decreasing porosities in natural rocks. This discrepancy shows that diagenetic changes that accompany compaction influence velocity more than solely compaction at increasing overburden pressure.The susceptibility of carbonates to diagenetic changes, that occur far more quickly than compaction, causes a special velocity distribution in carbonates and complicates velocity estimations. By assigning characteristic velocity patterns to the observed diagenetic processes, we are able to link sonic velocity to the diagenetic stage of the rock.  相似文献   

20.
For theoretically studying the intensity of the influence exerted by the polarization of the rocks on the results of direct current (DC) well logging, a solution is suggested for the direct inner problem of the DC electric logging in the polarizable model of plane-layered medium containing a heterogeneity by the example of the three-layer model of the hosting medium. Initially, the solution is presented in the form of a traditional vector volume-integral equation of the second kind (IE2) for the electric current density vector. The vector IE2 is solved by the modified iteration–dissipation method. By the transformations, the initial IE2 is reduced to the equation with the contraction integral operator for an axisymmetric model of electrical well-logging of the three-layer polarizable medium intersected by an infinitely long circular cylinder. The latter simulates the borehole with a zone of penetration where the sought vector consists of the radial Jr and Jz axial (relative to the cylinder’s axis) components. The decomposition of the obtained vector IE2 into scalar components and the discretization in the coordinates r and z lead to a heterogeneous system of linear algebraic equations with a block matrix of the coefficients representing 2x2 matrices whose elements are the triple integrals of the mixed derivatives of the second-order Green’s function with respect to the parameters r, z, r', and z'. With the use of the analytical transformations and standard integrals, the integrals over the areas of the partition cells and azimuthal coordinate are reduced to single integrals (with respect to the variable t = cos ? on the interval [?1, 1]) calculated by the Gauss method for numerical integration. For estimating the effective coefficient of polarization of the complex medium, it is suggested to use the Siegel–Komarov formula.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号