Hydrodynamic resistance of pore–throat structures and its effect on shale oil apparent permeability

Oil transport is greatly affected by heterogeneous pore–throat structures present in shale.It is therefore very important to accurately characterize pore–throat structures.Additionally,it remains unclear how pore–throat structures affect oil transport capacity.In this paper,using finite element(FE)simulation and mathematical modeling,we calculated the hydrodynamic resistance for four pore–throat structure.In addition,the influence of pore throat structure on shale oil permeability is analyzed.According to the results,the hydrodynamic resistance of different pore throat structures can vary by 300%.The contribution of additional resistance caused by streamline bending is also in excess of 40%,even without slip length.Fur-thermore,Pore–throat structures can affect apparent permeability by more than 60%on the REV scale,and this influence increases with heterogeneity of pore size distribution,organic matter content,and organic matter number.Clearly,modeling shale oil flow requires consideration of porous–throat structure and additional resistance,otherwise oil recovery and flow capacity may be overestimated.

Modeling of multiphase flow in low permeability porous media:Effect of wettability and pore structure properties

Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.

Postconditioning of sevoflurane and propofol is associated with mitochondrial permeability transition pore

Background: Sevoflurane and propofol are effective cardioprotective anaesthetic agents, though the cardioprotection of propofol has not been shown in humans. Their roles and underlying mechanisms in anesthetic postconditioning are unclear. Mitochondrial permeability transition pore (MPTP) opening is a major cause of ischemia-reperfusion injury. Here we investigated sevoflurane- and propofol-induced postconditioning and their relationship with MPTP. Methods: Isolated perfused rat hearts were exposed to 40 min of ischemia followed by 1 h of reperfusion. During the first 15 min of reperfusion, hearts were treated with either control buffer (CTRL group) or buffer containing 20 μmol/L atractyloside (ATR group), 3% (v/v) sevoflurane (SPC group), 50 μmol/L propofol (PPC group), or the combination of atractyloside with respective anesthetics (SPC+ATR and PPC+ATR groups). Infarct size was determined by dividing the total necrotic area of the left ventricle by the total left ventricular slice area (percent necrotic area). Results: Hearts treated with sevoflurane or propofol showed significantly better recovery of coronary flow, end-diastolic pressures, left ventricular developed pressure and derivatives compared with controls. Sevoflurane resulted in more protective alteration of hemodynamics at most time point of reperfusion than propofol. These improvements were paralleled with the reduction of lactate dehydrogenase release and the decrease of infarct size (SPC vs CTRL: (17.48±2.70)% vs (48.47±6.03)%, P<0.05; PPC vs CTRL: (35.60±2.10)% vs (48.47±6.03)%, P<0.05). SPC group had less infarct size than PPC group (SPC vs PPC: (17.48±2.70)% vs (35.60±2.10)%, P<0.05). Atractyloside coadministration attenuated or completely blocked the cardioprotective effect of postconditioning of sevoflurane and propofol. Conclusion: Postconditioning of sevoflurane and propofol has cardio-protective effect against ischemia-reperfusion injury of heart, which is associated with inhibition of MPTP opening. Compared to propofol, sevoflurane provides superior protection of functional recovery and infarct size.

Study of the permeability characteristics of porous media with methane hydrate by pore network model

The permeability in the methane hydrate reservoir is one of the key parameters in estimating the gas production performance and the flow behavior of gas and water during dissociation.In this paper,a three-dimensional cubic pore-network model based on invasion percolation is developed to study the effect of hydrate particle formation and growth habit on the permeability.The variation of permeability in porous media with different hydrate saturation is studied by solving the network problem.The simulation results are well consistent with the experimental data.The proposed model predicts that the permeability will reduce exponentially with the increase of hydrate saturation,which is crucial in developing a deeper understanding of the mechanism of hydrate formation and dissociation in porous media.

The calcium uniporter regulates the permeability transition pore in isolated cortical mitochondria

