Hydrofluoric Acid Reaction with a Porous Sandstone

1971 ◽  
Vol 11 (03) ◽  
pp. 306-314 ◽  
Author(s):  
B.B. Williams ◽  
M.E. Whiteley

Abstract A procedure is developed for predicting changes in the porosity distribution in a sandstone resulting from reaction with hydrofluoric acid. This procedure is based on a theory for slow heterogeneous reactions in a porous solid where the solid matrix is consumed in the reaction Process. Reaction-rate data for use in this theory are obtained from experiments where acid is injected through short cores and effluent concentration measured using a fluoride specific ion electrode. This rate is found to be first order in hydrofluoric acid concentration. Variations in rate with temperature and quantity of rock dissolved are shown. Introduction Mixtures of hydrofluoric and hydrochloric acid are used to stimulate gas and oil production from sandstone reservoirs by increasing formation porosity and permeability near the wellbore. This porosity and permeability near the wellbore. This acid will react with almost all constituents of naturally occurring sandstones, such as silica, feldspar, clays, and calcareous material. In order to utilize this acid effectively, it is necessary to predict where acid reacts and changes that occur predict where acid reacts and changes that occur with reaction. Chemical reactions between hydrofluoric acid and silica or calcite in the rock matrix are simple, well known reactions. However, reactions with silicates such as clays or feldspars are complex since these minerals occur as three-dimensional lattices with only average empirical formulas. Examples are kaolinite ([A Fe +3 Mg]-Si O1.8 0.1 0.1 2 5 [OH]. Ca), montmorillonite (A Mg Si4 0.05 1.67 0.33 4 O [OH] . NA), and feldspars such as albite10 2 0.33 ([NaSi A ]). In the reactions shown below the3 8 reaction of sodium silicate is used to represent hydrofluoric acid reaction with silicates found in the matrix. REACTION WITH SILICA SiO + 4HF SiF + 2H O2 4 2 SiF + 2HF H SiF4 2 6 REACTION WITH SILICATES (FELDSPAR OR CLAYS) Na SiO + 8HF SiF + 4NaF + 4H O4 4 4 2 2NaF + SiF Na SiF4 2 6 2HF + SiF H SiF4 2 6 REACTION WITH CALCITE CaCO + 2HF CaF + H O + CO3 2 2 2 Reaction of HF-HCl mixtures with silicates and quartz has been the subject of studies by Blumberg, Blumberg and Stavinou, Gatewood et al. and Smith and Hendrickson. These studies indicate that the reaction is first order in HF concentration and that reaction rate with silicates is at least 10 times faster than reaction with silica. To dare, a reliable method for relating reaction data taken on finely ground silica, dispersed clays, or glass slides to the acidization process in a sandstone formation has not been developed Figs. 1 and 2 are photomicrographs of a Berea sandstone core illustrating the system in which acid reaction occurs. In these photomicrographs silica grains are black, and a few feldspar grains are apparent because of their internal porosity, which gives a streaked appearance. Unfortunately, clay or calcite cannot be differentiated from the pore space since all appear as an area shading from pore space since all appear as an area shading from gray to white. It is apparent that the heterogeneous nature of the porous material greatly complicates the reaction problem. For this reason, a theory including mass problem. For this reason, a theory including mass transport, surface kinetics, and a statistical representation for the porous material is required to describe acid reaction. SPEJ P. 306

2009 ◽  
Vol 20 (04) ◽  
pp. 513-526 ◽  
Author(s):  
SHOUJU LI ◽  
YUEFANG WANG ◽  
YINGXI LIU ◽  
WEI SUN

The relationship between thermal conductivity of porous material and fractal dimension is numerically simulated by using the finite element method. The solid matrix and pore space are generated randomly according to material porosity. Material parameters and element properties are changed by using ANSYS parameter design language. The effective thermal conductivity is derived according to thermal flux through some sections computed by FEM and Fourier heat transform law. The investigation shows that the effective thermal conductivity decreases with increasing porosity. The effective thermal conductivity will decrease exponentially with increasing fractal dimension of porosity space and increase exponentially with increasing fractal dimension of solid matrix.


2014 ◽  
Vol 28 (2) ◽  
pp. 231-237 ◽  
Author(s):  
Lech W. Szajdak ◽  
Jerzy Lipiec ◽  
Anna Siczek ◽  
Artur Nosalewicz ◽  
Urszula Majewska

