包气带土壤对石油烃的截留作用研究

包气带土壤对地表污染物的截留作用是地下水免受污染的一道天然屏障.土壤截留污染物的能力对地下水石油烃污染的治理起着关键作用,在对实际污染场地进行调查研究的基础上,采用土柱淋滤、吸附/解吸、石油挥发等室内模拟实验,研究包气带土壤对石油烃的截留作用及其影响因素.结果表明,细砂、中砂和粗砂3种土壤对石油烃的截留率分别为81.0...     

Investigation of the vadose zone using barometric pressure cycles

This paper documents a technique for investigating one-dimensional airflow in the vadose zone. Variations in pore gas pressures resulting from barometric cycles were measured at depths as great as 180 m in several gas monitoring wells. The data were transformed to the frequency domain, enabling comparison with closed-form analytic expressions of one-dimensional pressure transport in layered porous media. The data reveal evidence for vertical fracture flow that was not apparent from in situ measurements of permeability. The data also reveal that the basalt underlying the site at depths greater than 100 m has permeability exceeding 1000 darcies, and is vented to the atmosphere at an estimated distance of a few kilometers from the site.     

Plant aided bioremediation in the vadose zone: model development and applications

Phytoremediation has the potential to enhance clean up of land contaminated by various pollutants. A mathematical model that includes a two-fluid phase flow model of water flow as well as a two-region soil model of contaminant reactions was developed and applied to various bioremediation scenarios in the unsaturated zone, especially to plant-aided bioremediation. To investigate model behavior and determine the main parameters and mechanisms that affect bioremediation in unplanted and planted soils, numerical simulations of theoretical scenarios were conducted before applying the model to field data. It is observed from the results that parameters affecting the contaminant concentration in the water phase, such as aqueous solubility, the octanol-water partition coefficient, and organic carbon content of the soil controlled the contaminant fate in the vadose zone. Simulation using the developed model also characterized the fate and transport of the contaminants both in planted and unplanted soils satisfactorily for field applications. Although phytoremediation has the potential for remediation of contaminated soils, results from both modeling and field studies suggested that plants may not always enhance the remediation efficiency when the soil already has a high microbial concentration, when the contaminant bioavailability is low, or when the overall reaction is mass transfer-limited. Therefore, other steps to increase contaminant bioavailability are needed in phytoremediation applications; natural purification mechanisms such as aging, volatilization, and natural bioremediation should be considered to maximize the plant effect and minimize the cost.     

Two-dimensional laboratory simulation of LNAPL infiltration and redistribution in the vadose zone

A quantitative two-dimensional laboratory experiment was conducted to investigate the immiscible flow of a light non-aqueous phase liquid (LNAPL) in the vadose zone. An image analysis technique was used to determine the two-dimensional saturation distribution of LNAPL, water and air during LNAPL infiltration and redistribution. Vertical water saturation variations were also continuously monitored with miniature resistivity probes. LNAPL and water pressures were measured using hydrophobic and hydrophilic tensiometers. This study is limited to homogeneous geological conditions, but the unique experimental methods developed will be used to examine more complex systems. The pressure measurements and the quantification of the saturation distribution of all the fluids in the entire flow domain under transient conditions provide quantitative data essential for testing the predictive capability of numerical models. The data are used to examine the adequacy of the constitutive pressure-saturation relations that are used in multiphase flow models. The results indicate that refinement of these commonly used hydraulic relations is needed for accurate model prediction. It is noted in particular that, in three-fluid phase systems, models should account for the existence of a residual NAPL saturation occurring after NAPL drainage. This is of notable importance because residual NAPL can act as a non negligible persistent source of contamination.     

An analytical solution for vertical transport of volatile chemicals in the vadose zone

