Modeling squirt dispersion and attenuation in fluid-saturated rocks using pressure dependency of dry ultrasonic velocities

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. WA157-WA168 ◽  
Author(s):  
Osni Bastos de Paula ◽  
Marina Pervukhina ◽  
Dina Makarynska ◽  
Boris Gurevich

Modeling dispersion and attenuation of elastic waves in fluid-saturated rocks due to squirt flow requires the knowledge of a number of geometrical parameters of the pore space, in particular, the characteristic aspect ratio of the pores. These parameters are usually inferred by fitting measurements on saturated rocks to model predictions. To eliminate such fitting and thus make the model more predictive, we propose to recover the geometrical parameters of the pore space from the pressure dependency of elastic moduli on dry samples. Our analysis showed that the pressure dependency of elastic properties of rocks (and their deviation from Gassmann’s prediction) at ultrasonic frequencies is controlled by the squirt flow between equant, stiff, and so-called intermediate pores (with aspect ratios between [Formula: see text]). Such intermediate porosity is expected to close at confining pressures of between 200 and 2000 MPa, and thus cannot be directly obtained from ultrasonic experiments performed at pressures below 50 MPa. However, the presence of this intermediate porosity is inferred from the significant linear trend in the pressure dependency of elastic properties of the dry rock and the difference between the bulk modulus of the dry rock computed for spherical pores and the measured modulus at 50 MPa. Moreover, we can infer the magnitude of the intermediate porosity and its characteristic aspect ratio. Substituting these parameters into the squirt model, we have computed elastic moduli and velocities of the water-saturated rock and compared these predictions against laboratory measurements of these velocities. The agreement is good for a number of clean sandstones, but not unexpectedly worse for a broad range of shaley sandstones. Our predictions showed that dispersion and attenuation caused by the squirt flow between compliant and stiff pores may occur in the seismic frequency band. Confirmation of this prediction requires laboratory measurements of elastic properties at these frequencies.

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Prasanta Kumar Mohanta ◽  
B. T. N. Sridhar ◽  
R. K. Mishra

Abstract Experiments and simulations were carried on C-D nozzles with four different exit geometry aspect ratios to investigate the impact of supersonic decay characteristics. Rectangular and elliptical exit geometries were considered for the study with various aspect ratios. Numerical simulations and Schlieren image study were studied and found the agreeable logical physics of decay and spread characteristics. The supersonic core decay was found to be of different length for different exit geometry aspect ratio, though the throat to exit area ratio was kept constant to maintain the same exit Mach number. The impact of nozzle exit aspect ratio geometry was responsible to enhance the mixing of primary flow with ambient air, without requiring a secondary method to increase the mixing characteristics. The higher aspect ratio resulted in better mixing when compared to lower aspect ratio exit geometry, which led to reduction in supersonic core length. The behavior of core length reduction gives the identical signature for both under-expanded and over-expanded cases. The results revealed that higher aspect ratio of the exit geometry produced smaller supersonic core length. The aspect ratio of cross section in divergent section of the nozzle was maintained constant from throat to exit to reduce flow losses.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Kamran Kardel ◽  
Ali Khoshkhoo ◽  
Andres L. Carrano

Purpose The purpose of this paper is to investigate the effects of layer thickness, aspect ratio, part thickness and build orientation on distortion to have a better understanding of its behavior in material jetting technology. Design/methodology/approach Specimens with two layer thicknesses (14 and 28 µm) were printed in two aspect ratios (2:1) and (10:1), four thickness values (1, 2, 3 and 4 mm) and three build orientations (45d, XY and YX) and scanned with a wide-area 3D surface scanner to quantify distortion. The material used to build the test specimens was a commercially available resin, VeroWhitePlus RGD835. Findings The results of this study showed that all printed specimens by material jetting 3D printers had some level of distortion. The 1-mm thickness specimens, for both layer thicknesses of 14 µm and 28 µm, showed a wide range of anomalies including reverse coil set (RCS), reverse cross bow (RCB), cross bow (CB), wavy edge (WE) and some moderate twisting (T). Similar occurrences were observed for the 2-mm thickness specimens as there were RCS, WE, RCB and T anomalies that show the difference between the thinner specimens (1- and 2-mm) with the thicker ones (3- and 4-mm). In both 3- and 4-mm thickness specimens, there was more consistency in terms of distortion with mainly RCS and RCB anomalies. In total, six different types of flatness anomalies were found to occur with the following incidences: reverse coil set (91 specimens, 63.19%), reverse cross bow (50 specimens, 34.72%), wavy edge (23 specimens, 15.97%), twist (19 specimens, 12.50%), coil set (11 specimens, 7.64%) and cross bow (7 specimens, 4.86%). Originality/value This study expands the research on how the preprocess parameters such as layer thickness and build orientation and the geometrical parameters such as part thickness and aspect ratio cause dimensional distortion. Distortion is a pervasive consequence of the curing process in photopolymerization and explores one of the most common defects that come across in polymeric-based additive manufacturing. In addition to the characterization of the type and magnitude of distortion, the contributions of this work also include establishing the foundation for design guidelines aiming at minimizing distortion in material jetting.