To investigate the influence of the mitochondrial calcium uniporter on the mitochondrial permeability transition pore, the present study observed mitochondrial morphology in cortical neurons isolated from adult rats using transmission electron microscopy, and confirmed the morphology and activity of isolated mitochondria by detecting succinic dehydrogenase and monoamine oxidase, two mitochondrial enzymes. Isolated mitochondria were treated with either ruthenium red, an inhibitor of the uniporter, spermine, an activator of the uniporter, or in combination with cyclosporin A, an inhibitor of the mitochondrial permeability transition pore. Results showed that ruthenium red inhibited CaCl2-induced mitochondrial permeability transition pore opening, spermine enhanced opening, and cyclosporin A attenuated the effects of spermine. Results demonstrated that the mitochondrial calcium uniporter plays a role in regulating the mitochondrial permeability transition pore in mitochondria isolated from the rat brain cortex.

Ethanol induced mitochondria injury and permeability transition pore opening: Role of mitochondria in alcoholic liver disease

AIM: To observe changes of mitochondria and investigate the effect of ethanol on mitochondrial perme- ability transition pore (PTP), mitochondrial membrane potential (MMP, ΔΨm) and intracellular calcium concentration in hepatocytes by establishing an animal model of alcoholic liver disease (ALD). METHODS: Fourty adult male Wistar rats were randomly divided into two groups, the model group (20) was administered alcohol intragastrically plus an Oliver oil diet to establish an ALD model, and the control group (20) was given an equal amount of normal saline. The ultramicrostructural changes of mitochondria were observed under electron microscopy. Mitochondria of liver was extracted, and patency of PTP, mitochondrial membrane potential (ΔΨm), mitochondrial mass and intracellular calcium concentration of isolated hepacytes were detected by flow cytometry using rhodamine123 (Rh123), Nonyl-Acridine Orange and calcium fluorescent probe Fluo-3/AM, respectively. RESULTS: Membrane and cristae were broken or disappeared in mitochondria in different shapes under electron microscopy. Some mitochondria showed U shape or megamitochondrion. In the model group, liver mitochondria PTP was broken, and mitochondria swelled, the absorbance at 450 nm, A540 decreased (0.0136 ± 0.0025 vs 0.0321 ± 0.0013, model vs control, P < 0.01); mitochondria transmembrane potential (239.4638 ± 12.7263 vs 377.5850 ± 16.8119, P < 0.01) was lowered; mitochondrial mass (17.4350 ± 1.9880 vs 31.6738 ± 3.4930, P < 0.01); and [Ca2+]i was increased in liver cells (7.0020 ± 0.5008 vs 10.2050 ± 0.4701, P < 0.01).CONCLUSION: Chronic alcohol intake might lead to broken mitochondria PTP, decreased mitochondria membrane potential and injury, and elevated intracellular Ca2+ production. Ethanol-induced chondriosome injury may be an important mechanism of alcoholic diseases.

An improved apparent permeability model considering full pore pressure range,variable intrinsic permeability and slippage coefficient

Although the slippage effect has been extensively studied,most of the previous studies focused on the impact of the slippage effect on apparent permeability within a low pore pressure range,resulting in the inability of matching the evolution of permeability in the remaining pressure range.In this paper,a new apparent permeability model that reveals the evolution of permeability under the combined action of effective stress and slippage in the full pore pressure range was proposed.In this model,both intrinsic permeability and slippage coefficient are stress dependent.Three experimental tests with pore pressure lower than 2 MPa and a test with pore pressure at about 10 MPa using cores from the same origin under constant confining stress and constant effective stress are conducted.By comparing experimental data and another apparent permeability model,we proved the fidelity of our newly developed model.Furthermore,the contribution factor of the slippage effect Rslip is used to determine the low pore pressure limit with significant slippage effect.Our results show that both narrow initial pore size and high effective stress increase the critical pore pressure.Finally,the evolutions of the slippage coefficient and the intrinsic permeability under different boundary conditions were analyzed.

Pore Structure and Permeability Characterization of Highrank Coal Reservoirs: A Case of the Bide-Santang Basin, Western Guizhou, South China

