With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can b...With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can be extended,and the prediction of the depth limit of oil accumulation(DLOA),are issues that have attracted significant attention in petroleum geology.Since it is difficult to characterize the evolution of the physical properties of the marine carbonate reservoir with burial depth,and the deepest drilling still cannot reach the DLOA.Hence,the DLOA cannot be predicted by directly establishing the relationship between the ratio of drilling to the dry layer and the depth.In this study,by establishing the relationships between the porosity and the depth and dry layer ratio of the carbonate reservoir,the relationships between the depth and dry layer ratio were obtained collectively.The depth corresponding to a dry layer ratio of 100%is the DLOA.Based on this,a quantitative prediction model for the DLOA was finally built.The results indicate that the porosity of the carbonate reservoir,Lower Ordovician in Tazhong area of Tarim Basin,tends to decrease with burial depth,and manifests as an overall low porosity reservoir in deep layer.The critical porosity of the DLOA was 1.8%,which is the critical geological condition corresponding to a 100%dry layer ratio encountered in the reservoir.The depth of the DLOA was 9,000 m.This study provides a new method for DLOA prediction that is beneficial for a deeper understanding of oil accumulation,and is of great importance for scientific guidance on deep oil drilling.展开更多
On 2019-03-04,the largest induced earthquake(ML4.18)occurred in the East Shale Basin,Alberta,and the underlying physical mechanisms have not been fully understood.This paper proposes a synthetical geoengineering metho...On 2019-03-04,the largest induced earthquake(ML4.18)occurred in the East Shale Basin,Alberta,and the underlying physical mechanisms have not been fully understood.This paper proposes a synthetical geoengineering methodology to comprehensively characterize this earthquake caused by hydraulic fracturing.Based on 3D structural,petrophysical,and geomechanical models,an unconventional fracture model is constructed by considering the stress shadow between adjacent hydraulic fractures and the interactions between hydraulic and natural fractures.Coupled poroelastic simulations are conducted to reveal the triggering mechanisms of induced seismicity.It is found that four vertical basement-rooted faults were identified via focal mechanisms analysis.The brittleness index(BI)along two horizontal wells has a high magnitude(BI>0.5),indicating the potential susceptibility of rock brittleness.Due to the presence of overpressure,pre-existing faults in the Duvernay Formation are highly susceptible to fault reactivation.The occurrence of the earthquake clusters has been attributed to the fracturing fluid injection during the west 38^(th)-39^(th) stage and east 38^(th) stage completions.Rock brittleness,formation overpressure,and large fracturing job size account for the nucleation of earthquake clusters,and unconventional natural-hydraulic fracture networks provide fluid flow pathways to cause fault reactivation.This workflow can be used to mitigate potential seismic risks in unconventional reservoirs in other fields.展开更多
A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distri...A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distribution and magnitude of hydraulic fracturing-induced seismicity. An integrated machine learning-based investigation is conducted to systematically evaluate multiple factors that contribute to induced seismicity. Feature importance indicates that a distance to fault, a distance to basement, minimum principal stress, cumulative fluid injection, initial formation pressure, and the number of fracturing stages are among significant model predictors. Our seismicity prediction map matches the observed spatial seismicity, and the prediction model successfully guides the fracturing job size of a new well to reduce seismicity risks. This study can apply to mitigating potential seismicity risks in other seismicity-frequent regions.展开更多
BACKGROUND Postoperative complications of triceps surae intramuscular hemangioma surgery with talipes equinovarus have rarely been described,and the evidence for treatment is limited.The purpose of this case study was...BACKGROUND Postoperative complications of triceps surae intramuscular hemangioma surgery with talipes equinovarus have rarely been described,and the evidence for treatment is limited.The purpose of this case study was to report the new application of the Ilizarov technique,which successfully treated talipes equinovarus in adults after triceps surae intramuscular hemangioma.CASE SUMMARY A 29-year-old woman treated with the Ilizarov technique for talipes equinovarus in the right leg after triceps surae intramuscular hemangioma surgery.The equinus deformity was roughly corrected after 2 years of follow-up,without significant secondary sequelae.CONCLUSION Talipes equinovarus caused by postoperative sequelae of intramuscular hemangioma was successfully corrected by the Ilizarov technique.The Ilizarov technique may be used for treating talipes equinovarus caused by various causes.展开更多
BACKGROUND Rotationplasty is often performed for malignant tumors,but type BIIIb rotationplasty is rarely reported,and there needs to be more evidence of the procedure and treatment.The purpose of this case study was ...BACKGROUND Rotationplasty is often performed for malignant tumors,but type BIIIb rotationplasty is rarely reported,and there needs to be more evidence of the procedure and treatment.The purpose of this case study was to report a new direction in the use of type BIIIb rotationplasty in treating patients with limb salvage and longterm non-healing infections.CASE SUMMARY Case 1:A 47-year-old man underwent radiotherapy for hemangioendothelioma in his left thigh,resulting in a femoral fracture.Despite the use of plates,intramedullary nailing,and external fixators,the femoral bone failed to unite due to infectious nonunion.Multiple operations were unable to control the infection,leaving the patient immobile.We performed a modified tibia-pelvic-constrained hip rotationplasty,utilizing a constrained prosthetic hip between the tibia and pelvis following a femur resection.Two years post-surgery,the patient was able to walk with the prosthetic device without any signs of recurring infection.The corresponding functional scores were 72 points for the Musculoskeletal Tumor Society(MSTS),53 for the Functional Mobility Assessment(FMA),93 for the Toronto Extremity Salvage Score(TESS),and 56 for the MOS 36-item short form health survey(SF-36).Case 2:A 59-year-old woman presented with liposarcoma in her left thigh.Magnetic resonance imaging revealed tumors in the medial,anterior,and posterior femur muscles,encircling the femoral vessels and nerves.Fortunately,there were no symptoms of sciatic dysfunction,and the tumor had not invaded the sciatic nucleus.After one year of follow-up,the patient expressed satisfaction with limb preservation post-type BIIIb rotationplasty.The corresponding functional scores were 63 points for the MSTS,47 for the FMA,88 for the TESS,and 52 for the SF-36.CONCLUSION Our study suggests that type BIIIb rotationplasty may be an alternative to amputation in patients with incurable infections.For malignant tumors of the lower extremities without invasion of the sciatic nerve,type BIIIb rotationplasty remains an excellent alternative to amputation.This surgical method may prevent amputation,improve functional outcomes,and facilitate biological reconstruction.展开更多
Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeabili...Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.展开更多
Hydrocarbon expulsion features and resource potential evaluation of source rocks are crucial for the petroleum exploration.High-maturity marine source rocks have not exhibited a hydrocarbon expulsion mode owing to the...Hydrocarbon expulsion features and resource potential evaluation of source rocks are crucial for the petroleum exploration.High-maturity marine source rocks have not exhibited a hydrocarbon expulsion mode owing to the lack of low-maturity source rocks in deep petroliferous basins.We considered the Ediacaran microbial dolomite in the Sichuan Basin,the largest high-maturity marine gas layer in China,to exhibit a method that quantitatively characterizes the hydrocarbon expulsion of high-maturity marine source rocks.The experiment of fluid inclusion,rock pyrolysis,and vitrinite reflectance(Ro)of 119 microbial dolomite core samples obtained from the Dengying Formation were performed.A hydrocarbon expulsion model of high-maturity source rock was established,and its resource potential was evaluated.The results showed that the Ediacaran microbial dolomite in the Sichuan Basin is a good source rock showing vast resource potential.The hydrocarbon expulsion threshold is determined to be vitrinite reflectance at 0.92%.The hydrocarbon expulsion intensities in the geologic history is high with maximum of 1.6×10^(7)t/km^(2).The Ediacaran microbial dolomite expelled approximately 1.008×10^(12)t of hydrocarbons,and the recoverable resource was 1.5×10^(12)m^(3).The region can be categorized into areasⅠ,Ⅱ,Ⅲ,andⅣ,in decreasing order of hydrocarbon expulsion intensity.Areas with a higher hydrocarbon expulsion intensity have a lower drilling risk and should be prioritized for exploration in the orderⅠ>Ⅱ>Ⅲ>Ⅳ.Two areas,northern and central parts of Ediacaran in the Sichuan Basin,were selected as prospects which had the drilling priority in the future gas exploration.The production data of 55 drilled wells verified the high reliability of this method.This model in this study does not require low-maturity samples and can be used for evaluating high-maturity marine source rocks,which has broad applicability in deep basins worldwide.展开更多
AIM: TO screen the differential expressed genes in colorectal cancer and polyp tissue samples. METHODS: Tissue specimens containing 16 cases of colorectal adenocarcinoma and colorectal polyp vs nor- mal mucosae were...AIM: TO screen the differential expressed genes in colorectal cancer and polyp tissue samples. METHODS: Tissue specimens containing 16 cases of colorectal adenocarcinoma and colorectal polyp vs nor- mal mucosae were collected and subjected to cDNA microarray and bioinformatical analyses. Quantitative reverse transcription-polymerase chain reaction (qRT- PCR) was used to confirm some of the cDNA microarray data.RESULTS: The experimental data showed that eight genes were differentially expressed, most of which were upregulated in adenomatous polyp lesions. Forty-six genes expressions were altered in colorectal cancers, of which 29 were upregulated and 17 downregulated, as compared to the normal mucosae. In addition, 18 genes were similarly altered in both adenomatous polyps and colorectal cancer, qRT-PCR analyses confirmed the cDNA microarray data for four of those 18 genes: MTA1, PDCD4, TSC1 and PDGFRA. CONCLUSION: These differentially expressed genes likely represent biomarkers for early detection of co- Iorectal cancer and may be potential therapeutic targets after confirmed by further studies.展开更多
In this paper we consider mixed finite element methods for second order elliptic problems. In the case of the lowest order Brezzi-Douglas-Marini elements (if d = 2) or Brezzi- Douglas-Duran-Fortin elements (if d = ...In this paper we consider mixed finite element methods for second order elliptic problems. In the case of the lowest order Brezzi-Douglas-Marini elements (if d = 2) or Brezzi- Douglas-Duran-Fortin elements (if d = 3) on rectangular parallelepipeds, we show that the mixed method system, by incorporating certain quadrature rules, can be written as a simple, cell-centered finite difference method. This leads to the solution of a sparse, positive semidefinite linear system for the scalar unknown. For a diagonal tensor coefficient, the sparsity pattern for the scalar unknown is a five point stencil if d = 2, and seven if d = 3. For a general tensor coefficient, it is a nine point stencil, and nineteen, respectively. Applications of the mixed method implementation as finite differences to nonisothermal multiphase, multicomponent flow in porous media are presented.展开更多
基金This work was supported by the Beijing Nova Program[Z211100002121136]Open Fund Project of State Key Laboratory of Lithospheric Evolution[SKL-K202103]+1 种基金Joint Funds of National Natural Science Foundation of China[U19B6003-02]the National Natural Science Foundation of China[42302149].We would like to thank Prof.Zhu Rixiang from the Institute of Geology and Geophysics,Chinese Academy of Sciences.
文摘With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can be extended,and the prediction of the depth limit of oil accumulation(DLOA),are issues that have attracted significant attention in petroleum geology.Since it is difficult to characterize the evolution of the physical properties of the marine carbonate reservoir with burial depth,and the deepest drilling still cannot reach the DLOA.Hence,the DLOA cannot be predicted by directly establishing the relationship between the ratio of drilling to the dry layer and the depth.In this study,by establishing the relationships between the porosity and the depth and dry layer ratio of the carbonate reservoir,the relationships between the depth and dry layer ratio were obtained collectively.The depth corresponding to a dry layer ratio of 100%is the DLOA.Based on this,a quantitative prediction model for the DLOA was finally built.The results indicate that the porosity of the carbonate reservoir,Lower Ordovician in Tazhong area of Tarim Basin,tends to decrease with burial depth,and manifests as an overall low porosity reservoir in deep layer.The critical porosity of the DLOA was 1.8%,which is the critical geological condition corresponding to a 100%dry layer ratio encountered in the reservoir.The depth of the DLOA was 9,000 m.This study provides a new method for DLOA prediction that is beneficial for a deeper understanding of oil accumulation,and is of great importance for scientific guidance on deep oil drilling.
基金This research was supported by Science Foundation of China University of Petroleum,Beijing(No.2462023BJRC001)National Natural Science Foundation of China Project(No.52204039).
文摘On 2019-03-04,the largest induced earthquake(ML4.18)occurred in the East Shale Basin,Alberta,and the underlying physical mechanisms have not been fully understood.This paper proposes a synthetical geoengineering methodology to comprehensively characterize this earthquake caused by hydraulic fracturing.Based on 3D structural,petrophysical,and geomechanical models,an unconventional fracture model is constructed by considering the stress shadow between adjacent hydraulic fractures and the interactions between hydraulic and natural fractures.Coupled poroelastic simulations are conducted to reveal the triggering mechanisms of induced seismicity.It is found that four vertical basement-rooted faults were identified via focal mechanisms analysis.The brittleness index(BI)along two horizontal wells has a high magnitude(BI>0.5),indicating the potential susceptibility of rock brittleness.Due to the presence of overpressure,pre-existing faults in the Duvernay Formation are highly susceptible to fault reactivation.The occurrence of the earthquake clusters has been attributed to the fracturing fluid injection during the west 38^(th)-39^(th) stage and east 38^(th) stage completions.Rock brittleness,formation overpressure,and large fracturing job size account for the nucleation of earthquake clusters,and unconventional natural-hydraulic fracture networks provide fluid flow pathways to cause fault reactivation.This workflow can be used to mitigate potential seismic risks in unconventional reservoirs in other fields.
基金This research has been made possible by contributions from the Natural Sciences and Engineering Research Council(NSERC)/Energi Simulation Industrial Research Chair in Reservoir Simulation and the Alberta Innovates(iCore)Chair in Reservoir ModelingThis research was supported by the Science Foundation of China University of Petroleum,Beijing(No.2462023BJRC001)the National Natural Science Foundation of China Joint Fund Key Support Project(No.U19B6003).
文摘A comprehensive dataset from 594 fracturing wells throughout the Duvernay Formation near Fox Creek, Alberta, is collected to quantify the influences of geological, geomechanical, and operational features on the distribution and magnitude of hydraulic fracturing-induced seismicity. An integrated machine learning-based investigation is conducted to systematically evaluate multiple factors that contribute to induced seismicity. Feature importance indicates that a distance to fault, a distance to basement, minimum principal stress, cumulative fluid injection, initial formation pressure, and the number of fracturing stages are among significant model predictors. Our seismicity prediction map matches the observed spatial seismicity, and the prediction model successfully guides the fracturing job size of a new well to reduce seismicity risks. This study can apply to mitigating potential seismicity risks in other seismicity-frequent regions.
基金Supported by Zhangzhou Natural Science Foundation Projects in 2019,No.ZZ2019J19。
文摘BACKGROUND Postoperative complications of triceps surae intramuscular hemangioma surgery with talipes equinovarus have rarely been described,and the evidence for treatment is limited.The purpose of this case study was to report the new application of the Ilizarov technique,which successfully treated talipes equinovarus in adults after triceps surae intramuscular hemangioma.CASE SUMMARY A 29-year-old woman treated with the Ilizarov technique for talipes equinovarus in the right leg after triceps surae intramuscular hemangioma surgery.The equinus deformity was roughly corrected after 2 years of follow-up,without significant secondary sequelae.CONCLUSION Talipes equinovarus caused by postoperative sequelae of intramuscular hemangioma was successfully corrected by the Ilizarov technique.The Ilizarov technique may be used for treating talipes equinovarus caused by various causes.
基金Supported by the Army Logistics Scientific Research Projects,No.CNJ16C013.
文摘BACKGROUND Rotationplasty is often performed for malignant tumors,but type BIIIb rotationplasty is rarely reported,and there needs to be more evidence of the procedure and treatment.The purpose of this case study was to report a new direction in the use of type BIIIb rotationplasty in treating patients with limb salvage and longterm non-healing infections.CASE SUMMARY Case 1:A 47-year-old man underwent radiotherapy for hemangioendothelioma in his left thigh,resulting in a femoral fracture.Despite the use of plates,intramedullary nailing,and external fixators,the femoral bone failed to unite due to infectious nonunion.Multiple operations were unable to control the infection,leaving the patient immobile.We performed a modified tibia-pelvic-constrained hip rotationplasty,utilizing a constrained prosthetic hip between the tibia and pelvis following a femur resection.Two years post-surgery,the patient was able to walk with the prosthetic device without any signs of recurring infection.The corresponding functional scores were 72 points for the Musculoskeletal Tumor Society(MSTS),53 for the Functional Mobility Assessment(FMA),93 for the Toronto Extremity Salvage Score(TESS),and 56 for the MOS 36-item short form health survey(SF-36).Case 2:A 59-year-old woman presented with liposarcoma in her left thigh.Magnetic resonance imaging revealed tumors in the medial,anterior,and posterior femur muscles,encircling the femoral vessels and nerves.Fortunately,there were no symptoms of sciatic dysfunction,and the tumor had not invaded the sciatic nucleus.After one year of follow-up,the patient expressed satisfaction with limb preservation post-type BIIIb rotationplasty.The corresponding functional scores were 63 points for the MSTS,47 for the FMA,88 for the TESS,and 52 for the SF-36.CONCLUSION Our study suggests that type BIIIb rotationplasty may be an alternative to amputation in patients with incurable infections.For malignant tumors of the lower extremities without invasion of the sciatic nerve,type BIIIb rotationplasty remains an excellent alternative to amputation.This surgical method may prevent amputation,improve functional outcomes,and facilitate biological reconstruction.
基金supported by the National Natural Science Foundation of China(11102237)Program for Changjiang Scholars and Innovative Research Team in University(IRT1294)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20110133120012)China Scholarship Council(CSC)
文摘Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.
基金supported by the Open Fund Project of State Key Laboratory of Lithospheric Evolution [SKL-K202103]support of the Exploration and Development Research Institute of Petro China Southwest Oil & Gas Field
文摘Hydrocarbon expulsion features and resource potential evaluation of source rocks are crucial for the petroleum exploration.High-maturity marine source rocks have not exhibited a hydrocarbon expulsion mode owing to the lack of low-maturity source rocks in deep petroliferous basins.We considered the Ediacaran microbial dolomite in the Sichuan Basin,the largest high-maturity marine gas layer in China,to exhibit a method that quantitatively characterizes the hydrocarbon expulsion of high-maturity marine source rocks.The experiment of fluid inclusion,rock pyrolysis,and vitrinite reflectance(Ro)of 119 microbial dolomite core samples obtained from the Dengying Formation were performed.A hydrocarbon expulsion model of high-maturity source rock was established,and its resource potential was evaluated.The results showed that the Ediacaran microbial dolomite in the Sichuan Basin is a good source rock showing vast resource potential.The hydrocarbon expulsion threshold is determined to be vitrinite reflectance at 0.92%.The hydrocarbon expulsion intensities in the geologic history is high with maximum of 1.6×10^(7)t/km^(2).The Ediacaran microbial dolomite expelled approximately 1.008×10^(12)t of hydrocarbons,and the recoverable resource was 1.5×10^(12)m^(3).The region can be categorized into areasⅠ,Ⅱ,Ⅲ,andⅣ,in decreasing order of hydrocarbon expulsion intensity.Areas with a higher hydrocarbon expulsion intensity have a lower drilling risk and should be prioritized for exploration in the orderⅠ>Ⅱ>Ⅲ>Ⅳ.Two areas,northern and central parts of Ediacaran in the Sichuan Basin,were selected as prospects which had the drilling priority in the future gas exploration.The production data of 55 drilled wells verified the high reliability of this method.This model in this study does not require low-maturity samples and can be used for evaluating high-maturity marine source rocks,which has broad applicability in deep basins worldwide.
基金Supported by Grant from Medical Technology Innovation Project of Nanjing Military,No. 09MA066
文摘AIM: TO screen the differential expressed genes in colorectal cancer and polyp tissue samples. METHODS: Tissue specimens containing 16 cases of colorectal adenocarcinoma and colorectal polyp vs nor- mal mucosae were collected and subjected to cDNA microarray and bioinformatical analyses. Quantitative reverse transcription-polymerase chain reaction (qRT- PCR) was used to confirm some of the cDNA microarray data.RESULTS: The experimental data showed that eight genes were differentially expressed, most of which were upregulated in adenomatous polyp lesions. Forty-six genes expressions were altered in colorectal cancers, of which 29 were upregulated and 17 downregulated, as compared to the normal mucosae. In addition, 18 genes were similarly altered in both adenomatous polyps and colorectal cancer, qRT-PCR analyses confirmed the cDNA microarray data for four of those 18 genes: MTA1, PDCD4, TSC1 and PDGFRA. CONCLUSION: These differentially expressed genes likely represent biomarkers for early detection of co- Iorectal cancer and may be potential therapeutic targets after confirmed by further studies.
文摘In this paper we consider mixed finite element methods for second order elliptic problems. In the case of the lowest order Brezzi-Douglas-Marini elements (if d = 2) or Brezzi- Douglas-Duran-Fortin elements (if d = 3) on rectangular parallelepipeds, we show that the mixed method system, by incorporating certain quadrature rules, can be written as a simple, cell-centered finite difference method. This leads to the solution of a sparse, positive semidefinite linear system for the scalar unknown. For a diagonal tensor coefficient, the sparsity pattern for the scalar unknown is a five point stencil if d = 2, and seven if d = 3. For a general tensor coefficient, it is a nine point stencil, and nineteen, respectively. Applications of the mixed method implementation as finite differences to nonisothermal multiphase, multicomponent flow in porous media are presented.