Reservoir heterogeneity analysis using multi-directional textural attributes from deep learning-based enhanced acoustic impedance inversion:A study from Poseidon,NW shelf Australia

Reservoir heterogeneities play a crucial role in governing reservoir performance and management.Traditionally,detailed and inter-well heterogeneity analyses are commonly performed by mapping seismic facies change in the seismic data,which is a time-intensive task.Many researchers have utilized a robust Grey-level co-occurrence matrix(GLCM)-based texture attributes to map reservoir heterogeneity.However,these attributes take seismic data as input and might not be sensitive to lateral lithology variation.To incorporate the lithology information,we have developed an innovative impedance-based texture approach using GLCM workflow by integrating 3D acoustic impedance volume(a rock propertybased attribute)obtained from a deep convolution network-based impedance inversion.Our proposed workflow is anticipated to be more sensitive toward mapping lateral changes than the conventional amplitude-based texture approach,wherein seismic data is used as input.To evaluate the improvement,we applied the proposed workflow to the full-stack 3D seismic data from the Poseidon field,NW-shelf,Australia.This study demonstrates that a better demarcation of reservoir gas sands with improved lateral continuity is achievable with the presented approach compared to the conventional approach.In addition,we assess the implication of multi-stage faulting on facies distribution for effective reservoir characterization.This study also suggests a well-bounded potential reservoir facies distribution along the parallel fault lines.Thus,the proposed approach provides an efficient strategy by integrating the impedance information with texture attributes to improve the inference on reservoir heterogeneity,which can serve as a promising tool for identifying potential reservoir zones for both production benefits and fluid storage.

A geochemical perspective on the genesis of Cenozoic basic volcanism in northeastern Turkey:an overview of metasomatism and heterogeneity of the sub-continental lithospheric mantle in a post-collisional setting

The post-collisional Cenozoic basic volcanic rocks in NE Turkey show temporal variations in whole-rock lithophile element and highly siderophile element(HSE)systematics that are mainly associated with the nature of sub-continental lithospheric mantle(SCLM)sources and parental melt generation.So far,the traditional whole-rock lithophile geochemical data of these basic volcanic rocks have provided important constraints on the nature of SCLM sources.Integrated lithophile element and HSE geochemical data of these basic volcanic rocks also reveal the heterogeneity of the SCLM source,which is principally related to variable metasomatism resulting from previous subduction(s)and post-collisional mantle-crust interactions in an extensional setting.Lithophile element geochemical features suggest that the parental magmas have derived from metasomatized spinel-to garnet-bearing SCLM sources for Eocene and Miocene basic volcanic rocks with subduction signatures whereas originated from spinel-to garnet-bearing SCLM sources for Mio-Pliocene and Plio-Quaternary basaltic volcanic rocks without the subduction signature.Lithophile element and HSE geo-chemistry also reveal that Eocene and Miocene basic vol-canic rocks were affected by more pronounced crustal contamination than the basaltic volcanic rocks of Mio-Pliocene and Quaternary.Furthermore,the integrated lithophile element and HSE compositions of these basic volcanic rocks,together with the regional asymmetric lithospheric delamination model,reveal that the compositional variation(especially due to metasomatism)was significant temporally in the heterogeneity of the SCLM sources from which parental magmas formed during the Cenozoic era.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.

Influence of Transportation Accessibility on Urban-rural Income Disparity and Its Spatial Heterogeneity

Transportation accessibility has been treated as an important means of reducing the urban-rural income disparity.However,only a few studies have examined the effects of different types of transportation accessibility on urban-rural income disparity and their spatial heterogeneity.Based on data from 285 prefecture-level(and above)Chinese cities in 2000,2005,2010,2015,and 2020,this study uses spatial econometric models to examine how highway accessibility and railway accessibility influence the urban-rural income disparity and to identify their spatial heterogeneity.The result reveals that highway accessibility and railway accessibility have‘coreperiphery’ring-like circle structures.The urban-rural income disparity exhibits strong spatial clustering effects.Both highway accessibility and railway accessibility are negatively associated with urban-rural income disparity,and the former having a greater effect size.Moreover,there is a substitution effect between highway accessibility and railway accessibility in the whole sample.Furthermore,these associations differ in geographic regions.In the central region,highway accessibility is more important in reducing the urban-rural income disparity,but its effect is weakened with the increase of railway accessibility.In the western region,railway accessibility has a larger effect on narrowing the urban-rural income disparity,and this effect is strengthened by the increase of highway accessibility.We conclude that improving transportation accessibility is conducive to reducing the urban-rural income disparity but its effect is spatial heterogenetic.Highways and railways should be developed in a coordinated manner to promote an integrated transport network system.