Geophysics ◽  
2012 ◽  
Vol 77 (5) ◽  
pp. D171-D183 ◽  
Author(s):  
Bernardo Moyano ◽  
Kyle T. Spikes ◽  
Tor Arne Johansen ◽  
Nazmul Haque Mondol

Modeling the elastic properties of clay-bearing rocks (shales) requires thorough knowledge of the mineral constituents, their elastic properties, pore space microstructure, and orientations of clay platelets. Information about these variables and their complex interrelationships is rarely available for real rocks. We theoretically modeled the elastic properties of synthetic clay-water composites compacted in the laboratory, including estimates of pore space topology and percolation behavior. The mineralogy of the samples was known exactly, and the focus was on two monomineralic samples comprised of kaolinite and smectite. We used differential effective medium theory (DEM) and analysis of scanning electron microscope (SEM) images of the compacted kaolinite and smectite samples. Percolation behavior was included through calculations of critical porosities from measurements of the liquid limits of the individual clay powders. Quantitative analysis of the SEM images showed that the large scale ([Formula: see text]) pore space of the smectite composite had more rounded pores (mean aspect ratio [Formula: see text]) than the kaolinite composite (mean pore’s aspect ratio [Formula: see text]). However, models that used only these large-scale pore shapes could not explain the compressional and shear velocity measurements. DEM simulations with a single pore aspect ratio showed that bulk and shear moduli are controlled by different pore shapes. Conversely, modeling results that combined critical porosity and dual porosity models into DEM theory compared well with the measured bulk and shear moduli of compacting kaolinite and smectite composites. The methods and results we used could be used to model unconsolidated clay-bearing rocks of more complex mineralogy.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Prasanta Kumar Mohanta ◽  
B. T. N. Sridhar ◽  
R. K. Mishra

Abstract Experiments and simulations were carried on C-D nozzles with four different exit geometry aspect ratios to investigate the impact of supersonic decay characteristics. Rectangular and elliptical exit geometries were considered for the study with various aspect ratios. Numerical simulations and Schlieren image study were studied and found the agreeable logical physics of decay and spread characteristics. The supersonic core decay was found to be of different length for different exit geometry aspect ratio, though the throat to exit area ratio was kept constant to maintain the same exit Mach number. The impact of nozzle exit aspect ratio geometry was responsible to enhance the mixing of primary flow with ambient air, without requiring a secondary method to increase the mixing characteristics. The higher aspect ratio resulted in better mixing when compared to lower aspect ratio exit geometry, which led to reduction in supersonic core length. The behavior of core length reduction gives the identical signature for both under-expanded and over-expanded cases. The results revealed that higher aspect ratio of the exit geometry produced smaller supersonic core length. The aspect ratio of cross section in divergent section of the nozzle was maintained constant from throat to exit to reduce flow losses.


2021 ◽  
Author(s):  
Michele Paulatto ◽  
Joanna Morgan ◽  
Kajetan Chrapkiewicz ◽  
Emilie Hooft ◽  
Doug Toomey ◽  
...  

