Search NASA⌕ Search

SEARCH · Search NASA

Results for “matched field processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

Predicting Decade-to-Century Climate Change: Prospects for Improving Models

Recent research has led to a greatly increased understanding of the uncertainties in today's climate models. In attempting to predict the climate of the 21st century, we must confront not only computer limitations on the affordable resolution of global models, but also a lack of physical realism in attempting to model key processes. Until we are able to incorporate adequate treatments of critical elements of the entire biogeophysical climate system, our models will remain subject to these uncertainties, and our scenarios of future climate change, both anthropogenic and natural, will not fully meet the requirements of either policymakers or the public. The areas of most-needed model improvements are thought to include air-sea exchanges, land surface processes, ice and snow physics, hydrologic cycle elements, and especially the role of aerosols and cloud-radiation interactions. Of these areas, cloud-radiation interactions are known to be responsible for much of the inter-model differences in sensitivity to greenhouse gases. Recently, we have diagnostically evaluated several current and proposed model cloud-radiation treatments against extensive field observations. Satellite remote sensing provides an indispensable component of the observational resources. Cloud-radiation parameterizations display a strong sensitivity to vertical resolution, and we find that vertical resolutions typically used in global models are far from convergence. We also find that newly developed advanced parameterization schemes with explicit cloud water budgets and interactive cloud radiative properties are potentially capable of matching observational data closely. However, it is difficult to evaluate the realism of model-produced fields of cloud extinction, cloud emittance, cloud liquid water content and effective cloud droplet radius until high-quality measurements of these quantities become more widely available. Thus, further progress will require a combination of theoretical and modeling research, together with intensified emphasis on both in situ and space-based remote sensing observations.

Somerville, Richard C. J.↗

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of free-stream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Re(ex)=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

LES-LDKM↗

New Developments in MadGraph/MadEvent

We here present some recent developments of MadGraph/MadEvent since the latest published version, 4.0. These developments include: Jet matching with Pythia parton showers for both Standard Model and Beyond the Standard Model processes, decay chain functionality, decay width calculation and decay simulation, process generation for the Grid, a package for calculation of quarkonium amplitudes, calculation of Matrix Element weights for experimental events, automatic dipole subtraction for next-to-leading order calculations, and an interface to FeynRules, a package for automatic calculation of Feynman rules and model files from the Lagrangian of any New Physics model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Non-perturbative effects in $μ → eγ$

We compute the non-perturbative contribution of semileptonic tensor operators $(\bar{q}σ^{μν}q)(\bar{ℓ}σ_{μv}ℓ)$ to the purely leptonic process μ → eγ and to the electric and magnetic dipole moments of charged leptons by matching onto chiral perturbation theory at low energies. This matching procedure has been used extensively to study semileptonic and leptonic weak decays of hadrons. In this paper, we apply it to observables that contain no strongly interacting external particles. The non-perturbative contribution to μ → e processes is used to extract the best current bound on lepton-flavor-violating semileptonic tensor operators, Λ BSM ≳ 450 TeV. We briefly discuss how the same method applies to dark-matter interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

From underlying event sensitive to insensitive: factorization and resummation

In this study we study the transverse energy spectrum for the Drell-Yan process. The transverse energy is measured within the central region defined by a (pseudo-) rapidity cutoff. Soft-collinear effective theory (SCET) is used to factorize the cross section and resum large logarithms of the rapidity cutoff and ratios of widely separated scales that appear in the fixed order result. We develop a framework which can smoothly interpolate between various regions of the spectrum and eventually match onto the fixed order result. This way a reliable calculation is obtained for the contribution of the initial state radiation to the measurement. By comparing our result for Drell-Yan against Pythia we obtain a simple model that describes the contribution from multiparton interactions (MPI). A model with little or no dependence on the primary process gives results in agreement with the simulation. Based on this observation we propose MPI insensitive measurements.

74 ATOMIC AND MOLECULAR PHYSICS↗

Matching next-to-leading order predictions to parton showers in supersymmetric QCD

We present a fully automated framework based on the FeynRules and MadGraph5_aMC@NLO programs that allows for accurate simulations of supersymmetric QCD processes at the LHC. Starting directly from a model Lagrangian that features squark and gluino interactions, event generation is achieved at the next-to-leading order in QCD, matching short-distance events to parton showers and including the subsequent decay of the produced supersymmetric particles. As an application, we study the impact of higher-order corrections in gluino pair-production in a simplified benchmark scenario inspired by current gluino LHC searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measuring Low-Order Aberrations in a Segmented Telescope