Island biogeography theory and the habitat heterogeneity jointly explain global patterns of Rhododendron diversity

Mountain biodiversity is of great importance to biogeography and ecology.However,it is unclear what ecological and evolutionary processes best explain the generation and maintenance of its high levels of species diversity.In this study,we determined which of six common hypotheses(e.g.,climate hypotheses,habitat heterogeneity hypothesis and island biogeography theory)best explain global patterns of species diversity in Rhododendron.We found that Rhododendron diversity patterns were most strongly explained by proxies of island biogeography theory(i.e.,mountain area)and habitat heterogeneity(i.e.,elevation range).When we examined other relationships important to island biogeography theory,we found that the planimetric area and the volume of mountains were positively correlated with the Rhododendron diversity,whereas the‘mountains-to-mainland’distance was negatively correlated with Rhododendron diversity and shared species.Our findings demonstrate that Rhododendron diversity can be explained by island biogeography theory and habitat heterogeneity,and mountains can be regarded as islands which supported island biogeography theory.

Impact of Reservoir Heterogeneity on Bitumen Content in the Mackay River Oil Sands,Athabasca(Canada)

The Lower Cretaceous Manville Group of Upper Mc Murray Formation is one of the main bitumen reservoirs in Athabasca.In this study,the relationship between reservoirs heterogeneity and bitumen geochemical characteristics were analyzed through core and microscopic observation,lab analysis,petrophysics and logging data.Based on the sedimentology framework,the formation environment of high-quality oil sand reservoirs and their significance for development were discussed.The results indicate that four types lithofacies were recognized in the Upper Mc Murray Formation based on their depositional characteristics.Each lithofacies reservoirs has unique physical properties,and is subject to varying degrees of degradation,resulting in diversity of bitumen content and geochemical composition.The tidal bar(TB)or tidal channel(TC)facies reservoir have excellent physical properties,which are evaluated as gas or water intervals due to strong degradation.The reservoir of sand bar(SB)facies was evaluated as oil intervals,due to its poor physical properties and weak degradation.The reservoir of mixed flat(MF)facies is composed of sand intercalated with laminated shale,which is evaluated as poor oil intervals due to its poor connectivity.The shale content in oil sand reservoir is very important for the reservoir physical properties and bitumen degradation degree.In the context of regional biodegradation,oil sand reservoirs with good physical properties will suffer from strong degradation,while oil sand reservoirs with relatively poor physical properties are more conducive to the bitumen preservation.

Spatial Heterogeneity and Risk Factors of Dental Caries in 12-Year-Old Children in Shanxi Province,China

Objective This study aimed to explore the spatial heterogeneity and risk factors for dental caries in 12-year-old children in Shanxi province,China.Methods The data encompassed 3,721 participants from the two most recent oral health surveys conducted across 16 districts in Shanxi Province in 2015 and 2018.Eighteen specific variables were analyzed to examine the interplay between socioeconomic factors,medical resources and environmental conditions.The Geo-detector model was employed to assess the impacts and interactions of these ecological factors.Results Socioeconomic factors(Q=0.30,P<0.05)exhibited a more substantial impact compared to environmental(Q=0.19,P<0.05)and medical resource factors(Q=0.25,P<0.05).Notably,the urban population percentage(UPP)demonstrated the most significant explanatory power for the spatial heterogeneity in caries prevalence,as denoted by its highest q-value(q=0.51,P<0.05).Additionally,the spatial distribution’s heterogeneity of caries was significantly affected by SO2 concentration(q=0.39,P<0.05)and water fluoride levels(q=0.27,P<0.05)among environmental factors.Conclusion The prevalence of caries exhibited spatial heterogeneity,escalating from North to South in Shanxi Province,China,influenced by socioeconomic factors,medical resources,and environmental conditions to varying extents.