The methods of nuclear magnetic resonance(NMR)spectroscopy,mercury injection porosimetry(MIP),and gas-water relative permeability(GWRP)were used to reveal the pore structure and permeability characteristics of high-rank coal reservoirs in the Bide-Santang basin,western Guizhou,South China,to provide guidance for coalbed methane(CBM)exploration and exploitation and obtain direct insights for the development of CBM wells.The results indicate that the coal reservoirs in the study area are characterized by well-developed adsorption pores and poorly developed seepage pores.The bimodal NMR transverse relaxation time(T2)spectra and the mutation in the fractal characteristic of the MIP pore volume indicate poor connectivity between the adsorption pores and the seepage pores.As a result,the effective porosity is relatively low,with an average of 1.70%.The irreducible water saturation of the coal reservoir is relatively high,with an average of 66%,leading to a low gas relative permeability under irreducible water saturation.This is the main reason for the low recovery of high-rank CBM reservoirs,and effective enhanced CBM recovery technology urgently is needed.As a nondestructive and less time-consuming technique,the NMR is a promising method to quantitatively characterize the pores and fractures of coals.

Effects of dry-wet cycles on three-dimensional pore structure and permeability characteristics of granite residual soil using X-ray micro computed tomography

Due to seasonal climate alterations,the microstructure and permeability of granite residual soil are easily affected by multiple dry-wet cycles.The X-ray micro computed tomography(micro-CT)acted as a nondestructive tool for characterizing the microstructure of soil samples exposed to a range of damage levels induced by dry-wet cycles.Subsequently,the variations of pore distribution and permeability due to drywet cycling effects were revealed based on three-dimensional(3D)pore distribution analysis and seepage simulations.According to the results,granite residual soils could be separated into four different components,namely,pores,clay,quartz,and hematite,from micro-CT images.The reconstructed 3D pore models dynamically demonstrated the expanding and connecting patterns of pore structures during drywet cycles.The values of porosity and connectivity are positively correlated with the number of dry-wet cycles,which were expressed by exponential and linear functions,respectively.The pore volume probability distribution curves of granite residual soil coincide with the χ^(2)distribution curve,which verifies the effectiveness of the assumption of χ^(2)distribution probability.The pore volume distribution curves suggest that the pores in soils were divided into four types based on their volumes,i.e.micropores,mesopores,macropores,and cracks.From a quantitative and visual perspective,considerable small pores are gradually transformed into cracks with a large volume and a high connectivity.Under the action of dry-wet cycles,the number of seepage flow streamlines which contribute to water permeation in seepage simulation increases distinctly,as well as the permeability and hydraulic conductivity.The calculated hydraulic conductivity is comparable with measured ones with an acceptable error margin in general,verifying the accuracy of seepage simulations based on micro-CT results.

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study wave- induced pore water pressure in high permeability seabed. In the model, the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes （RANS） equations with a k-s closure, and Forchheimer equations are adopted for flow within the porous media. By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow, a highly efficient coupling between the two flows is implemented. The numerical tests are conducted to study the effects of seabed thickness, porosity, particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response. The results indicate that, as compared with regular wave-induced, solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters. The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed. The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.

Determination of full-scale pore size distribution of Gaomiaozi bentonite and its permeability prediction

Gaomiaozi(GMZ)bentonite is a potential buffer/backfill material for a deep geological disposal of highlevel radioactive waste.It has a wide pore size distribution(PSD)with sizes ranging from several nanometers to more than one hundred microns.Thus,properly characterizing the pore structures of GMZ bentonite is a challenging issue.In this study,pressure-controlled porosimetry(PCP),ratecontrolled porosimetry(RCP),and scanning electron microscopy(SEM)were used to investigate the PSD of GMZ bentonite,The results indicate that each method has its limitation,and a combined use of PCP and RCP is suitable to obtain the full-scale PSD of GMZ bentonite.Moreover,we also compared the full-scale PSD with nuclear magnetic resonance(NMR)result.It is found that there is no significant difference in the range of PSD characterization between NMR and mercury intrusion method(PCP and RCP).However,in a ce rtain range,the detection accuracy of NMR is higher than that of mercury injection method.Finally,permeability prediction based on PCP and SEM data was conducted,and both of the two methods were found to be able to predict the permeability.The combined method is effective to obtain the full-scale PSD of GMZ bentonite,which is the key to estimation of the sealing ability of bentonite buffer.