The in-focus PSF optimizer (IPO) is an algorithm for use in monitoring and controlling the alignment of the segments of a segmented-mirror astronomical telescope. IPO is so named because it computes wave-front aberrations of the telescope from digitized pointspread functions (PSFs) measured in infocus images. Inasmuch as distant astronomical objects that behave optically as point sources can typically be seen in almost any astronomical image, the main benefit afforded by IPO may be to enable maintenance of mirror-segment alignments without detracting from valuable scientific-observation time. IPO evolved from prescription-retrieval type algorithms. Prescription retrieval uses in-focus and out-of-focus PSFs to infer the state of an imaging optical system. The state, in this context, refers to the positions, orientations, and low-order figure errors of the optical elements in the system. Both prescription- retrieval and IPO use an iterative, nonlinear, least-squares optimizer to compute the optimal state parameters such that a digital computer-generated model image matches the digitized image acquired from the real system. The difference between IPO and prescription- retrieval algorithms is that IPO is specifically designed to utilize infocus images only. Although the restriction to in-focus images limits IPO to calculating only the lowest-order wave front aberrations, it also causes the resulting computation to take much less time because fewer degrees of freedom are included in the optimization process. In the prescription retrieval software developed at JPL, the model images are generated using the ray-trace/physical optics program, MACOS. IPO, on the other hand, uses a linear sensitivity matrix to compute the exit-pupil wave front from the system parameters; the wave front is then converted into a complex pupil field, which is then propagated to the image plane via a fast Fourier transform. This approach is computationally faster and requires less computer memory than is needed for prescription retrieval.

Ohara, Catherine↗

Precise predictions for same-sign W-boson scattering at the LHC

Vector-boson scattering processes are of great importance for the current run-II and future runs of the Large Hadron Collider. The presence of triple and quartic gauge couplings in the process gives access to the gauge sector of the Standard Model (SM) and possible new-physics contributions there. To test any new-physics hypothesis, sound knowledge of the SM contributions is necessary, with a precision which at least matches the experimental uncertainties of existing and forthcoming measurements. In this article we present a detailed study of the vector-boson scattering process with two positively-charged leptons and missing transverse momentum in the final state. In particular, we first carry out a systematic comparison of the various approximations that are usually performed for this kind of process against the complete calculation, at LO and NLO QCD accuracy. Such a study is performed both in the usual fiducial region used by experimental collaborations and in a more inclusive phase space, where the differences among the various approximations lead to more sizeable effects. Afterwards, we turn to predictions matched to parton showers, at LO and NLO: we show that on the one hand, the inclusion of NLO QCD corrections leads to more stable predictions, but on the other hand the details of the matching and of the parton-shower programs cause differences which are considerably larger than those observed at fixed order, even in the experimental fiducial region. We conclude with recommendations for experimental studies of vector-boson scattering processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Findings on Subtask 3.1 - Bakken Rich Gas Enhanced Oil Recovery Project