Abstract The aim of this study was to verify first-order kinetic reaction rate model performance in predicting of leaching of atrazine and inorganic compounds (K+1, Fe+3, Mg+2, Mn+2, NH4 +, NO3 - and PO4 -3) from tilled and orchard silty loam soils. This model provided an excellent fit to the experimental concentration changes of the compounds vs. time data during leaching. Calculated values of the first-order reaction rate constants for the changes of all chemicals were from 3.8 to 19.0 times higher in orchard than in tilled soil. Higher first-order reaction constants for orchard than tilled soil correspond with both higher total porosity and contribution of biological pores in the former. The first order reaction constants for the leaching of chemical compounds enables prediction of the actual compound concentration and the interactions between compound and soil as affected by management system. The study demonstrates the effectiveness of simultaneous chemical and physical analyses as a tool for the understanding of leaching in variously managed soils.


2021 ◽  
Vol 11 (6) ◽  
pp. 2495
Author(s):  
Belén Ferrer ◽  
María-Baralida Tomás ◽  
David Mas

Some materials undergo hygric expansion when soaked. In porous rocks, this effect is enhanced by the pore space, because it allows water to reach every part of its volume and to hydrate most swelling parts. In the vicinity, this enlargement has negative structural consequences as adjacent elements support some compressions or displacements. In this work, we propose a normalized cross-correlation between rock surface texture images to determine the hygric expansion of such materials. We used small porous sandstone samples (11 × 11 × 30 mm3) to measure hygric swelling. The experimental setup comprised an industrial digital camera and a telecentric objective. We took one image every 5 min for 3 h to characterize the whole swelling process. An error analysis of both the mathematical and experimental methods was performed. The results showed that the proposed methodology provided, despite some limitations, reliable hygric swelling information by a non-contact methodology with an accuracy of 1 micron and permitted the deformation in both the vertical and horizontal directions to be explored, which is an advantage over traditional linear variable displacement transformers.


2019 ◽  
Vol 292 ◽  
pp. 01063
Author(s):  
Lubomír Macků

An alternative method of determining exothermic reactor model parameters which include first order reaction rate constant is described in this paper. The method is based on known in reactor temperature development and is suitable for processes with changing quality of input substances. This method allows us to evaluate the reaction substances composition change and is also capable of the reaction rate constant (parameters of the Arrhenius equation) determination. Method can be used in exothermic batch or semi- batch reactors running processes based on the first order reaction. An example of such process is given here and the problem is shown on its mathematical model with the help of simulations.


1976 ◽  
Vol 29 (2) ◽  
pp. 443 ◽  
Author(s):  
MA Haleem ◽  
MA Hakeem

Kinetic data are reported for the decarboxylation of β-resorcylic acid in resorcinol and catechol for the first time. The reaction is first order. The observation supports the view that the decomposition proceeds through an intermediate complex mechanism. The parameters of the absolute reaction rate equation are calculated.


2015 ◽  
Vol 782 ◽  
pp. 260-299 ◽  
Author(s):  
Preyas N. Shah ◽  
Eric S. G. Shaqfeh

Surfaces that include heterogeneous mass transfer at the microscale are ubiquitous in nature and engineering. Many such media are modelled via an effective surface reaction rate or mass transfer coefficient employing the conventional ansatz of kinetically limited transport at the microscale. However, this assumption is not always valid, particularly when there is strong flow. We are interested in modelling reactive and/or porous surfaces that occur in systems where the effective Damköhler number at the microscale can be $O(1)$ and the local Péclet number may be large. In order to expand the range of the effective mass transfer surface coefficient, we study transport from a uniform bath of species in an unbounded shear flow over a flat surface. This surface has a heterogeneous distribution of first-order surface-reactive circular patches (or pores). To understand the physics at the length scale of the patch size, we first analyse the flux to a single reactive patch. We use both analytic and boundary element simulations for this purpose. The shear flow induces a 3-D concentration wake structure downstream of the patch. When two patches are aligned in the shear direction, the wakes interact to reduce the per patch flux compared with the single-patch case. Having determined the length scale of the interaction between two patches, we study the transport to a periodic and disordered distribution of patches again using analytic and boundary integral techniques. We obtain, up to non-dilute patch area fraction, an effective boundary condition for the transport to the patches that depends on the local mass transfer coefficient (or reaction rate) and shear rate. We demonstrate that this boundary condition replaces the details of the heterogeneous surfaces at a wall-normal effective slip distance also determined for non-dilute patch area fractions. The slip distance again depends on the shear rate, and weakly on the reaction rate, and scales with the patch size. These effective boundary conditions can be used directly in large-scale physics simulations as long as the local shear rate, reaction rate and patch area fraction are known.


SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 280-292 ◽  
Author(s):  
John Lyons ◽  
Hadi Nasrabadi ◽  
Hisham A. Nasr-El-Din