Anisotropic Rock Poroelasticity Evolution in Ultra-low Permeability Sandstones under Pore Pressure,Confining Pressure,and Temperature:Experiments with Biot's Coefficient

This study aimed to show anisotropic poroelasticity evolution in ultra-low permeability reservoirs under pore pressure,confining pressure,and temperature.Several groups of experiments examining Biot's coefficient under different conditions were carried out.Results showed that Biot's coefficient decreased with increased pore pressure,and the variation trend is linear,but the decreasing rate is variable between materials.Biot's coefficient increased with increased confining pressure;the variation trend is linear,but the increasing rate varies by material as well.Generally,Biot's coefficient remains stable with increased temperature.Lithology,clay mineral content,particle arrangement,and pore arrangement showed impacts on Biot's coefficient.For strong hydrophilic clay minerals,expansion in water could result in a strong surface adsorption reaction,which could result in an increased fluid bulk modulus and higher Biot's coefficient.For skeleton minerals with strong lipophilicity,such as quartz and feldspar,increased oil saturation will also result in an adsorption reaction,leading to increased fluid bulk modulus and a higher Biot's coefficient.The study's conclusions provide evidence of poroelasticity evolution of ultra-low permeability and help the enhancing oil recovery(EOR)process.

Pore-scale modeling of pore structure properties and wettability effect on permeability of low-rank coal

Permeability is a key parameter for coalbed methane development.Although the absolute permeability of coal has been extensively studied,wettability and pore structure properties continue to challenge the microscopic description of water-gas flow in coal.For this purpose,we reconstructed the microstructures of low-rank coal using micro-computed tomography(micro-CT)images.Pore geometry and pore-throat parameters are introduced to establish a relationship with absolute permeability.A dual-porosity pore network model is developed to study water-gas displacement under different wetting and pore structure properties.Results show that absolute permeability is significantly affected by pore geometry and can be described using a binary quadratic function of porosity and fractal dimension.Water-gas relative permeability varies significantly and the residual gas saturation is lower;the crossover saturation first decreased and then increased with increasing porosity under hydrophobic conditions.While the water relative permeability is lower and a certain amount of gas is trapped in complex pore-throat networks;the crossover saturation is higher under hydrophilic conditions.Models with large percolating porosity and well-developed pore networks have high displacement efficiency due to low capillary resistance and avoidance of trapping.This work provides a systematic description of absolute permeability and water-gas relative permeability in coal microstructure for enhanced gas recovery.

Pore network modeling of water block in low permeability reservoirs

A pore network model was used in this paper to investigate the factors, in particular, throat radius, wettability and initial water saturation, causing water block in low permeability reservoirs. A new term - ＇relative permeability number＇ （RPN） was firstly defined, and then used to describe the degree of water block. Imbibition process simulations show that the RPN drops in accordance with the extension of the averaged pore throat radius from 0.05 to 1.5 μm, and yet once beyond that point of 1.5 μm, the RPN reaches a higher value, indicating the existence of a critical pore throat radius where water block is the maximum. When the wettability of the samples changes from water-wet to weakly water-wet, weakly gas-wet, or gas（oil）-wet, the gas RPN increases consistently, but this consistency is disturbed by the RPN dropping for weakly water-wet samples for water saturations less than 0.4, which means weakly waterwet media are more easily water blocked than water-wet systems. In the situation where the initial water saturation exceeds 0.05, water block escalates along with an increase in initial water saturation.

Pore-scale simulation of gas-water flow in low permeability gas reservoirs