An analytical solution is presented for one-dimensional vertical transport of volatile chemicals through the vadose zone to groundwater. The solution accounts for the important transport mechanisms of the steady advection of water and gas, diffusion and dispersion in water and gas, as well as adsorption, and first-order degradation. By assuming a linear, equilibrium partitioning between water, gas and the adsorbed chemical phases, the dependent variable in the mathematical model becomes the total resident concentration. The general solution was derived for cases having a constant initial total concentration over a discrete depth interval and a zero initial total concentration elsewhere. A zero concentration gradient is assumed at the groundwater table. Examples are given to demonstrate the application of the new solution for calculating the case of a non-uniform initial source concentration, and estimating the transport of chemicals to the groundwater and the atmosphere. The solution was also used to verify a numerical code called VLEACH. We discovered an error in VLEACH, and found that the new solution agreed very well with the numerical results by corrected VLEACH. A simplified solution to predict the migration of volatile organic chemical due to the gas density effect has shown that a high source concentration may lead to significant downward advective gas-phase transport in a soil with a high air-permeability.     

基于DRASTIC的包气带阻滞污染物能力研究

对包气带阻滞污染物能力的研究,对保护地下水脆弱性具有重要意义.将传统DRASTIC方法结合层次分析法,对各评价指标赋以新的权重,并依托地理信息系统(GIS)平台,对河南省土壤包气带阻滞污染物的能力进行评价,得到了河南省包气带阻滞能力分区图.结果表明,针对现状所进行的包气带阻滞污染物能力的评价对实际情况具有一定的代表性,...     

Transport and degradation of propylene glycol in the vadose zone: model development and sensitivity analysis

Transport and degradation of de-icing chemical (containing propylene glycol, PG) in the vadose zone were studied with a lysimeter experiment and a model, in which transient water flow, kinetic degradation of PG and soil chemistry were combined. The lysimeter experiment indicated that aerobic as well as anaerobic degradation occurs in the vadose zone. Therefore, the model included both types of degradation, which was made possible by assuming advection-controlled (mobile) and diffusion-controlled (immobile) zones. In the mobile zone, oxygen can be transported by diffusion in the gas phase. The immobile zone is always water-saturated, and oxygen only diffuses slowly in the water phase. Therefore, the model is designed in a way that the redox potential can decrease when PG is degraded, and thus, anaerobic degradation can occur. In our model, manganese oxide (MnO₂, which is present in the soil) and NO \(_{3}^{-}\) (applied to enhance biodegradation) can be used as electron acceptors for anaerobic degradation. The application of NO \(_{3}^{-}\) does not result in a lower leaching of PG nor in a slower depletion of MnO₂. The thickness of the snowcover influences the leached fraction of PG, as with a high infiltration rate, transport is fast, there is less time for degradation and thus more PG will leach. The model showed that, in this soil, the effect of the water flow dominates over the effect of the degradation parameters on the leaching at a 1-m depth.     

Assessing tungsten transport in the vadose zone: from dissolution studies to soil columns

This study investigates the dissolution, sorption, leachability, and plant uptake of tungsten and alloying metals from canister round munitions in the presence of model, well characterized soils. The source of tungsten was canister round munitions, composed mainly of tungsten (95%) with iron and nickel making up the remaining fraction. Three soils were chosen for the lysimeter studies while four model soils were selected for the adsorption studies. Lysimeter soils were representatives of the typical range of soils across the continental USA; muck-peat, clay-loamy and sandy-quartzose soil. Adsorption equilibrium data on the four model soils were modeled with Langmuir and linear isotherms and the model parameters were obtained. The adsorption affinity of soils for tungsten follows the order: Pahokee peat > kaolinite > montmorillonite > illite. A canister round munition dissolution study was also performed. After 24 d, the measured dissolved concentrations were: 61.97, 3.56, 15.83 mg L⁻¹ for tungsten, iron and nickel, respectively. Lysimeter transport studies show muck peat and sandy quartzose soils having higher tungsten concentration, up to 150 mg kg⁻¹ in the upper layers of the lysimeters and a sharp decline with depth suggesting strong retardation processes along the soil profile. The concentrations of tungsten, iron and nickel in soil lysimeter effluents were very low in terms of posing any environmental concern; although no regulatory limits have been established for tungsten in natural waters. The substantial uptake of tungsten and nickel by ryegrass after 120 d of exposure to soils containing canister round munition suggests the possibility of tungsten and nickel entering the food chain.     