Total in-place oil for the Bakken petroleum system (BPS) (which includes the Bakken and Three Forks Formations) has been estimated to be 600 billion barrels (bbl). However, BPS wells have decline rates as high as 85% over the first 3 years of their lives, and primary recovery factors typically range from 3% to 10% of original oil in place. Given the low initial recovery rates, even small incremental productivity improvements could dramatically increase technically recoverable oil in the BPS. One potential solution is enhanced oil recovery (EOR) using gas injection, such as carbon dioxide (CO2) or hydrocarbon (HC) gases. While commonly used in conventional reservoirs, CO2 EOR in unconventional tight oil reservoirs has been limited to pilot tests. EOR using rich gas (mixture of methane, ethane, and propane) has also been employed in numerous pilots in several unconventional plays and has recently been successfully applied in the Eagle Ford play. If successful, large-scale gas-based EOR in the BPS could dramatically increase oil productivity and recovery factors and extend the life of the play for decades. While CO2 may be a technically suitable working fluid for EOR in the BPS, supplies are limited and costs for using CO2 in EOR pilots are prohibitively high. Meanwhile, produced gas flaring has presented challenges for BPS operators in North Dakota. Analysis conducted by the North Dakota Pipeline Authority indicates that the current gas-gathering infrastructure in North Dakota is insufficient to accommodate all of the associated gas that is produced from the BPS. The geographically isolated location of North Dakota relative to large natural gas markets, combined with suppressed natural gas prices, has made it economically challenging for industry to invest capital in expanding gas-gathering infrastructure in the state. These circumstances led to a research program conducted by the Energy & Environmental Research Center (EERC) in partnership with Liberty Resources Management Company LLC (LR) to examine the potential to use rich gas injection for EOR and mitigate flaring. A rich gas EOR pilot test was designed and executed by LR at its Stomping Horse development area in Williams County, North Dakota. From July 2018 through May 2019, a total of 160 million standard cubic feet (MMscf) of rich produced gas was injected into the BPS using five different wells in a sequential injection strategy. LR’s Leon–Gohrick drill spacing unit (DSU) was used as the test site. Regulatory oversight was provided by the North Dakota Industrial Commission (NDIC). Technical support was provided by the EERC through a series of laboratory, modeling, and field-based activities, and additional post-pilot research activities incorporated learnings from the test, developed new laboratory data, improved fracture modeling methods, and developed machine learning and big data analytics. The results from the Stomping Horse rich gas EOR pilot activities indicate that developing an effective, economical EOR approach for the BPS will require more field tests. Another key lesson learned from the Stomping Horse tests is that detailed pre- and posttest data on reservoir conditions and fluids production are essential. Robust reservoir characterization provides information that is crucial to creating realistic geomodels and conducting valid dynamic simulations of potential EOR scenarios. A detailed understanding of the completions and production history of offset wells is also necessary for valid test result interpretations. This knowledge is essential to designing the operational parameters of injectivity tests and interpreting the results. A conformance control strategy is also essential to success. Laboratory-based examinations of rich gas interactions with reservoir fluids and rocks were conducted, with an emphasis on determining the ability to mobilize oil in the tight reservoir rocks and shales of the BPS. Injection fluid composition was shown to have a positive impact on reducing reservoir oil minimum miscibility pressure (MMP), reducing interfacial tension (IFT), and altering wettability. IFT and contact angle measurements demonstrated that wettability can be altered in the presence of rich gas, suggesting the potential to improve oil recovery. Iterative modeling of surface infrastructure and reservoir performance using data generated by the various project activities was conducted. A geologic model of the Stomping Horse area was built; history-matched oil, gas, and water production was used in simulations of various EOR scenarios. Early programmatic modeling results were used to support LR’s design and operation of the EOR pilot and to provide insight regarding optimization of future commercial-scale BPS EOR design and operations. Post-pilot modeling focused on alternative methods of understanding complex fracture networks and accelerating simulation time. These led to improved simulation run times and provide excellent history-matching results. Several of these iterative models were used as the bases for developing algorithms into machine learning and big data analytics. History matching in reservoir simulation is time-consuming and computer processing-intensive. Machine learning algorithms were created, and an automated history-matching tool was developed. A large set of synthetic reservoir simulations were created to generate well responses (oil, gas, and water production, well bottomhole pressure [BHP], and tracer or propane breakthrough) for a set of EOR operating parameters that included offset well status (open or closed), injectate (rich gas or propane), injection rate, and injection well BHP. A user interface was developed to provide real-time visualization. Machine learning-based models were developed to provide rapid forecasting of well performance given a set of user-defined EOR operating parameters. These predictive models allow the user to modify the offset well status, injection rate, and injection well BHP and rapidly forecast future production performance. The combination of real-time visualization tools with real-time forecasting tools provides a framework for real-time control—operational changes that the EOR site operator can enact (e.g., changing gas injection rates) to affect the observed performance and potentially improve the EOR outcome. There is great reason to be optimistic about the future of EOR in the Bakken. The results of the laboratory studies suggest significant potential for high rates of oil mobilization using produced field gas injection under the right conditions. The results of the lab studies, combined with rigorous statistical analysis of well production data and associated modeling efforts, confirm the notion that fluid mobility within the reservoir is controlled by fractures. As more knowledge is gained about the nature and distribution of fracture networks in the Bakken, the industry will be in a better position to predict and, ultimately, influence fluid mobility. New field tests are necessary to develop a more complete understanding of those conditions. Thoughtful and creatively engineered field tests within a well-characterized geologic setting will yield the fundamental knowledge needed to take Bakken oil production to the next level. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission’s Oil and Gas Research Program and Computer Modelling Group.