A novel method was developed to establish a realistic three dimensional(3D) network model representing pore space in low permeability sandstone.Digital core of rock sample was established by the combination of micro-CT scanning and image processing,then 3D pore-throat network model was extracted from the digital core through analyzing pore space topology,calculating pore-throat parameters and simplifying the shapes of pores and throats.The good agreements between predicted and measured porosity and absolute permeability verified the validity of this new network model.Gas-water flow mechanism was studied by using pore-scale simulations,and the influence of pore structure parameters,including coordination number,aspect ratio and shape factor,on gas-water flow,was investigated.The present simulation results show that with the increment of coordination number,gas flow ability in network improves and the effect of invading water on blocking gas flow weakens.The smaller the aspect ratio is,the stronger the anisotropy of the network is,resulting in the increase of seepage resistance.It is found that the shape factor mainly affects the end points in relative permeability curves,and for a highly irregular pore or throat with a small shape factor,the irreducible water saturation(Swi) and residual gas saturation(Sgr) are relatively high.

Pore structure differences of the extra-low permeability sandstone reservoirs and the causes of low resistivity oil layers: A case study of Block Yanwumao in the middle of Ordos Basin, NW China

The influence of pore structure difference on rock electrical characteristics of reservoir and oil reservoir was analyzed taking Triassic Chang 6 reservoir in Block Yanwumao in the middle of Ordos Basin as an example. The relationship between the pore structure difference and the low resistivity oil layer was revealed and demonstrated through core observation, lab experiments, geological research, well log interpretation and trial production etc. The results show that there were two kinds of oil layers in Chang 6 oil layer set, normal oil layer and low resistivity oil layer in the region, corresponding to two types of pore structures, pore type mono-medium and micro-fracture-pore type double-medium; the development of micro-fracture changed greatly the micro-pore structure of the reservoir, and the pore structure difference had an important influence on the rock electrical characteristics of the extra-low permeability sandstone reservoir and oil reservoir; the normal oil layers had obvious characteristics of pore-type mono-medium, and were concentrated in Chang 61, Chang 6232 and Chang 62; the low resistivity oil layers had obvious characteristics of micro-fracture-pore type double-medium, which were mainly distributed in Chang 612 and Chang 63. The mud filtrate penetrated deep into the oil layers along the micro-cracks, leading to sharp reduction of resistivity, and thus low resistivity of the oil layer; the low resistivity oil layers had better storage capacity and higher productivity than the normal oil layers.

A new permeability calculation method using nuclear magnetic resonance logging based on pore sizes: A case study of bioclastic limestone reservoirs in the A oilfield of the Mid-East

It is difficult to accurately obtain the permeability of complex lithologic reservoirs using conventional methods because they have diverse pore structures and complex seepage mechanisms.Based on in-depth analysis of the limitation of classical nuclear magnetic resonance(NMR)permeability calculation models,and the understanding that the pore structure and porosity are the main controlling factors of permeability,this study provides a new permeability calculation method involving classifying pore sizes by using NMR T_2 spectrum first and then calculating permeability of different sizes of pores.Based on this idea,taking the bioclastic limestone reservoir in the A oilfield of Mid-East as an example,the classification criterion of four kinds of pore sizes:coarse,medium,fine and micro throat,was established and transformed into NMR T_2 standard based on shapes and turning points of mercury intrusion capillary pressure curves.Then the proportions of the four kinds of pore sizes were obtained precisely based on the NMR logging data.A new NMR permeability calculation model of multicomponent pores combinations was established based on the contributions of pores in different sizes.The new method has been used in different blocks.The results show that the new method is more accurate than the traditional ones.

Experimental analysis of the pore structure, relative permeability, and water flooding characteristics of the Yan'an Formation sandstone, southwestern Ordos Basin