Natural attenuation/phytoremediation in the vadose zone of a former industrial sludge basin

The natural attenuation of polyaromatic hydrocarbons (PAHs) in the vadose zone of a naturally revegetated former industrial sludge basin (0.45 ha) was examined. This was accomplished by comparing the concentration of 16 PAH contaminants present in sludge collected below the root zone of plants with contaminants present at 3 shallower depths within the root zone. Chemical analysis of 240 samples from 60 cores showed the average concentration of total and individual PAHs in the 0-30 cm, 30-60 cm, and bottom of the root zone strata were approximately 10, 20, and 50%, respectively, of the 16, 800 ppm average total PAH concentration in deep non-rooted sludge. Statistically significant differences in average PAH concentrations were observed between each strata studied and the non-rooted sludge except for the concentrations of acenaphthene and chrysene present at the bottom of the root zone in comparison to sludge values. The rooting depth of the vegetation growing in the basin was dependent on both vegetation type and plant age. Average rooting depths for trees, forbs (herbaceous non-grasses), and grasses were 90, 60, and 50 cm, respectively. The deepest root systems observed (100-120 cm) were associated with the oldest (12-14 year-old) mulberry trees. Examination of root systems and PAH concentrations at numerous locations and depths within the basin indicated that plant roots and their microbially active rhizospheres fostered PAH disappearance; including water insoluble, low volatility compounds, i.e. benzo(a)pyrene and benzo(ghi)perylene. The reduced concentration of PAHs in the upper strata of this revegetated former sludge basin indicated that natural attenuation had occurred. This observation supports the concept that through appropriate planting and management practices (phytoremediation) it will be possible to accelerate, maximize, and sustain natural processes, whereby even the most recalcitrant PAH contaminants (i.e. benzo(a)pyrene) can be remediated over time.     

Effects of soil moisture and physical-chemical properties of organic pollutants on vapor-phase transport in the vadose zone

Vapor-phase transport of organic pollutants is one of the important pathways in the distribution and attenuation of volatile organic compounds in the vadose zone. In this study, the impact of vapor-phase partitioning and of the physical-chemical properties of organic pollutants on vapor-phase transport was assessed. An experimentally derived relationship to predict vapor sorption for a variety of soil types under varying soil moisture conditions was incorporated into the two-dimensional finite-element model, Vocwaste. The revised model was then used to simulate the transport of volatile organics. Vapor-phase partitioning in the model accounted for vapor uptake by sorption onto moist mineral surfaces as well as sorption at the liquid-solid interface and dissolution into soil water. Under dry conditions, vapor-phase sorption of volatile organic pollutants was shown to have a retarding effect on transport of organic vapors. However, for shallow, contaminated soils, volatilization was controlled by vapor diffusion, even under dry conditions where vapor-phase sorption was high. The influence of Henry's law constant and of the aqueous-phase (solid-liquid) partition coefficient for volatile organic pollutants was considered in the simulations. Volatilization of organic vapors was shown to be favored for contaminants with high Henry's law constants and low aqueous-phase partitioning coefficients. Because of the interdependence of these two physical-chemical properties, individual properties of the contaminant should not be considered in isolation in the evaluation of vapor transport.     

Estimation of primary drainage three-phase relative permeability for organic liquid transport in the vadose zone

The modeling of transport of organic liquid contaminants through the vadose zone often requires three-phase relative permeabilities. Since these are difficult to measure, predictive models are usually used. The objective of this study is to assess the ability of eight common models to predict the drainage relative permeability to oil in a three-phase system (water-oil-air). A comparison of the models' estimates using data set from Oak [Oak, M.J., 1990. Three-phase relative permeability of water-wet Berea. In: Seventh Symposium on Enhanced Oil Recovery, Paper SPE/Doe 20183. Tulsa, OK, April 22-25] showed that they provide very different predictions for the same system. The goodness of the models does not increase with the amount of data or computation that the models require. Also, the calculations showed how different interpretations of the models and of the terminology associated with them can significantly impact the predictions. Thus, considerable error may be introduced into the simulations of organic liquid transport in the vadose zone depending on the selection and interpretation of the three-phase relative permeability model.     