Spatial Heterogeneity Modeling Using Machine Learning Based on a Hybrid of Random Forest and Convolutional Neural Network (CNN)

Spatial heterogeneity refers to the variation or differences in characteristics or features across different locations or areas in space. Spatial data refers to information that explicitly or indirectly belongs to a particular geographic region or location, also known as geo-spatial data or geographic information. Focusing on spatial heterogeneity, we present a hybrid machine learning model combining two competitive algorithms: the Random Forest Regressor and CNN. The model is fine-tuned using cross validation for hyper-parameter adjustment and performance evaluation, ensuring robustness and generalization. Our approach integrates Global Moran’s I for examining global autocorrelation, and local Moran’s I for assessing local spatial autocorrelation in the residuals. To validate our approach, we implemented the hybrid model on a real-world dataset and compared its performance with that of the traditional machine learning models. Results indicate superior performance with an R-squared of 0.90, outperforming RF 0.84 and CNN 0.74. This study contributed to a detailed understanding of spatial variations in data considering the geographical information (Longitude & Latitude) present in the dataset. Our results, also assessed using the Root Mean Squared Error (RMSE), indicated that the hybrid yielded lower errors, showing a deviation of 53.65% from the RF model and 63.24% from the CNN model. Additionally, the global Moran’s I index was observed to be 0.10. This study underscores that the hybrid was able to predict correctly the house prices both in clusters and in dispersed areas.

Type 2 diabetes in children and adolescents:Exploring the disease heterogeneity and research gaps to optimum management

Over the past 20 years,the incidence and prevalence of type 2 diabetes mellitus(T2DM)in children and adolescents have increased,particularly in racial and ethnic minorities.Despite the rise in T2DM in children and adolescents,the pathophysiology and progression of disease in this population are not clearly understood.Youth-onset T2DM has a more adverse clinical course than is seen in those who develop T2DM in adulthood or those with T1DM.Furthermore,the available therapeutic options are more limited for children and adolescents with T2DM compared to adult patients,mostly due to the challenges of implementing clinical trials.A better understanding of the mechanisms underlying the development and aggressive disease phenotype of T2DM in youth is important to finding effective prevention and management strategies.This review highlights the key evidence about T2DM in children and adolescents and its current burden and challenges both in clinical care and research activities.

Effect of Modulus Heterogeneity on the Equilibrium Shape and Stress Field ofαPrecipitate in Ti-6Al-4V

For media with inclusions(e.g.,precipitates,voids,reinforcements,and others),the difference in lattice parameter and the elastic modulus between the matrix and inclusions cause stress concentration at the interfaces.These stress fields depend on the inclusions’size,shape,and distribution and will respond instantly to the evolving microstructure.This study develops a phase-field model concerningmodulus heterogeneity.The effect of modulus heterogeneity on the growth process and equilibrium state of theαplate in Ti-6Al-4V during precipitation is evaluated.Theαprecipitate exhibits strong anisotropy in shape upon cooling due to the interplay of the elastic strain and interfacial energy.The calculated orientation of the habit plane using the homogeneous modulus ofαphase shows the smallest deviation fromthat of the habit plane observed in the experiment,compared to the case where the homogeneous modulus ofβphase is adopted.In addition,the equilibrium volume ofαphase within the systemusing homogeneousβmodulus exhibits the largest dependency on the applied stresses.The stress fields across theα/βinterface are further calculated under the assumption of modulus heterogeneity and compared to those using homogeneous modulus of eitherαorβphase.This study provides an essential theoretical basis for developing mechanics models concerning systems with heterogeneous structures.

Quantifying the impact of earthquakes and geological factors on spatial heterogeneity of debris-flow prone areas:A case study in the Hengduan Mountains