04 OIL SHALES AND TAR SANDS↗

Lightweight jet reconstruction and identification as an object detection task

We apply object detection techniques based on deep convolutional blocks to end-to-end jet identification and reconstruction tasks encountered at the CERN large hadron collider (LHC). Collision events produced at the LHC and represented as an image composed of calorimeter and tracker cells are given as an input to a Single Shot Detection network. The algorithm, named PFJet-SSD performs simultaneous localization, classification and regression tasks to cluster jets and reconstruct their features. This all-in-one single feed-forward pass gives advantages in terms of execution time and an improved accuracy w.r.t. traditional rule-based methods. A further gain is obtained from network slimming, homogeneous quantization, and optimized runtime for meeting memory and latency constraints of a typical real-time processing environment. We experiment with 8-bit and ternary quantization, benchmarking their accuracy and inference latency against a single-precision floating-point. We show that the ternary network closely matches the performance of its full-precision equivalent and outperforms the state-of-the-art rule-based algorithm. Finally, we report the inference latency on different hardware platforms and discuss future applications.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Thick-film materials for silicon photovoltaic cell manufacture

Thick film technology is applicable to three areas of silicon solar cell fabrication; metallization, junction formation, and coating for protection of screened ohmic contacts, particularly wrap around contacts, interconnection and environmental protection. Both material and process parameters were investigated. Printed ohmic contacts on n- and p-type silicon are very sensitive to the processing parameters of firing time, temperature, and atmosphere. Wrap around contacts are easily achieved by first printing and firing a dielectric over the edge and subsequently applying a low firing temperature conductor. Interconnection of cells into arrays can be achieved by printing and cofiring thick film metal pastes, soldering, or with heat curing conductive epoxies on low cost substrates. Printed (thick) film vitreous protection coatings do not yet offer sufficient optical uniformity and transparency for use on silicon. A sprayed, heat curable SiO2 based resin shows promise of providing both optical matching and environmental protection.

Field, M. B.↗

Toward a Unified Parameterization of the Boundary Layer and Moist Convection: Simulations of Clear and Cloudy Convection - Part 3

A model that employs a new form of mass-flux closure (described in Part I of this paper) is applied to a variety of clear and cloudy planetary boundary layers (PBLs) including dry convection from the Wangara Experiment, trade wind cumulus from the Barbados Oceanographic and Meteorological Experiment (BOMEX), and marine stratocumulus from the Atlantic Stratocumulus Experiment (ASTEX). For Wangara, the simulated variances and fluxes match that expected from similarity arguments, while the mean state is a little less mixed than the observations. In the BOMEX simulation. the shape and magnitude of the fluxes and the turbulence kinetic energy budget agree with LES results and observations. However, the liquid water mixing ratio is too large. This is attributed to an underprediction of the skewness. In agreement with observations from the ASTEX experiment, many of the model-simulated fields distinctly reflect a regime in transition between the trade wind cumulus and the classic stratocumulus-topped boundary layers. In general, the simulated entrainment rate tends to be a little underpredicted in regimes where there is little cloud-top radiative cooling (Wangara and BOMEX), while it is overpredicted in regimes where this process is more critical (e.g., ASTEX). Prior work suggests that this may be related to the manner in which the pressure terms are parameterized in the model. Overall, the model is able to capture some key physical features of these PBL regimes, and appears to have the potential to represent both cloud and boundary layer processes. Thus, this approach is a first step toward unifying these processes in large-scale models.

Lappen, Cara-Lyn↗

Comparisons of Observations with Results from 3D Simulations and Implications for Predictions of Ozone Recovery

Although chemistry and transport models (CTMs) include the same basic elements (photo- chemical mechanism and solver, photolysis scheme, meteorological fields, numerical transport scheme), they produce different results for the future recovery of stratospheric ozone as chlorofluorcarbons decrease. Three simulations will be contrasted: the Global Modeling Initiative (GMI) CTM driven by a single year\'s winds from a general circulation model; the GMI CTM driven by a single year\'s winds from a data assimilation system; the NASA GSFC CTM driven by a winds from a multi-year GCM simulation. CTM results for ozone and other constituents will be compared with each other and with observations from ground-based and satellite platforms to address the following: Does the simulated ozone tendency and its latitude, altitude and seasonal dependence match that derived from observations? Does the balance from analysis of observations? Does the balance among photochemical processes match that expected from observations? Can the differences in prediction for ozone recovery be anticipated from these comparisons?