<p>The lack of direct seismological evidence for large molten magma chambers is considered to be one of the most important arguments in support of the mush paradigm. However, most published melt fraction estimates based on interpretation of seismological data are associated with large uncertainties because of two limitations: i) inherent limits to resolution of seismic tomography and, ii) trade-offs in the constitutive relationships that tie seismic properties to melt fraction. Low-velocity volumes associated with magma storage are particularly difficult to image with conventional travel-time tomography due to limited resolution and wavefront healing, resulting in blurred images and a high velocity bias. We tackle these limitation by applying full waveform inversion to active source seismic data collected over the Kolumbo submarine volcano (Greece). We recover a previously undetected Vp anomaly of –50% beneath the volcano and interpret this as a shallow magmatic intrusion. Extension of this approach to the wider Santorini volcanic system is ongoing. Concurrently, we are tackling the second limitation, which is the result of the dependence of elastic properties on the microgeometry of the melt. Seismological melt estimates rely on the assumption that the melt pore space can be represented by simple geometrical shapes, usually ellipsoids, with a given aspect ratio. Since the aspect ratio is poorly constrained, this results in a trade-off between melt fraction and melt geometry. We have adapted a method for the homogenisation of the elastic properties of multi-phase composites and applied it to calculating the elastic properties of partially molten rocks starting from the melt microstructure determined by X-ray CT scanning. The microgeometry of the mush can be inferred from the study of glomerocrysts: crystal mush inclusions with quenched interstitial melt that are carried to the surface by erupted lava. After the sample is digitized and segmented into its constitutive phases (crystals, melt, vesicles), the average elastic properties are determined by numerical homogenisation which consists of numerically simulating the deformation of the sample under load and predicting its elastic response. The results are compared to a semi analytical solution for ellipsoidal inclusions. We apply this approach to a plutonic nodule from St Kitts and show that the melt microstructure leads to an elastic response equivalent to that of ellipsoidal melt inclusions with an aspect ratio of 0.1 (oblate spheroids). This equivalent aspect ratio is used to refine melt estimates for Montserrat, Santorini and Kolumbo volcano.</p>


2021 ◽  
Vol 312 ◽  
pp. 08001
Author(s):  
Micol Pucci ◽  
Stefania Zanforlin ◽  
Debora Bellafiore ◽  
Stefano Deluca ◽  
Georg Umgiesser

A MATLAB routine, based on a Double Multiple Stream Tube model, developed to quickly predict the performance of cross-flow hydrokinetic turbine, here is presented. The routine evaluate flow data obtained with the open-source marine circulation code SHYFEM. The tool can establish the best locations to place tidal devices taking into account bathymetric constraints and the hydrokinetic potential. Hence, it can be used to decide the best set of geometrical parameters. The geometrical variables of our analysis are turbine frontal area, aspect ratio and solidity. Several sub-models, validated with 3D and 2D CFD simulations, reproduce phenomena such as dynamic stall, fluid dynamic tips losses and the lateral deviation of streamlines approaching the turbine. As a case study, the tool is applied to an area of the northern Adriatic Sea. After having identified some suitable sites to exploit the energy resource, we have compared behaviours of different turbines. The set of geometrical parameters that gives the best performance in terms of power coefficient can vary considering several locations. Conversely, the power production is always greater for turbine with low aspect ratio (for a fixed solidity and area). Indeed, shorter devices benefit from higher hydrokinetic potentials at the top of the water column.


Author(s):  
Nishit J. Mehta ◽  
Dilipkumar Bhanudasji Alone ◽  
Harish S. Choksi

While the effects of axisymmetric casing treatment on performance of an axial compressor stage have been extensively studied numerically as well as experimentally, the major geometrical parameters which govern these effects have been identified. Studies are now focused on understanding how each of these parameters individually impacts the performance of a casing treatment. The present work aims to study the impact on performance of casing treatment geometry when aspect ratio of the grooves is varied in a circumferential groove casing treatment. The compressor geometry chosen for this study has design characteristics of a transonic compressor stage. Flow field solutions were derived for baseline model by solving steady state 3-D Reynolds-Averaged Navier-Stokes (RANS) equations for three grid densities and the grid independence was proved. The basic casing treatment geometry has 10 circumferential grooves of width 4mm and axial spacing of 2mm between each groove. The aspect ratio was varied by changing the depth of the grooves in each case. These casing treatment geometries were superimposed over the rotor domain with the grooves extending over the entire blade tip chord and flow field solutions were again obtained for various aspect ratios of grooves. These results depict improvement in the range of operation in terms of mass flow rate. Results also show that the aspect ratio of the grooves significantly influences the overall effectiveness of casing treatment on the performance of compressor stage. Improvement in overall compressor efficiency is noted with lower aspect ratio casing treatments when compared to those with higher aspect ratios, however, the range improvement is higher with higher aspect ratios. It is also observed that, after a certain depth of grooves is reached, there is no significant improvement in performance on further increasing the depth and hence the aspect ratio. Post processing results of the flow solutions are presented which confirm the trends and show that the flow behavior near rotor tip governs this effect.