Understanding the spatial heterogeneity of debris-flow-prone areas holds significant implications for regional risk management, particularly in seismically active regions with geological faults. Despite the significance of this knowledge, a comprehensive quantification of the influence of regional topographical and geological factors on the spatial heterogeneity of debris-flow-prone areas has been lacking. This study selected the Hengduan Mountains, an earthquake-prone region characterized by diverse surface conditions and complex landforms, as a representative study area. An improved units zoning and objective factors identification methodology was employed in earthquake and fault analysis to assess the impact of seismic activity and geological factors on spatial heterogeneity of debrisflow prone areas. Results showed that the application of GIS technology with hydrodynamic intensity and geographical units analysis can effectively analyze debris-flow prone areas. Meanwhile, earthquake and fault zones obviously increase the density of debrisflow prone catchments and make them unevenly distributed. The number of debris-flow prone areas shows a nonlinear variation with the gradual increase of geomorphic factor value. Specifically, the area with 1000 m-2500 m elevation difference, 25°-30° average slope, and 0.13-0.15 land use index is the most favorable conditions for debris-flow occurrence;The average annual rainfall from 600 to 1150 mm and landslides gradient from 16° to 35° are the main causal factors to trigger debris flow. Our study sheds light on the quantification of spatial heterogeneity in debris flow-prone areas in earthquake-prone regions, which can offer crucial support for post-debris flow risk management strategies.

Age-specific heterogeneity of genetic susceptibility to cardiovascular disease might have opposite outcomes depending on the presence of prediabetes

Examining age-specific heterogeneity of susceptibility to cardiovascular disease is also essential in individuals without prediabetes to determine its relative size and direction compared to those with prediabetes.Of particular interest,age-specific heterogeneity in genetic susceptibility may exhibit opposite directions depending on the presence or absence of prediabetes.

Prediction of Porous Media Fluid Flow with Spatial Heterogeneity Using Criss-Cross Physics-Informed Convolutional Neural Networks

Recent advances in deep neural networks have shed new light on physics,engineering,and scientific computing.Reconciling the data-centered viewpoint with physical simulation is one of the research hotspots.The physicsinformedneural network(PINN)is currently the most general framework,which is more popular due to theconvenience of constructing NNs and excellent generalization ability.The automatic differentiation(AD)-basedPINN model is suitable for the homogeneous scientific problem;however,it is unclear how AD can enforce fluxcontinuity across boundaries between cells of different properties where spatial heterogeneity is represented bygrid cells with different physical properties.In this work,we propose a criss-cross physics-informed convolutionalneural network(CC-PINN)learning architecture,aiming to learn the solution of parametric PDEs with spatialheterogeneity of physical properties.To achieve the seamless enforcement of flux continuity and integration ofphysicalmeaning into CNN,a predefined 2D convolutional layer is proposed to accurately express transmissibilitybetween adjacent cells.The efficacy of the proposedmethodwas evaluated through predictions of several petroleumreservoir problems with spatial heterogeneity and compared against state-of-the-art(PINN)through numericalanalysis as a benchmark,which demonstrated the superiority of the proposed method over the PINN.

Optimal Insurance with Background Risk under the Ambiguity and Belief Heterogeneity Structure

In this paper, we discuss the optimal insurance in the presence of background risk while the insured is ambiguity averse and there exists belief heterogeneity between the insured and the insurer. We give the optimal insurance contract when maxing the insured’s expected utility of his/her remaining wealth under the smooth ambiguity model and the heterogeneous belief form satisfying the MHR condition. We calculate the insurance premium by using generalized Wang’s premium and also introduce a series of stochastic orders proposed by [1] to describe the relationships among the insurable risk, background risk and ambiguity parameter. We obtain the deductible insurance is the optimal insurance while they meet specific dependence structures.

Insight into the spoilage heterogeneity of meat-borne bacteria isolates with high-producing collagenase

Chilled chicken is inevitably contaminated by microorganisms during slaughtering and processing,resulting in spoilage.Cutting parts of chilled chicken,especially wings,feet,and other skin-on products,are abundant in collagen,which may be the primary target for degradation by spoilage microorganisms.In this work,a total of 17 isolates of spoilage bacteria that could secrete both collagenase and lipase were determined by raw-chicken juice agar(RJA)method,and the results showed that 7 strains of Serratia,Aeromonas,and Pseudomonas could significantly decompose the collagen ingredients.The gelatin zymography showed that Serratia liquefaciens(F5)and(G7)had apparent degradation bands around 50 kDa,and Aeromonas veronii(G8)and Aeromonas salmonicida(H8)had a band around.65 and 95 kDa,respectively.The lipase and collagenase activities were detected isolate-by-isolate,with F5 showing the highest collagenase activity.For spoilage ability on meat in situ,F5 performed strongest in spoilage ability,indicated by the total viable counts,total volatile basic nitrogen content,sensory scores,lipase,and collagenase activity.This study provides a theoretical basis for spoilage heterogeneity of strains with high-producing collagenase in meat.