Douglass, Anne R.↗

Ultra-High Temperature Materials Characterization for Propulsion Applications

Propulsion system efficiency increases as operating temperatures are increased. Some very high-temperature materials are being developed, including refractory metal alloys, carbides, borides, and silicides. System design requires data for materials properties at operating temperatures. Materials property data are not available for many materials of interest at the desired operating temperatures (up to approx. 3000 K). The objective of this work is to provide important physical property data at ultra-high temperatures. The MSFC Electrostatic levitation (ESL) facility can provide measurements of thermophysical properties which include: creep strength, density and thermal expansion for materials being developed for propulsion applications. The ESL facility uses electrostatic fields to position samples between electrodes during processing and characterization studies. Because the samples float between the electrodes during studies, they are free from any contact with a container or test apparatus. This provides a high purity environment for the study of high-temperature, reactive materials. ESL can be used to process a wide variety of materials including metals, alloys, ceramics, glasses and semiconductors. The MSFC ESL has provided non-contact measurements of properties of materials up to 3400 C. Density and thermal expansion are measured by analyzing digital images of the sample at different temperatures. Our novel, non-contact method for measuring creep uses rapid rotation to deform the sample. Digital images of the deformed samples are analyzed to obtain the creep properties, which match those obtained using ASTM Standard E-139 for Nb at 1985 C. Data from selected ESL-based characterization studies will be presented. The ESL technique could support numerous propulsion technologies by advancing the knowledge base and the technology readiness level for ultra-high temperature materials. Applications include non-eroding nozzle materials and lightweight, high-temperature alloys for turbines and structures.

Rogers, Jan↗

RECONSTRUCTING THE INITIAL DENSITY FIELD OF THE LOCAL UNIVERSE: METHODS AND TESTS WITH MOCK CATALOGS

Our research objective in this paper is to reconstruct an initial linear density field, which follows the multivariate Gaussian distribution with variances given by the linear power spectrum of the current cold dark matter model and evolves through gravitational instabilities to the present-day density field in the local universe. For this purpose, we develop a Hamiltonian Markov Chain Monte Carlo method to obtain the linear density field from a posterior probability function that consists of two components: a prior of a Gaussian density field with a given linear spectrum and a likelihood term that is given by the current density field. The present-day density field can be reconstructed from galaxy groups using the method developed in Wang et al. Using a realistic mock Sloan Digital Sky Survey DR7, obtained by populating dark matter halos in the Millennium simulation (MS) with galaxies, we show that our method can effectively and accurately recover both the amplitudes and phases of the initial, linear density field. To examine the accuracy of our method, we use N-body simulations to evolve these reconstructed initial conditions to the present day. The resimulated density field thus obtained accurately matches the original density field of the MS in the density range 0.3{approx}<{rho}/ {rho}-bar {approx}<20 without any significant bias. In particular, the Fourier phases of the resimulated density fields are tightly correlated with those of the original simulation down to a scale corresponding to a wavenumber of {approx}1 h Mpc{sup -1}, much smaller than the translinear scale, which corresponds to a wavenumber of {approx}0.15 h Mpc{sup -1}.

79 ASTRONOMY AND ASTROPHYSICS↗

Underway Sampling of Marine Inherent Optical Properties on the Tara Oceans Expedition as a Novel Resource for Ocean Color Satellite Data Product Validation