Geophysics ◽  
2012 ◽  
Vol 77 (2) ◽  
pp. R95-R104 ◽  
Author(s):  
Yang Zhang ◽  
M. Nafi Toksöz

With the current developments in imaging/computational techniques and resources, computational rock physics has been emerging as a new field of study. Properties of rocks are examined by carrying out extensive numerical simulations on rocks that have been digitized using high-resolution X-ray CT scans. The ultimate goal of computational rock physics is to supplement the traditional laboratory measurements, which are time consuming, with faster numerical simulations that allow the parameter space to be explored more thoroughly. We applied the finite-element method to compute the static effective elastic properties from 3D microtomographic images of Berea sandstone saturated with different fluids. From the computations, we found discrepancies between the numerical results and the laboratory measurements. The reason for such a problem is the loss of small features, such as fine cracks and micropores, in the digitized matrix during the imaging and aggregation process. We used a hybrid approach, combining the numerical computation and the effective media theories — the differential effective medium model and the Kuster-Toksöz model — to deduce the lost cracks by a very fast simulated annealing method. We analyzed the sensitivity of the inverted results — the distributions of crack aspect ratios and concentrations — to the clay content. We found that the inverted crack distribution is not so sensitive to clay content. Compared with the effect of cracks on the computed effective elastic properties, clay has only a secondary effect. Our approach can recover the lost cracks and is capable of predicting the effective elastic properties of the rocks from the microtomographic images for different fluid saturations. Compared with the traditional inversion schemes, based only on the effective media theories, this hybrid scheme has the advantage of utilizing the complex microstructures that are resolved in the imaging process, and it helps define the inversion space for crack distribution.


2014 ◽  
Vol 745 ◽  
pp. 180-207 ◽  
Author(s):  
Tyler Van Buren ◽  
Edward Whalen ◽  
Michael Amitay

AbstractThe formation and evolution of flow structures of a finite-span synthetic jet issuing into a quiescent flow were investigated experimentally using stereoscopic particle image velocimetry (SPIV). The effect of two geometrical parameters, the orifice aspect ratio and the neck length, were explored at a Strouhal number of 0.115 and a Reynolds number of 615. Normalized orifice neck lengths of 2, 4 and 6 and aspect ratios of 6, 12, and 18 were examined. It was found that the effect of the aspect ratio is much larger than the effect of the neck length, and as the aspect ratio increases the size of the edge vortices decreases and the presence of secondary structures is more evident. Moreover, axis switching was observed and its streamwise location increases as the aspect ratio increases. The effect of the neck length on the flow structures and the evolution of the synthetic jet was found to be secondary, where the effect was only in the very near field (i.e. close to the jet’s orifice).


2001 ◽  
Vol 80 (3-4) ◽  
pp. 315-321 ◽  
Author(s):  
J.F. Cadorin ◽  
D. Jongmans ◽  
A. Plumier ◽  
T. Camelbeeck ◽  
S. Delaby ◽  
...  

AbstractTo provide quantitative information on the ground acceleration necessary to break speleothems, laboratory measurements on samples of stalagmite have been performed to study their failure in bending. Due to their high natural frequencies, speleothems can be considered as rigid bodies to seismic strong ground motion. Using this simple hypothesis and the determined mechanical properties (a minimum value of 0.4 MPa for the tensile failure stress has been considered), modelling indicates that horizontal acceleration ranging from 0.3 m/s2 to 100 m/s2 (0.03 to 10g) are necessary to break 35 broken speleothems of the Hotton cave for which the geometrical parameters have been determined. Thus, at the present time, a strong discrepancy exists between the peak accelerations observed during earthquakes and most of the calculated values necessary to break speleothems. One of the future research efforts will be to understand the reasons of the defined behaviour. It appears fundamental to perform measurements on in situ speleothems.


Sign in / Sign up

Export Citation Format

Share Document