Các phản ứng hóa học có thể xảy ra trong quá trình dịch chuyển nước bơm ép trong vỉa (như sa lắng, thành đá, hòa tan, trao đổi cation…) thường không được tính đến trong quá trình mô phỏng khai thác bằng các phần mềm chuyên dụng. Bài báo giới thiệu nguyên lý và phần mềm kết nối giữa mô hình địa hóa mã nguồn mở (PHREEQC) và công cụ mô phỏng khai thác, để xây dựng mô hình mô phỏng tối ưu nhằm đánh giá quá trình dịch chuyển nước bơm ép từ giếng bơm ép đến giếng khai thác. Sự thay đổi dọc đường dịch chuyển của nước vỉa/nước bơm ép do các tương tác địa hóa được tính toán và mô hình hóa dưới dạng 2D/3D. Mô hình kết nối cũng được sử dụng để đánh giá thử nghiệm sự dịch chuyển của lưu thể vỉa giữa giếng bơm ép và giếng khai thác tại một mỏ thực tế.
THĂM DỊ - KHAI THÁC DẦU KHÍ TẠP CHÍ DẦU KHÍ Số 12 - 2019, trang 20 - 27 ISSN-0866-854X KẾT NỐI THƯ VIỆN DỮ LIỆU MƠ HÌNH ĐỊA HĨA MÃ NGUỒN MỞ (PHREEQC) VỚI MƠ HÌNH MƠ PHỎNG KHAI THÁC ĐỂ MƠ PHỎNG CÁC TƯƠNG TÁC ĐỊA HĨA VÀ SỰ DỊCH CHUYỂN CỦA NƯỚC BƠM ÉP TRONG VỈA DẦU KHÍ Kiều Anh Trung1, Đồn Huy Hiên1, Phạm Q Ngọc1, Hà Thu Hương1, Hoàng Long1 Lê Thị Thu Hường1, Nguyễn Minh Quý1, Ngô Hồng Anh1, Phạm Thị Thúy2 Viện Dầu khí Việt Nam Đại học Mỏ - Địa chất Email: trungka.epc@vpi.pvn.vn Tóm tắt Các phản ứng hóa học xảy q trình dịch chuyển nước bơm ép vỉa (như sa lắng, thành đá, hòa tan, trao đổi cation…) thường khơng tính đến q trình mơ khai thác phần mềm chuyên dụng Bài báo giới thiệu nguyên lý phần mềm kết nối mơ hình địa hóa mã nguồn mở (PHREEQC) công cụ mô khai thác, để xây dựng mơ hình mơ tối ưu nhằm đánh giá trình dịch chuyển nước bơm ép từ giếng bơm ép đến giếng khai thác Sự thay đổi dọc đường dịch chuyển nước vỉa/nước bơm ép tương tác địa hóa tính tốn mơ hình hóa dạng 2D/3D Mơ hình kết nối sử dụng để đánh giá thử nghiệm dịch chuyển lưu thể vỉa giếng bơm ép giếng khai thác mỏ thực tế Từ khóa: PHREEQC, phản ứng địa hóa, bơm ép nước Giới thiệu Bơm ép nước giải pháp hiệu để tăng lưu lượng giếng khai thác, ổn định tỷ số khí - dầu, nâng cao hệ số thu hồi dầu thơng qua việc trì áp suất vỉa nâng cao hiệu ứng đẩy quét Một số công nghệ bơm ép như: bơm ép thông minh, bơm ép nước kết hợp bơm ép khí ngày nghiên cứu ứng dụng nhiều Tuy nhiên, mô trình bơm ép nước, thay đổi thành phần hóa học nước bơm ép tương tác nước bơm ép với nước vỉa thường không tính đến Các phần mềm mơ khai thác thông dụng (như ECLIPSE, UTCHEM, CMG) mô túy trình dịch chuyển chất lưu mà khơng tính tốn đến trao đổi tương tác hệ chất lưu thông qua phản ứng hóa học [1] Các phản ứng hóa học xảy trình dịch chuyển nước vỉa, nước bơm ép ảnh hưởng lớn đến thay đổi dòng chảy, thay đổi độ rỗng, góc thấm ướt lưu lượng khai thác dầu khí Việc tính tốn ảnh hưởng phản ứng hóa học q trình mơ dòng chảy lưu thể vỉa cho phép xác hóa q trình mơ khai thác, đặc Ngày nhận bài: 11/11/2019 Ngày phản biện đánh giá sửa chữa: 11 - 14/11/2019 Ngày báo duyệt đăng: 6/12/2019 20 DẦU KHÍ - SỐ 12/2019 biệt mỏ áp dụng biện pháp giúp gia tăng hệ số thu hồi dầu như: bơm ép nước, bơm ép nước thông minh, bơm ép ASP, bơm ép khí nước ln phiên [2] Bên cạnh đó, mơ khai thác dầu khí, việc khớp lịch sử thường quan tâm thông số khai thác như: lưu lượng khí (gas rate), lưu lượng dầu (oil rate) áp suất đáy giếng (bottom-hole pressure) Việc cố gắng khớp lịch sử thành phần nước vỉa khai thác giải thích biến đổi thành phần nước khai thác vỉa ảnh hưởng phản ứng hóa học thường khơng ý đến Việc nghiên cứu chuyên sâu nước khai thác ứng dụng giải pháp công nghệ thông tin cho phép tái tạo tranh toàn cảnh ảnh hưởng trình bơm ép giếng khai thác, từ góp phần nâng cao hiệu khai thác dầu khí Các phần mềm mơ thường dựa phương trình dòng chảy (như ECLIPSE, CMG) mơ tốt q trình dịch chuyển lưu thể cấu trúc rỗng không tính đến việc thay đổi thành phần hóa học lưu thể, đặc biệt nước tương tác hóa học Ngược lại, mơ hình nhiệt động học mơ phản ứng địa hóa (như Multiscale, Scalechem, Geochemist) mơ tả tốt phản ứng hóa học xảy điều kiện bề mặt mà điều kiện vỉa PETROVIETNAM mô biến đổi pha, lưu lượng dịch chuyển hệ thống rỗng Do vậy, việc tìm phương pháp kết hợp phần mềm mơ vừa mơ q trình khai thác dầu khí, dịch chuyển chất lưu hệ thống rỗng vừa mô tương tác địa hóa điều kiện vỉa nhiệm vụ cần giải để xác hóa ảnh hưởng nước bơm ép giếng khai thác dầu khí Một mơ hình mơ dịch chuyển tương tác hóa học (reactive transport modelling) nội dung nghiên cứu trình bày báo Phương trình vận chuyển ion hòa tan theo dòng chảy =Buckley-Leverett Phương trình đối ,với trường hợp , − , , − nước chứa hàm lượng định ion chất hòa tan (1) − , − , [4 - 6]) sau:, = = ,, − , − , , = , , (1) =0 Trong đó: u: Vận tốc Darcy định nghĩa lưu lượng − tích; �, fw, Sw độ rỗng, dòng chảy đơn vị=diện tỷ phần dòng chảy nước độ bão hòa nước tương ứng Tỷ phần dòng chảy tỷ số lưu lượng nước (qw) lưu lượng tổng (qt) = − − , , , Lời giải số phương trình (1) trình bày + = Hình Với lưới thơng qua lưới tính tốn thể + i, lấy tích =0 phân vế, phương tính tốn này, ơ, lưới + = , − = − − , , trình (1) trở thành: =− = =− − (2) Rời rạc hóa tích phân vế phương trình (2) theo chiều khơng gian, ta có: = − − , , , = − − (3) = − − , , ,, , , Phương trình đương với: − tương , (3) , = − − , − ,, − ,, = − − = − − ,, , (4) ,, , (6) =, − = Sử dụng khái niệm + velocity) +vận tốc hạt (particle = phương trình (6) viết + lạiℎnhư sau: = + = (7) + ℎ = + ℎ = + là: Phương trình dòng chảy chiều pha Buckley-Leverett (BL) [3] sau: + −+ + + + , = , , = , ℎ = = = = (8) − , − , + = Phương trình vi phân đạo hàm riêng (8) = ,ẩn sau: − , , − , giải dạng hàm = , − , , − , − , − (9) , Phương trình (9) trình bày trình vận chuyển ion nước Kết hợp phản ứng hóa học q trình mơ dòng chảy lưu thể vỉa Khi ion hóa học dung dịch tương tác với nhau, tương tác, trao đổi cation với đất đá vỉa xung quanh, trộn lẫn với tạo thành dung dịch có nồng độ khác trạng thái cân (Hình 3) Việc tích hợp phương trình dòng chảy thực q trình hóa học quan trọng việc đánh giá định lượng ảnh hưởng giếng bơm ép lên giếng khai thác phương pháp phân tích hóa học [7] Trong trình dịch chuyển phản ứng sau thường xảy nước bơm ép đất đá vỉa cần tính đến mơ [9] Phương trình hàm độ bão hòa nước: , = , − , − (5) , L Δx = NL Dòng vào xi - xi + N + N Dòng = Đá Nước Calcite KV sét CO2 Ranh giới cell Bước đệm x0 x1 x2 xi-1 xi xi + Hình Lưới tính tốn + choℎlời giải số=phương trình (1) xN - xN Hình Mơ hình tương tác đất đá cell tính tốn [8] DẦU KHÍ - SỐ 12/2019 + = 21 THĂM DỊ - KHAI THÁC DẦU KHÍ • Các phản ứng hóa học mơi trường nước: H2 O H+ + OH-; K = 10-14,0 Ba2+ + CO32CO32- + H+ Ba2+ + HCO3- BaCO3; K = 102,71 Ba2+ + H2O HCO3-; K = 1010,329 BaOH+ + H+; K = 10-13,47 Mg2+ + CO32-+ H+ Ca2+ + CO32- + H+ CaHCO3+; K = 1011,435 Na+ + HCO3- Sr2+ + CO32- + H+ SrHCO3+; K = 1011,509 Na+ + CO32- Ca2+ + CO32- CaCO3; K = 103,244 Mg2+ + CO32- Sr2+ + CO32- SrCO3; K = 102,81 Ca2+ + H2O Mg2+ + H2O MgOH+ + H+; K = 10-11,44 Na+ + OH- Sr2+ + H2O BaHCO3+; K = 100,982 MgHCO3+; K = 1011,399 NaHCO3; K = 10-0,25 NaCO3-; K = 101,27 MgCO3; K = 102,98 CaOH+ + H+; K = 10-12,78 NaOH; K = 10-10 SrOH+ + H+; K = 10-13,29 • Các phản ứng hóa học hòa tan hình thành đá/sa lắng: Aragonite Strontianite Ca2+ + CO32- CaCO3 KSP = 10-8,336 Calcite KSP = 10-9,271 Sylvite Ca2+ + CO32- CaCO3 KSP = 10-8,48 Celestite Dolomite KSP = 10 NaCl K+ + Cl- CaMg(CO3)2 Ca2+ + Mg2+ + CO32- KSP = 10 -6,119 Halite KCl KSP = 100,9 Sr2+ + SO42- SrSO4 Sr2+ + CO32- SrCO3 -17,09 Witherite Na+ + Cl- KSP = 101,57 BaCO3 Ba2+ + CO32- KSP = 10-8,562 • Các phản ứng trao đổi cation nước đất đá vỉa: Na+ + X- NaX Kex = 100,9 Ca2+ + 2X- KX Kex = 100,7 CaX2 Kex = 100,8 Ba2+ + 2X- K + + X- Mg2+ + 2X- MgX2 Kex = 100,6 BaX2 Kex = 100,91 Sr2+ + 2X- SrX2 Kex = 100,8 Một số phần mềm địa hóa mơ tương tác địa hóa MINEQL+ (Schecher McAvoy, 1992), Geochemist’s Workbench (Bethke Yeakel, 2009), PHREEQC [10] Tuy nhiên, nghiên cứu sử dụng mơ hình địa hóa mã nguồn mở PHREEQC cung cấp khả tính tốn cho mơ hình hóa vận chuyển kết hợp tương tác hóa học, đặc biệt mơi trường rỗng với có mặt pha dầu, khí nước PHREEQC (pH-Redoc-Equilibrium in C programming language) phát triển Hiệp hội Địa chất Hoa Kỳ 22 DẦU KHÍ - SỐ 12/2019 (United States Geological - Survey USGS) Đây công cụ hữu ích để mô nghiên cứu vận chuyển tương tác hóa học với thư viện liệu phong phú PHREEQC có khả tính tốn số trạng thái bão hòa (saturation index) đặc trưng hóa học hợp chất ion nước, tính tốn cân phản ứng (reactions balance) tính tốn dịch chuyển mơ hình chiều (1D reactive transport) hợp phần kết hợp với phản ứng thuận nghịch bất thuận nghịch hòa tan nước, kết tủa, trao đổi ion, tương tác bề mặt, chuyển khối, phối trộn mơ hình nghịch đảo [10] Mơ hình tốn học sử dụng để xác định đặc trưng hóa học PHREEQC phương trình cân phản ứng dạng mole (mole balance equations), mơ hình số hoạt độ (activity coefficient model) giá trị tích số tan Mơ hình số hoạt độ mơ hình sử dụng để mô tả mối quan hệ số hoạt độ hợp chất lực ion dung dịch Đặc điểm bật PHREEQC là: PETROVIETNAM - Có thể tính tốn phản ứng hóa học dung dịch có tiếp xúc với pha khí; - Có thể tính tốn cân động học phản ứng tương tác rắn/lỏng; - Có thư viện đầy đủ mơ hình số hoạt độ nhiều phản ứng nhiệt độ, áp suất khác Sử dụng PHREEQC mơ tốt tương tác địa hóa mơ hình vận chuyển chiều (1D), nhiên áp dụng cho mơ hình dòng chảy đa pha 2D/3D (dầu khí) có số hạn chế Do đó, cần kết nối PHREEQC với mơ hình mơ vỉa chun biệt để xác hóa tương tác địa hóa xảy vỉa độ sâu lớn, nhiệt độ áp suất lớn Sử dụng Matlab công cụ dạng mã nguồn mở (như MRST - Matlab Reservoir Simulation Toolbox [11]) kết nối với PHREEQC giải tốn mơ hình dòng chảy đa pha 2D/3D mơi trường vi rỗng, đồng thời xuất hiển thị kết cách dễ dàng Kết nối PHREEQC với mơ hình mơ khai thác Để kết nối PHREEQC với mơ hình mơ khai thác, nghiên cứu sử dụng công cụ IPhreeC, phát triển dựa PHREEQC IPhreeQC [12] công cụ giúp tích hợp hàm, thư viện PhreeQC vào ngơn ngữ lập trình (như C, C2+) gọi từ phần mềm ứng dụng windows (Microsoft Excel, Matlab ) thông qua cổng giao tiếp COM module [13] Điều cho phép việc kết hợp phương trình dòng chảy ion với q trình địa hóa thực dễ dàng Quy trình kết nối PHREEQC với mơ hình mô khai thác sau: Khi kết nối PHREEQC với mơ hình mơ khai thác, biến đổi ion nước dịch chuyển có khác biệt lớn với mơ hình mơ túy khơng tính đến tương tác hóa học Hình thể kết mô dịch chuyển nước bơm ép trường hợp mơ hình có đặc trưng rỗng, thấm Áp dụng mơ q trình dịch chuyển nước bơm ép vỉa Sau kết nối thành công PHREEQC với công cụ MRST, nhóm nghiên cứu áp dụng để mơ dịch chuyển nước bơm ép tương tác với đất đá vỉa, nước vỉa dọc đường dịch chuyển đến giếng khai thác Các thông số vỉa liệu đầu vào thể Bảng Kết mô dạng lát cắt 2D thể Hình (độ bão hòa nước dịch chuyển ion Cl-) kết khớp lịch sử thành phần hóa học số ion đặc trưng nước khai thác giếng khai thác thể Hình Đọc thơng số đầu vào vỉa, chất lưu (độ rỗng, độ thấm, áp suất, nhiệt độ, độ nhớt…), giếng khoan Khởi tạo mơ hình với thông số đầu vào ban đầu Đọc database (các số phản ứng tốc độ phản ứng), hàm lượng nước bơm ép, nước vỉa, khống vật phản ứng hóa học Tiến hành tính tốn để xác định thành phần hóa học nước Khởi tạo hệ hóa học tất ô lưới Thực mô với vài bước thời gian điều kiện ban đầu có khơng có giếng khai thác/bơm ép để đảm bảo hợp lý dòng chảy hệ hóa học Thực mô bước theo thời gian Vận chuyển thành phần hóa học pha nước Tại ô lưới gọi hàm PhreeQC để thực phản ứng hóa học (khơng thay đổi nhiệt độ, áp suất độ bão hòa nước) No Cập nhật lại, tổng tất thành phần hóa học Kết thúc mơ Hình Quy trình kết nối PHREEQC với mơ hình mơ khai thác DẦU KHÍ - SỐ 12/2019 23 THĂM DỊ - KHAI THÁC DẦU KHÍ Phản ứng Không phản ứng Phản ứng Không phản ứng Phản ứng Khơng phản ứng Phản ứng Khơng phản ứng Hình So sánh kết mơ q trình vận chuyển ion dòng chảy với tương tác khơng tương tác hóa học suốt q trình bơm ép Bảng Thông số đầu vào cho mô Độ rỗng trung bình 0,17 Hệ số Corey - nước, - dầu Độ nhớt (Pa.s) 1e-3 - nước, 2e-3 - dầu Độ bão hòa dầu dư 0,2 Độ bão hòa nước dư 0,2 Độ bão hòa nước ban đầu 0,2 Độ thấm trung bình (mD) 200 Diện tích vỉa (m2) 600 × 2.200m Độ lớn lưới 10m Độ thấm tương đối cuối dầu Độ thấm tương đối cuối nước Hàm lượng ion nước (mol/kgw) Na+ Nước vỉa 0,46 0,006 Ca 0,01 0,006 Mg2+ 0,053 0,0006 2+ 24 0,25 Nước bơm ép Cl- 0,54 0,07 HCO3- 0,001 0,0001 SO42- 0,027 0,0002 DẦU KHÍ - SỐ 12/2019 PETROVIETNAM Hình Mặt cắt 2D mơ độ bão hòa nước theo thời gian mơ Hình Mặt cắt 2D mô thay đổi hàm lượng Cl- nước bơm ép theo thời gian mô DẦU KHÍ - SỐ 12/2019 25 Mol/kg nước Mol/kg nước THĂM DỊ - KHAI THÁC DẦU KHÍ Số liệu mô Thực tế Mol/kg nước Mol/kg nước Số liệu mô Thực tế Số liệu mô Thực tế Số liệu mơ Thực tế Hình Khớp lịch sử thành phần ion mẫu nước khai thác giếng khai thác Sự kết nối PHREEQC với mơ hình mô khai thác cho phép mô tả thay đổi thành phần hóa học nước khai thác dọc đường dịch chuyển vỉa Bên cạnh phương pháp bơm ép nước truyền thống, khai thác dầu khí cần áp dụng phương pháp khác liên quan đến bơm ép có tương tác địa hóa bơm ép nước thơng minh, bơm ép khí nước luân phiên, bơm ép ASP, bơm ép surfactant Cần thiết phải mơ tả tương tác địa hóa nước bơm ép đất đá vỉa lưu thể khác vỉa, từ cho phép đánh giá ảnh hưởng trình bơm ép hiệu tăng cường thu hồi dầu Kết luận Nghiên cứu giới thiệu bước để kết nối mô hình địa hóa với mơ hình mơ khai thác Một phần mềm hoàn chỉnh dựa kết nối PHREEQC MRST cho phép mơ tả xác dịch chuyển dòng nước bơm 26 DẦU KHÍ - SỐ 12/2019 ép vỉa có tính đến tương tác địa hóa xảy nước đất đá vỉa ảnh hưởng áp suất nhiệt độ vỉa Mơ hình kết nối PHREEQC với mơ hình mơ khai thác khắc phục nhược điểm mơ hình mơ khai thác sử dụng mở hướng nghiên cứu chuyên sâu cho việc nghiên cứu trình bơm ép hóa học khai thác dầu khí Lời cám ơn Nhóm tác giả chân thành cảm ơn USGS cung cấp công cụ mã nguồn mở PHREEQC để sử dụng nghiên cứu Viện Dầu khí Việt Nam tài trợ kinh phí để thực nghiên cứu Tài liệu tham khảo Marco De Lucia, Michael Kühn Coupling R and PHREEQC: Efficient programming of geochemical models Energy Procedia, 2013; 40: p 464 - 471 PETROVIETNAM Ji Ho Lee, Kun Sang Lee Geochemical evaluation of low salinity hot water injection to enhance heavy oil recovery from carbonate reservoirs Petroleum Science 2019; 16: p 366 - 381 L.Dake Fundamentals of reservoir engineering (19 edition) Elsevier 1978 th Aruoture Voke Omekeh, Helmer A.Friis, Ingebret Fjelde, Steianar Evje Modeling of Ion-Exchange and solubility in low salinity water flooding SPE Improved Oil Recovery Symposium, Tulsa, Oklahoma, USA 14 - 18 April 2012 Eli L.Isaacson Global solution of a riemann problem for non strictly hyperbolic system of conservation laws arising in enhanced oil recovery Enhanced Oil Recovery Institute, University of Wyoming 1989 Deniz M.Dindoruk, Birol Dindoruk Analytical solution of nonisothermal Buckley-Leverett flow including tracers SPE Reservoir Evaluation & Engineering 2008: p 555 - 564 Tina Puntervold, Tor Austad Injection of seawater and mixtures with produced water into North Sea chalk formation: Impact of fluid-rock interactions on wettability and scale formation Journal of Petroleum Science and Engineering 2008; 63: p 23 - 33 J.B.Wouter Simulation of geochemical processes during low salinity water flooding by coupling multiphase Buckley-Leverett flow to the geochemical package PHREEQC MSc thesis TUDeft 2006 Aboulghasem Kazemi Nia Korrani, Kamy Sepehrnoori, Mojdeh Delshad Coupling IPhreeqc with UTCHEM to model reactive flow and transport Computers & Geosciences 2015; 82: p 152 - 169 10 D.L.Parkhurst, C.A.J Appelo User’s guide to PHREEQC (Version2) - A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations U.S Geological Survey WaterResources Investigations Report 99-4259 1999 11 Knut-Andreas Lie An Introduction to Reservoir Simulation Using MATLAB/GNU Octave: User guide for the MATLAB reservoir simulation toolbox (MRST) Cambridge University Press 2019 12 S.R.Charlton, D.L.Parkhurst Modules based on the geochemical model PHREEQC for use in scripting and programming languages Computer & Geosciences 2011; 37: p 1653 - 1664 13 Haishan Luo, Emad W.Al-Shalabi, Mojdeh Delshad, Krishna Panthi, Kamy Sepehrnoori A robust geochemical simulator to model improved oil recovery methods SPE Reservoir Simulation Symposium, Houston, Texas, USA 23 - 25 February, 2015 COUPLING BETWEEN PHREEQC AND RESERVOIR PRODUCTION MODEL TO SIMULATE THE REACTIVE TRANSPORT OF INJECTED WATER IN OIL AND GAS RESERVOIRS Kieu Anh Trung1, Doan Huy Hien1, Pham Quy Ngoc1, Ha Thu Huong1, Hoang Long1 Le Thi Thu Huong1, Nguyen Minh Quy1, Ngo Hong Anh1, Pham Thi Thuy2 Vietnam Petroleum Institute Hanoi University of Mining and Geology Email: trungka.epc@vpi.pvn.vn Summary Many chemical reactions occur during transport of fluid in oil reservoirs such as scaling, dissolution, evaporation and cation exchange However, these parameters are usually not considered in reservoir simulators This paper presents a method of coupling between geochemical code (PHREEQC) and a reservoir modelling tool This coupling creates a useful software to simulate the reactive transport of injected water from injection well to production well The change of water chemical composition due to interaction between fluid and rock, and between fluid and fluid will be calculated and visualised under 2D/3D graphic The coupling model is also tested with real data Key words: PHREEQC, geochemical reaction, injected water DẦU KHÍ - SỐ 12/2019 27 ... đến bơm ép có tương tác địa hóa bơm ép nước thơng minh, bơm ép khí nước luân phiên, bơm ép ASP, bơm ép surfactant Cần thiết phải mô tả tương tác địa hóa nước bơm ép đất đá vỉa lưu thể khác vỉa, ... để mô dịch chuyển nước bơm ép tương tác với đất đá vỉa, nước vỉa dọc đường dịch chuyển đến giếng khai thác Các thông số vỉa liệu đầu vào thể Bảng Kết mô dạng lát cắt 2D thể Hình (độ bão hòa nước. .. phép mô tả xác dịch chuyển dòng nước bơm 26 DẦU KHÍ - SỐ 12/2019 ép vỉa có tính đến tương tác địa hóa xảy nước đất đá vỉa ảnh hưởng áp suất nhiệt độ vỉa Mơ hình kết nối PHREEQC với mơ hình mơ khai