The application of ground penetrating radar attenuation tomography in a vadose zone infiltration experiment

Cross-borehole ground penetrating radar (XBGPR) is used in monitoring a long-term vadose zone infiltration experiment at a test site in Socorro, NM in order to examine contaminant transport in the vadose zone. XBGPR attenuation tomography is conducted in order to test the ability of using images of electromagnetic attenuation for hydrogeologic investigations. The results of four pre-infiltration attenuation inversions shows standard deviations below 0.1 Np/m, and demonstrate the consistency of the XBGPR tomography technique for making time-lapse observations. Correlation to the core records indicates that XBGPR attenuation tomograms provide high-resolution images of clay distribution in the vadose zone. Water infiltration at the ground surface was initiated in February 1999 at a constant rate of 2.7 cm/day, and continued at this rate throughout the data collection experiment. Time-lapse attenuation tomograms show that attenuation increases by approximately 0.3 Np/m during the water infiltration, and indicate a snowplow effect may be occurring where salts are dissolved by the water and concentrated at the front of the plume. Seasonal temperature changes may also cause changes in electromagnetic attenuation images, and masking the evidence of water infiltration. Thus caution must be taken when using time-lapse attenuation images to interpret the movement of a water plume during a long-term experiment as temperature changes.     

Coupled effects of solution chemistry and hydrodynamics on the mobility and transport of quantum dot nanomaterials in the vadose zone

To investigate the coupled effects of solution chemistry and hydrodynamics on the mobility of quantum dot (QD) nanoparticles in the vadose zone, laboratory scale transport experiments involving single and/or sequential infiltrations of QDs in unsaturated and saturated porous media, and computations of total interaction and capillary potential energies were performed. As ionic strength increased, QD retention in the unsaturated porous media increased; however, this retention was significantly suppressed in the presence of a non-ionic surfactant in the infiltration suspensions as indicated by surfactant enhanced transport of QDs. In the vadose zone, the non-ionic surfactant limited the formation of QD aggregates, enhanced QD mobility and transport, and lowered the solution surface tension, which resulted in a decrease in capillary forces that not only led to a reduction in the removal of QDs, but also impacted the vadose zone flow processes. When chemical transport conditions were favorable (ionic strength of 5 × 10(-4)M and 5 × 10(-3)M, or ionic strengths of 5 × 10(-2)M and 0.5M with surfactant), the dominating phenomena controlling the mobility and transport of QDs in the vadose zone were meso-scale processes, where infiltration by preferential flow results in the rapid transport of QDs. When chemical transport conditions were unfavorable (ionic strength of 5 × 10(-2)M and 0.5M) the dominating phenomena controlling the mobility and transport of QDs in the vadose zone were pore-scale processes governed by gas-water interfaces (GWI) that impact the mobility of QDs. The addition of surfactant enhanced the transport of QDs both in favorable and unfavorable chemical transport conditions. The mobility and retention of QDs was controlled by interaction and capillary forces, with the latter being the most influential. GWI were found to be the dominant mechanism and site for QD removal compared with solid-water interfaces (SWI) and pore straining. Additionally, ripening phenomena were demonstrated to enhance QDs removal or retention in porous media and to be attenuated by the presence of surfactant.     

Metabolism of the explosive hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in a contaminated vadose zone