Spatial-temporal Heterogeneity of Green Development Efficiency and Its Influencing Factors in Growing Metropolitan Area: A Case Study for the Xuzhou Metropolitan Area

This study analyzed the spatial-temporal heterogeneity of green development efficiency and its influencing factors in the growing Xuzhou Metropolitan Area for the period 2000–2015.The slacks-based measure(SBM)model,spatial autocorrelation,and the geographically weighted regression(GWR)model were used to conduct the analysis.The conclusions were as follows:first,the overall efficiency of green development of the Xuzhou Metropolitan Area decreased,the regional differences and spatial agglomeration shrunk and differences within the region were the main contributors to the regional differences of green development efficiency.Second,the counties with high-efficiency green development were distributed along the coast,and along the routes of the Beijing-Shanghai and the Eastern Longhai railways.A developing axis of the high-efficiency counties was the main feature of the spatial pattern for green development efficiency.Third,regarding spatial correlation and green development efficiency,the High-High type counties in the Xuzhou Metropolitan Area formed a centralized distribution corridor along the inter-provincial border areas of Henan and Jiangsu,whereas the Low-Low type counties were concentrated in the external,marginal parts of the metropolitan area.Fourth,the major factors(ranked in decreasing order of impact)influencing green development efficiency were innovation,government regulations,the economic development level,energy consumption,and industrial structure.These factors exerted their influence to varying extents;the influence of the same factor had different effects in different regions and obvious spatial differences were observed for the different regions.

Inversion of walkaway VSP data in the presence of lateral velocity heterogeneity

Multi-azimuth walkaway vertical seismic profiling is an established technique for the estimation of in situ slowness surfaces and inferring anisotropy parameters.Normally,this technique requires the assumption of lateral homogeneity,which makes the horizontal slowness components at depths of downhole receivers equal to those measured at the surface.Any violations of this assumption,such as lateral heterogeneity or nonzero dip of intermediate interfaces,lead to distortions in reconstructed slowness surfaces and,consequently,to errors in estimated anisotropic parameters.In this work,we relax the assumption of lateral homogeneity and discuss how to correct vertical seismic profile data for weak lateral heterogeneity.We describe a procedure of downward continuation of recorded traveltimes that accounts for the presence of both vertical inhomogeneity and weak lateral heterogeneity,which produces correct slowness surfaces at depths of downhole receivers,noticing that sufficiently dense receiver coverage along a borehole is required to separate influences of vertical and lateral heterogeneity on measured traveltimes and obtain accurate estimates of the slowness surfaces.Once the slowness surfaces are found and a desired type of anisotropic model to be inverted is selected,the corresponding anisotropic parameters,providing the best fit to the estimated slownesses,can be obtained.We invert the slowness surfaces of P-waves for parameters of the simplest anisotropic model describing dipping fractures(transversely isotropic medium with a tilted symmetry axis).Five parameters of this model,namely,the P-wave velocity V0 in the direction of the symmetry axis,Thomsen's anisotropic coefficients e and d,the tilt n,and the azimuth b of the symmetry axis,can be estimated in a stable manner when maximum source offset is greater than half of receiver depth.

Effects of microscopic pore structure heterogeneity on the distribution and morphology of remaining oil