The oil and gas potential of the Yan'an Formation in the Ordos Basin has yet to be fully tapped. In this study, the pore structure, mobile fluid saturation, and water flooding micro-mechanism of the Yan'an Formation sandstone are systematically studied through the application of a series of rock physics and fluid experiments. The results show that there is a good positive correlation between porosity and permeability, and the reservoirs are divided into types Ⅰ, Ⅱ, and Ⅲ. Mercury injection tests show that the average pore throat radius of the oil-bearing reservoir ranges from 1 to 7 μm. The displacement pressure of the Yan'an Formation is also relatively low, and it decreases from 0.1 MPa to 0.01 MPa as the rock porosity increases from 11% to 18%. NMR tests show that small (diameter <0.5 μm) and medium pores (diameter ranging from 0.5 to 2.5 μm) are predominant in the reservoir. Different types of reservoirs have different characteristics of relative permeability curve. In addition, when the average oil recovery rate is less than 1 ml/min, the oil displacement efficiency increases faster. However, when the average oil recovery rate is between 1–3.5 ml/min, the oil displacement efficiency is maintained at around 27%–30%. Physical properties of the reservoir, pore-throat structure, experimental pressure difference, and pore volume injected — all have significant effects on oil displacement efficiency. For Type Ⅰ and Type Ⅱ reservoirs, the increase of the pore volume injected has a significant effect on oil displacement efficiency. However, for Type Ⅲ reservoirs, the change of pore volume injected has insignificant effect on oil displacement efficiency. This study provides a reference for the formulation of estimated ultimate recovery (EUR) measures for similar sandstone reservoirs.

Pore-Throat Combination Types and Gas-Water Relative Permeability Responses of Tight Gas Sandstone Reservoirs in the Zizhou Area of East Ordos Basin, China

With the aim of better understanding the tight gas reservoirs in the Zizhou area of east Ordos Basin,a total of 222 samples were collected from 50 wells for a series of experiments.In this study,three pore-throat combination types in sandstones were revealed and confirmed to play a controlling role in the distribution of throat size and the characteristics of gas-water relative permeability.The type-I sandstones are dominated by intercrystalline micropores connected by cluster throats,of which the distribution curves of throat size are narrow and have a strong single peak(peak ratio>30%).The pores in the type-II sandstones dominantly consist of secondary dissolution pores and intercrystalline micropores,and throats mainly occur as slice-shaped throats along cleavages between rigid grain margins and cluster throats in clay cement.The distribution curves of throat size for the type-II sandstones show a bimodal distribution with a substantial low-value region between the peaks(peak ratio<15%).Primary intergranular pores and secondary intergranular pores are mainly found in type-III samples,which are connected by various throats.The throat size distribution curves of type-III sandstones show a nearly normal distribution with low kurtosis(peak ratio<10%),and the micro-scale throat radii(>0.5μm)constitute a large proportion.From type-I to type-III sandstones,the irreducible water saturation(Swo)decreased;furthermore,the slope of the curves of Krw/Krg in two-phase saturation zone decreased and the two-phase saturation zone increased,indicating that the gas relative flow ability increased.Variations of the permeability exist in sandstones with different porethroat combination types,which indicate the type-III sandstones are better reservoirs,followed by type-II sandstones and type-I sandstones.As an important factor affecting the reservoir quality,the pore-throat combination type in sandstones is the cumulative expression of lithology and diagenetic modifications with strong heterogeneity.

Influence mechanism of pore-scale anisotropy and pore distribution heterogeneity on permeability of porous media

Based on micro-CT scanning experiments, three-dimensional digital cores of tight sandstones were established to quantitatively evaluate pore-scale anisotropy and pore-distribution heterogeneity. The quartet structure generation set method was used to generate three-dimensional anisotropic, heterogeneous porous media models. A multi-relaxation-time lattice Boltzmann model was applied to analyze relationships of permeability with pore-scale anisotropy and pore distribution heterogeneity, and the microscopic influence mechanism was also investigated. The tight sandstones are of complex pore morphology, strong anisotropy and pore distribution heterogeneity, while anisotropy factor has obvious directivity. The obvious anisotropy influences the orientation of long axis of pores and fluid flow path, making tortuosity smaller and flowing energy loss less in the direction with the greater anisotropy factor. The strong correlation of tortuosity and anisotropy is the inherent reason of anisotropy acting on permeability. The influence of pore distribution heterogeneity on permeability is the combined effects of specific surface area and tortuosity, while the product of specific surface area and tortuosity shows significantly negative correlation with heterogeneity. The stronger the pore distribution heterogeneity, the smaller the product and the greater the permeability. In addition, the permeability and tortuosity of complex porous media satisfy a power relation with a high fitting precision, which can be applied for approximate estimation of core permeability.