The aim of this study was to explore biodegradation potential of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in a deep contaminated unsaturated zone over Israel's coastal aquifer. While anaerobic biodegradation potential was observed throughout the profile down to the water table at a depth of 45 m, aerobic biodegradation was limited to the surface of the unsaturated zone. Traces of nitroso-RDX intermediates were detected in the soil samples, indicating possible in situ activity. Polymerase chain reaction and denaturing gradient gel electrophoresis analysis revealed that the microbial population in the soil consisted of protobacteria, but no known RDX degraders were detected. However, a 16S rRNA gene sequence most similar to Sphingomonas sp. was detected at all depths. Biodegradation rates were faster in the surface (0 and 1m) versus deeper soil samples (22 and 45 m) and were not affected under anaerobic conditions by the presence of nitrate, indicating a concurrent reduction of both compounds. RDX half-life in the surface soil was mostly dependent on carbon content and to lesser extent on soil moisture. Biomineralization of RDX to CO(2) was confirmed by incubating surface soil with (14)C-labeled RDX. An aerobic RDX-degrading bacterium, identified as Gordonia sp., was isolated from the soil: it degraded RDX aerobically and produced 4-nitro-2,4-diazabutanal. This study, the first to explore RDX biodegradation in the deep vadoze zone, indicates biodegradation potential throughout the profile, which is likely to support natural attenuation.     

Assessing in situ mineralization of recalcitrant organic compounds in vadose zone sediments using delta13C and 14C measurements

Few techniques exist to measure the biodegradation of recalcitrant organic compounds such as chlorinated hydrocarbons (CHC) in situ, yet predictions of biodegradation rates are needed for assessing monitored natural attenuation. Traditional techniques measuring O2, CO2, or chemical concentrations (in situ respiration, metabolite and soil air monitoring) may not be sufficiently sensitive to estimate biodegradation rates for these compounds. This study combined isotopic measurements (14C and delta13C of CO2 and delta13C of CHCs) in conjunction with traditional methods to assess in situ biodegradation of perchloroethylene (PCE) and its metabolites in PCE-contaminated vadose zone sediments. CHC, ethene, ethane, methane, O2, and CO2 concentrations were measured over 56 days using gas chromatography (GC). delta13C of PCE, trichloroethylene (TCE) and cis-1,2-dichloroethylene (DCE), delta13C and 14C of vadose zone CO2 and sediment organic matter, and delta13C, 14C, and deltaD of methane were measured using a GC-isotope ratio mass spectrometer or accelerator mass spectrometer. PCE metabolites accounted for 0.2% to 18% of CHC concentration suggesting limited reductive dechlorination. Metabolites TCE and DCE were significantly enriched in (13)C with respect to PCE indicating metabolite biodegradation. Average delta13C-CO2 in source area wells (-23.5 per thousand) was significantly lower compared to background wells (-18.4 per thousand) indicating CHC mineralization. Calculated CHC mineralization rates were 0.003 to 0.01 mg DCE/kg soil/day based on lower 14C values of CO2 in the contaminated wells (63% to 107% modern carbon (pMC)) relative to the control well (117 pMC). Approximately 74% of the methane was calculated to be derived from in situ CHC biodegradation based on the 14C measurement of methane (29 pMC). 14C-CO2 analyses was a sensitive measurement for quantifying in situ recalcitrant organic compound mineralization in vadose zone sediments for which limited methodological tools exist.     

Accelerated transport of (90)Sr following a release of high ionic strength solution in vadose zone sediments