Summary Fracture acidizing is a well-stimulation technique used to improve the productivity of low-permeability reservoirs and to bypass deep formation damage. The reaction of injected acid with the rock matrix forms etched channels through which oil and gas can then flow upon production. The properties of these etched channels depend on the acid-injection rate, temperature, reaction chemistry, mass-transport properties, and formation mineralogy. As the acid enters the formation, it increases in temperature by heat exchange with the formation and the heat generated by acid reaction with the rock. Thus, the reaction rate, viscosity, and mass transfer of acid inside the fracture also increase. In this study, a new thermal-fracture-acidizing model is presented that uses the lattice Boltzmann method to simulate reactive transport. This method incorporates both accurate hydrodynamics and reaction kinetics at the solid/liquid interface. The temperature update is performed by use of a finite-difference technique. Furthermore, heterogeneity in rock properties (e.g., porosity, permeability, and reaction rate) is included. The result is a model that can accurately simulate realistic fracture geometries and rock properties at the pore scale and that can predict the geometry of the fracture after acidizing. Three thermal-fracture-acidizing simulations are presented here, involving injection of 15 and 28 wt% of hydrochloric acid into a calcite fracture. The results clearly show an increase in the overall fracture dissolution because of the addition of temperature effects (increasing the acid-reaction and mass-transfer rates). It has also been found that by introducing mineral heterogeneity, preferential dissolution leads to the creation of uneven etching across the fracture surfaces, indicating channel formation.


2008 ◽  
Vol 5 (3) ◽  
pp. 598-606
Author(s):  
Vandana Sharma ◽  
K. V. Sharma ◽  
V. W. Bhagwat

The kinetics and mechanism of cetyltrimethylammonium bromide catalyzed oxidation of tetraethylene glycol [2,2'-(oxibis(ethylenoxy)diethanol)] byN-chlorosaccharin in aqueous acetic acid medium in presence of perchloric acid have been investigated at 323K. The reaction is first order dependence on Nchlorosaccharin. The reaction rate follows first order kinetics with respect to [tetraethylene glycol] with excess concentration of other reactants. The miceller effect due to cetyltrimethylammonium bromide, a cationic surfactant has been studied. The change in ionic strength shows negligible salt effect. The dielectric effect is found to be positive. Addition of one of the products (saccharin) retards the reaction rate. Activation parameters are calculated from the Arrhenious plot. A possible mechanism consistent with the experimental results has been proposed.


2016 ◽  
Vol 16 (3) ◽  
pp. 1245-1254 ◽  
Author(s):  
T. P. Riedel ◽  
Y.-H. Lin ◽  
Z. Zhang ◽  
K. Chu ◽  
J. A. Thornton ◽  
...  

Abstract. Isomeric epoxydiols from isoprene photooxidation (IEPOX) have been shown to produce substantial amounts of secondary organic aerosol (SOA) mass and are therefore considered a major isoprene-derived SOA precursor. Heterogeneous reactions of IEPOX on atmospheric aerosols form various aerosol-phase components or "tracers" that contribute to the SOA mass burden. A limited number of the reaction rate constants for these acid-catalyzed aqueous-phase tracer formation reactions have been constrained through bulk laboratory measurements. We have designed a chemical box model with multiple experimental constraints to explicitly simulate gas- and aqueous-phase reactions during chamber experiments of SOA growth from IEPOX uptake onto acidic sulfate aerosol. The model is constrained by measurements of the IEPOX reactive uptake coefficient, IEPOX and aerosol chamber wall losses, chamber-measured aerosol mass and surface area concentrations, aerosol thermodynamic model calculations, and offline filter-based measurements of SOA tracers. By requiring the model output to match the SOA growth and offline filter measurements collected during the chamber experiments, we derive estimates of the tracer formation reaction rate constants that have not yet been measured or estimated for bulk solutions.


2011 ◽  
Vol 312-315 ◽  
pp. 364-369 ◽  
Author(s):  
Seyed Amir Bahrani ◽  
Catherine Loisel ◽  
Jean Yves Monteau ◽  
Sid Ahmed Rezzoug ◽  
Zoulikha Maache-Rezzoug

Two hydrothermal treatment processes (DV-HMT and DIC treatment) were investigated on standard maize starch for three processing temperatures; 100, 110 and 120°C. The gravimetric change of starch powder during the treatment was analyzed by a simultaneous water diffusion and starch reaction model. The effective diffusivity coefficient (Deff) and reaction rate constant (k) were estimated by minimizing the error between experimental and analytical results. The values of Deff and k clearly increased with temperature. The degree of starch melting was evaluated for the two treatments using the first-order reaction model as a function of processing time. The results suggest that the absorption process is controlled by water–starch reactivity that induces melting phenomenon of starch crystallites, which progresses when temperature increases. The two hydrothermal treatments considerably differ: DIC being more prone to water absorption as demonstrated by the values of Deff and k.


Sign in / Sign up

Export Citation Format

Share Document