Developing and validating data records from operational ocean color satellite instruments requires substantial volumes of high quality in situ data. In the absence of broad, institutionally supported field programs, organizations such as the NASA Ocean Biology Processing Group seek opportunistic datasets for use in their operational satellite calibration and validation activities. The publicly available, global biogeochemical dataset collected as part of the two and a half year Tara Oceans expedition provides one such opportunity. We showed how the inline measurements of hyperspectral absorption and attenuation coefficients collected onboard the R/V Tara can be used to evaluate near-surface estimates of chlorophyll-a, spectral particulate backscattering coefficients, particulate organic carbon, and particle size classes derived from the NASA Moderate Resolution Imaging Spectroradiometer onboard Aqua (MODISA). The predominant strength of such flow-through measurements is their sampling rate-the 375 days of measurements resulted in 165 viable MODISA-to-in situ match-ups, compared to 13 from discrete water sampling. While the need to apply bio-optical models to estimate biogeochemical quantities of interest from spectroscopy remains a weakness, we demonstrated how discrete samples can be used in combination with flow-through measurements to create data records of sufficient quality to conduct first order evaluations of satellite-derived data products. Given an emerging agency desire to rapidly evaluate new satellite missions, our results have significant implications on how calibration and validation teams for these missions will be constructed.

remote sensing↗

Application of real rock pore-threat statistics to a regular pore network model

This work reports the application of real rock statistical data to a previously developed regular pore network model in an attempt to produce an accurate simulation tool with low computational overhead. A core plug from the St. Peter Sandstone formation in Indiana was scanned with a high resolution micro CT scanner. The pore-throat statistics of the three-dimensional reconstructed rock were extracted and the distribution of the pore-throat sizes was applied to the regular pore network model. In order to keep the equivalent model regular, only the throat area or the throat radius was varied. Ten realizations of randomly distributed throat sizes were generated to simulate the drainage process and relative permeability was calculated and compared with the experimentally determined values of the original rock sample. The numerical and experimental procedures are explained in detail and the performance of the model in relation to the experimental data is discussed and analyzed. Petrophysical properties such as relative permeability are important in many applied fields such as production of petroleum fluids, enhanced oil recovery, carbon dioxide sequestration, ground water flow, etc. Relative permeability data are used for a wide range of conventional reservoir engineering calculations and in numerical reservoir simulation. Two-phase oil water relative permeability data are generated on the same core plug from both pore network model and experimental procedure. The shape and size of the relative permeability curves were compared and analyzed and good match has been observed for wetting phase relative permeability but for non-wetting phase, simulation results were found to be deviated from the experimental ones. Efforts to determine petrophysical properties of rocks using numerical techniques are to eliminate the necessity of regular core analysis, which can be time consuming and expensive. So a numerical technique is expected to be fast and to produce reliable results. In applied engineering, sometimes quick result with reasonable accuracy is acceptable than the more time consuming results. Present work is an effort to check the accuracy and validity of a previously developed pore network model for obtaining important petrophysical properties of rocks based on cutting-sized sample data.

58 GEOSCIENCES↗

Application of real rock pore-throat statistics to a regular pore network model

This work reports the application of real rock statistical data to a previously developed regular pore network model in an attempt to produce an accurate simulation tool with low computational overhead. A core plug from the St. Peter Sandstone formation in Indiana was scanned with a high resolution micro CT scanner. The pore-throat statistics of the three-dimensional reconstructed rock were extracted and the distribution of the pore-throat sizes was applied to the regular pore network model. In order to keep the equivalent model regular, only the throat area or the throat radius was varied. Ten realizations of randomly distributed throat sizes were generated to simulate the drainage process and relative permeability was calculated and compared with the experimentally determined values of the original rock sample. The numerical and experimental procedures are explained in detail and the performance of the model in relation to the experimental data is discussed and analyzed. Petrophysical properties such as relative permeability are important in many applied fields such as production of petroleum fluids, enhanced oil recovery, carbon dioxide sequestration, ground water flow, etc. Relative permeability data are used for a wide range of conventional reservoir engineering calculations and in numerical reservoir simulation. Two-phase oil water relative permeability data are generated on the same core plug from both pore network model and experimental procedure. The shape and size of the relative permeability curves were compared and analyzed and good match has been observed for wetting phase relative permeability but for non-wetting phase, simulation results were found to be deviated from the experimental ones. Efforts to determine petrophysical properties of rocks using numerical techniques are to eliminate the necessity of regular core analysis, which can be time consuming and expensive. So a numerical technique is expected to be fast and to produce reliable results. In applied engineering, sometimes quick result with reasonable accuracy is acceptable than the more time consuming results. Present work is an effort to check the accuracy and validity of a previously developed pore network model for obtaining important petrophysical properties of rocks based on cutting-sized sample data. Introduction

42 ENGINEERING↗