Numerical simulation of cation exchange and mineral precipitation/dissolution reactions using the multiphase reactive geochemical transport code TOUGHREACT has provided important insight into the distribution of (90)Sr among layers of geologic strata in a complex vadose zone at the U. S. Department of Energy's Idaho National Laboratory. During a transfer operation in November 1972, 70.4 m(3) of acidic, high ionic strength liquid containing 15,900 Ci of (90)Sr was released over five days into alluvial gravels 137 m above the Snake River Plain Aquifer. Sampling data from perched water zones 33 m below the release contain very high levels of (90)Sr as do soil samples obtained nearer the point of release. Use of traditional simulation approaches using laboratory-measured constant partitioning coefficients (K(d)) cannot simultaneously explain perched water and soil concentrations. To address the discrepancy, a reactive transport approach was adopted to include competitive cation exchange, dissolution/precipitation of calcite, carbon dioxide gas production and transport, and gibbsite precipitation. Simulation results using this model suggest that some of the (90)Sr could have been transported very rapidly immediately after the release with the acceleration facilitated by competition for cation exchange sites with high sodium concentrations in the released liquid and calcium dissolved from calcite, and to a lesser extent by formation of aqueous complexes with nitrate. Once the leading edge of the liquid assemblage was flushed from the alluvium, the mobility of the remaining (90)Sr decreased significantly in the absence of the competing cations. Calculations indicate that there should be a net increase in calcite, suggesting that (90)Sr could be entrained in the mineral lattice, but insufficient field data exist for confirmation. Sensitivity studies show that the cation exchange selectivity coefficients were the most sensitive individual parameters determining the (90)Sr distribution. However, the most sensitive overall quantity was the total ion exchange capacity which is a function of the moles of exchange sites per volume of pore water, the cation exchange capacity, and the total volume wetted by the infiltrating solution. In contrast, the future mobility of (90)Sr was found to be relatively insensitive to the normal range in the composition of influxing precipitation and anthropogenic waters.     

In situ bioremediation of nitrate and perchlorate in vadose zone soil for groundwater protection using gaseous electron donor injection technology.

When present in the vadose zone, potentially toxic nitrate and perchlorate anions can be persistent sources of groundwater contamination. Gaseous electron donor injection technology (GEDIT), an anaerobic variation of petroleum hydrocarbon bioventing, involves injecting electron donor gases, such as hydrogen or ethyl acetate, into the vadose zone, to stimulate biodegradation of nitrate and perchlorate. Laboratory microcosm studies demonstrated that hydrogen and ethanol promoted nitrate and perchlorate reduction in vadose zone soil and that moisture content was an important factor. Column studies demonstrated that transport of particular electron donors varied significantly; ethyl acetate and butyraldehyde were transported more rapidly than butyl acetate and ethanol. Nitrate removal in the column studies, up to 100%, was best promoted by ethyl acetate. Up to 39% perchlorate removal was achieved with ethanol and was limited by insufficient incubation time. The results demonstrate that GEDIT is a promising remediation technology warranting further validation.     

A proton transfer reaction mass spectrometry based system for determining plant uptake of volatile organic compounds

In order to evaluate the contribution that higher plants make to the removal of volatile organic compounds from the atmosphere, a measurement system consisting of a proton transfer reaction mass spectrometer (PTR-MS), CO₂ analyzer, diffusion devise and leaf enclosure was established. The uptake of VOCs by Golden Pothos (Epipremnum aureum) was investigated. The overall relative error associated with measurements made using this system was <2.2% when a Golden Pothos leaf was exposed to 75–750 ppbv of methyl isobutyl ketone (MIBK). Even at the lowest MIBK concentration, more than 2.2% of the inflowing VOC was lost to the leaf, representing a detectable and positive MIBK uptake rate by the plant. The results of the investigation were compared with a measurement system based on gas chromatography analysis and it was shown that the use of a PTR-MS based system can significantly increase the certainties in determining the rate of VOC uptake by plants.     

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.     

Heat transfer analysis of cylindrical anaerobic reactors with different sizes: a heat transfer model

The temperature is the essential factor that influences the efficiency of anaerobic reactors. During the operation of the anaerobic reactor, the fluctuations of ambient temperature can cause a change in the internal temperature of the reactor. Therefore, insulation and heating measures are often used to maintain anaerobic reactor’s internal temperature. In this paper, a simplified heat transfer model was developed to study heat transfer between cylindrical anaerobic reactors and their surroundings. Three cylindrical reactors of different sizes were studied, and the internal relations between ambient temperature, thickness of insulation, and temperature fluctuations of the reactors were obtained at different reactor sizes. The model was calibrated by a sensitivity analysis, and the calibrated model was well able to predict reactor temperature. The Nash-Sutcliffe model efficiency coefficient was used to assess the predictive power of heat transfer models. The Nash coefficients of the three reactors were 0.76, 0.60, and 0.45, respectively. The model can provide reference for the thermal insulation design of cylindrical anaerobic reactors.