Waterflooding experiments were performed using Micro-CT on four cores of different pore structures from Donghe sandstone reservoirs in the Tarim Basin. The water, oil and grains were accurately separated by the advanced image processing technology, the pore network model was established, and parameters such as the number of throats and the throat size distribution were calculated to characterize the microscopic heterogeneity of pore structure, the flow of oil phase during displacement, and the morphology and distribution of remaining oil after displacement. The cores with the same macroscopic porosity-permeability have great differences in microscopic heterogeneity of pore structure. Both macro porosity-permeability and micro heterogeneity of pore structure have an influence on the migration of oil phase and the morphology and distribution of remaining oil. When the heterogeneity is strong, the water phase will preferentially flow through the dominant paths and the remaining oil clusters will be formed in the small pores. The more the number of oil clusters(droplets) formed during displacement process, the smaller the average volume of cluster is, and the remaining oil is dominated by the cluster continuous phase with high saturation. The weaker the heterogeneity, the higher the pore sweep efficiency is, and the remaining oil clusters are mainly trapped in the form of non-continuous phase. The distribution and morphology of micro remaining oil are related to the absolute permeability, capillary number and micro-heterogeneity. So, the identification plate of microscopic residual oil continuity distribution established on this basis can describe the relationship between these three factors and distribution of remaining oil and identify the continuity of the remaining oil distribution accurately.

Pore Structure Heterogeneity of the Xiamaling Formation Shale Gas Reservoir in the Yanshan Area of China: Evaluation of Geological Controlling Factors

Micro-heterogeneity is an integral parameter of the pore structure of shale gas reservoir and it forms an essential basis for setting and adjusting development parameters.In this study,scanning electron microscopy,high-pressure mercury intrusion and low-temperature nitrogen adsorption experiments were used to qualitatively and quantitatively characterize the pore structure of black shale from the third member of the Xiamaling Formation in the Yanshan area.The pore heterogeneity was studied using fractal theory,and the controlling factors of pore development and heterogeneity were evaluated in combination with geochemical parameters,mineral composition,and geological evolution history.The results show that the pore structure of the reservoir was intricate and complicated.Moreover,various types of micro-nano scale pores such as dissolution pores,intergranular pores,interlayer pores,and micro-cracks are well developed in member 3 of the Xiamaling Formation.The average porosity was found to be 6.30%,and the mean value of the average pore size was 4.78 nm.Micropores and transition pores provided most of the storage space.Pore development was significantly affected by the region and was mainly related to the total organic carbon content,vitrinite reflectance and mineral composition.The fractal dimension,which characterizes the heterogeneity,is 2.66 on average,indicating that the pore structure is highly heterogeneous.Fractal dimension is positively correlated with maturity and clay mineral content,while it is negatively correlated with brittle mineral content and average pore size.These results indicate that pore heterogeneity is closely related to thermal history and material composition.Combined with the geological background of this area,it was found that the pore heterogeneity was mainly controlled by the Jurassic magmatism.The more intense the magma intrusion,the stronger the pore heterogeneity.The pore structure and its heterogeneity characteristics present today are a general reflection of the superimposed geological processes of sedimentary-diagenetic-late transformation.The influence of magmatic intrusion on the reservoir is the main geological factor that should be considered for detailed evaluation of the Xiamaling Formation shale gas reservoir in the Yanshan area.

HPV18 E6 and E7 Intratumour Heterogeneity in Esophageal Cancer

The development of esophageal cancer accompanied by the presence of human papillomavirus (HPV) DNA into the host genome. By evaluating the expression of this virus for tumor cell origin and also their cell grows and migrations, we examined esophageal cancer clonality in the context of intra-tumor heterogeneity. In this research, we have checked the expression of HPV18 E6 and E7 in different single cell clones by the manual cell picking method in the HPV positive esophageal cancer (EC109), EC109 cell line used as a negative control, and Hela cell line used as the positive control. Quantitative real-time PCR (QRT-PCR) was run to detect the expression levels of HPV E6 and E7, Cell Counting Kit-8 (CCK-8) assay was used to examine cell proliferation, invasion assays performed using Costar chambers and wounding assay to study cell migrations in vitro. We investigated the intra-tumor heterogeneity of HPV E6 and E7 in esophageal cancer and the evaluation of the growth and migrations at the clonal level, using 10 single cell clones. In particular clones, C7 & C10 displayed a highly variable expression in both HPV E6 and E7 and weak in four clones (C1, C3, C4, and C9) consequently, the cell invasion, proliferation, and migration increase with increasing the level of HPV expression and inverse. In conclusion, the resulting based on single cell cloning showed the relationship between HPV and cell growth and migration in esophageal cancer. Future study in HPV DNA integration needed to explore the mains specific integration site of HPV DNA in esophageal cancer and molecular monitoring of the HPV for future prevention researches and also effective therapeutic strategies.