Cột dầu là bài toán hai gradient áp suất
Một áp suất mao dẫn 60 kPa không tự cho biết chiều cao cột dầu nếu chưa có khối lượng riêng của cả nước vỉa và dầu tại điều kiện đang xét. Lấy 60 kPa chia trọng lượng riêng nước chỉ tái hiện một cột nước so với pha có khối lượng riêng bỏ qua được; dầu trong lỗ rỗng không đáp ứng giả thiết đó.
Bài Fundamentals of Wettability định nghĩa mực nước tự do, FWL, tại nơi áp suất mao dẫn dầu–nước bằng không và phân biệt mốc này với tiếp xúc dầu–nước thể hiện qua phân bố bão hòa. [1] Bài dưới đây tự dựng các bộ số, giả định cân bằng tĩnh và hệ ưa nước. Kết quả là minh họa quan hệ vật lý, không phải xác định tiếp xúc hay khả năng giữ dầu của một cấu tạo thực.
Định nghĩa dấu trước khi viết công thức
Đặt z đo lên trên từ FWL, P_c=p_o−p_w, trong đó o là dầu không ưa bề mặt và w là nước ưa bề mặt. Ở trạng thái thủy tĩnh, dp_o/dz=−ρ_o g và dp_w/dz=−ρ_w g. Trừ hai phương trình:
dP_c/dz = (ρ_w−ρ_o)g.
Nếu khối lượng riêng không đổi và P_c=0 tại z=0, tích phân được P_c=Δρgh, hay h=P_c/(Δρg). P_c tính bằng Pa, Δρ=ρ_w−ρ_o bằng kg/m³, g bằng m/s², h bằng m. Giá trị h là cao độ theo phương trọng lực, không phải chiều dài dọc một vỉa nghiêng hoặc dọc lỗ khoan.
Điểm quan trọng của phép trừ là hai pha đều có gradient. Dầu càng gần khối lượng riêng nước thì cùng P_c cần một h càng lớn. Khi Δρ tiến về không, quan hệ h trở nên rất nhạy; đó là giới hạn của việc suy cao độ từ chênh áp, không phải bằng chứng có thể giữ một cột vô hạn trong địa chất thực.
Ví dụ với khối lượng riêng tại vỉa
Cho ρ_w=1.050 kg/m³, ρ_o=800 kg/m³, g=9,81 m/s². Gradient chênh áp là 2.452,5 Pa/m. Với P_c=60.000 Pa, ta được h=24,464832 m. Nếu nhầm dùng riêng ρ_w, kết quả chỉ còn 5,824960 m, nhỏ hơn đúng 4,2 lần.
Ngược lại, tại h=15 m, P_c=36.787,5 Pa. Hai phép tính thuận và nghịch là một kiểm tra đơn vị đơn giản. Nếu P_c được ghi bằng kPa nhưng chưa nhân 1.000, kết quả cao độ sẽ sai một nghìn lần. Nếu dữ liệu khối lượng riêng lấy ở bồn bề mặt trong khi áp suất, nhiệt độ vỉa khác, phép tính có thể nhất quán về đơn vị nhưng chưa đúng cơ sở vật lý.
Nối áp suất với hình học một ống mao dẫn
Với ống tròn lý tưởng bán kính trong r và góc tiếp xúc θ đo qua pha ưa bề mặt, cân bằng sức căng cho P_c=2σcosθ/r. σ là sức căng liên pha, N/m; r là m. OpenStax trình bày dạng chiều cao mao dẫn trong ống, cho thấy vai trò của góc tiếp xúc và bán kính. [2] Ở hệ hai chất lỏng, thay gradient trọng lượng bằng Δρg để nối với cao độ vừa suy ra.
Giả định σ=0,025 N/m, θ=30° và r=2 µm. Khi đổi 2 µm thành 2×10⁻⁶ m, P_c=21.650,635095 Pa và h=8,827986 m với hai khối lượng riêng ở trên. Đây là bán kính ống tương đương trong mô hình, không mặc nhiên là bán kính trung bình của mọi lỗ rỗng trong mẫu đá.
Nếu một báo cáo đưa đường kính cổ lỗ 4 µm, phải chia hai trước khi đưa vào r. Dùng 4 µm làm bán kính sẽ làm P_c và h giảm một nửa. Cũng cần phân biệt bán kính cổ lỗ với bán kính cong của mặt phân cách: quan hệ hình học chứa cosθ nối hai đại lượng, không nên thay chúng bằng nhau trong mọi góc tiếp xúc.
Độ nhạy đối với bề mặt và kích thước
Giữ r, σ và khối lượng riêng nhưng tăng θ từ 30° lên 60°, h giảm từ 8,827986 còn 5,096840 m. Tỷ số bằng cos60°/cos30°=0,577350. Với thay đổi nhỏ, d ln h=d ln σ−d ln r−d ln Δρ−tanθ dθ, trong đó dθ phải dùng radian.
Tăng θ thêm 5° từ điểm 30° làm h giảm chính xác 5,412469%; tuyến tính hóa cho khoảng 5,0383%. Khi góc tiến gần 90°, độ nhạy tương đối tăng mạnh. Qua 90°, dấu cos thay đổi; cần quay lại quy ước pha ưa bề mặt và chiều áp suất, thay vì lấy trị tuyệt đối để giữ một chiều cao dương thuận mắt.
Tăng r 20% làm h còn 1/1,2 giá trị cũ, giảm 16,6667%. Tăng sức căng liên pha 20% làm h tăng 20%. Các phép so này chỉ giữ những đại lượng khác cố định; trong thực tế, thay thành phần lưu chất có thể đồng thời đổi sức căng, góc tiếp xúc và khối lượng riêng.
Chuyển một phép đo phòng thí nghiệm sang hệ lưu chất khác
Giả định cùng hình học mao dẫn và cùng nhánh xâm nhập, hệ thí nghiệm có σ_lab=0,072 N/m, θ_lab=0° và đo P_c,lab=120 kPa. Hệ vỉa giả định có σ_res=0,025 N/m, θ_res=30°. Chia hai phương trình 2σcosθ/r, bán kính triệt tiêu:
P_c,res = P_c,lab × (σ_res cosθ_res)/(σ_lab cosθ_lab).
Kết quả là 36.084,391824 Pa, tương ứng h=14,713310 m. Không chuyển trực tiếp 120 kPa thành chiều cao vỉa trước khi xử lý khác biệt hệ lưu chất. Tuy nhiên, phép tỷ lệ này chỉ là suy luận từ hình học giữ nguyên; nếu chuẩn bị mẫu làm thay bề mặt, ứng suất làm đổi cổ lỗ hoặc lịch sử xâm nhập khác, cần dữ liệu thích hợp hơn để chuyển đường cong.
Tính thấm ướt mô tả ưu tiên tiếp xúc bề mặt, không phải chỉ báo rằng pha nào đang chiếm nhiều thể tích. Schlumberger cũng phân biệt bề mặt trung gian với hệ có nhiều trạng thái thấm ướt cùng tồn tại. [1] Vì vậy một góc đại diện duy nhất là giả thiết cần xem xét khi mẫu có cấu trúc phức tạp.
Từ một cổ lỗ sang đường bão hòa theo chiều cao
Để minh họa phân bố, chọn một đường Brooks–Corey dạng P_c=P_e S_e^(−1/λ), với S_e=(S_w−S_wr)/(1−S_wr). MOOSE trình bày quan hệ này và ý nghĩa áp suất vào P_e. [3] Chọn P_e=8.000 Pa, λ=2 và S_wr=0,20; mọi tham số là giả định, không lấy từ phép thử lõi.
Ở vùng P_c≥P_e, đảo hàm cho S_w=0,20+0,80(P_e/P_c)². Với P_c=Δρgh, cao độ bắt đầu xâm nhập theo đường giả định là h_e=3,261978 m. Dưới h_e, ta đặt S_w=1 cho nhánh minh họa này, không ngoại suy lũy thừa ra giá trị lớn hơn một.
Tại h=5; 15; 30 m, S_w lần lượt bằng 0,540496; 0,237833; 0,209458. Một mốc P_c=0, một ngưỡng dầu bắt đầu xâm nhập và một mốc bão hòa được chọn để gọi tiếp xúc là ba định nghĩa khác nhau. Không nên ép chúng trùng nhau chỉ để đơn giản hóa hình vẽ phân lớp.
Khi gradient thay đổi theo chiều cao
Nếu Δρ thay đổi, công thức đúng là P_c(h)=g∫₀ʰΔρ(z)dz. Chọn riêng một mô hình thử Δρ(z)=250+0,2z, với hệ số 0,2 mang đơn vị kg/m⁴. Khi đó P_c/g=250h+0,1h². Với 60 kPa, nghiệm dương h=24,229995 m, khác 24,464832 m của mật độ không đổi.
Ví dụ cho thấy phải tích phân gradient khi tính chất biến thiên đáng kể. Ngoài ra, nếu đang có dòng chảy, gradient nhớt và không cân bằng có thể làm chênh áp không còn chỉ do trọng lực. Trước khi chuyển đường mao dẫn thành phân bố dầu, cần xác định điều kiện tĩnh, nhánh lịch sử, nhóm đá, tính chất lưu chất và mốc cao độ. Những thông tin ấy quyết định ý nghĩa của đường cong nhiều hơn số chữ số thập phân được in ra.
Nguồn tham khảo
[1] Fundamentals of Wettability — Oilfield Review, Summer 2007, pp.44–61 — Wael Abdallah et al.; Schlumberger Oilfield Review. Truy cập ngày 09/11/2026.
[2] College Physics — 11.8 Cohesion and Adhesion in Liquids: Surface Tension and Capillary Action — Paul Peter Urone and Roger Hinrichs; OpenStax, Rice University. Truy cập ngày 09/11/2026.
[3] Capillary pressure — PorousFlow module — Idaho National Laboratory, MOOSE Framework. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một áp suất mao dẫn 60 kPa không tự cho biết chiều cao cột dầu nếu chưa có khối lượng riêng của cả nước vỉa và dầu tại điều kiện đang xét. Lấy 60 kPa chia trọng lượng riêng nước chỉ tái hiện một cột nước so với pha có khối lượng riêng bỏ qua được; dầu trong lỗ rỗng không đáp ứng giả thiết đó.
Bài Fundamentals of Wettability định nghĩa mực nước tự do, FWL, tại nơi áp suất mao dẫn dầu–nước bằng không và phân biệt mốc này với tiếp xúc dầu–nước thể hiện qua phân bố bão hòa. [1] Bài dưới đây tự dựng các bộ số, giả định cân bằng tĩnh và hệ ưa nước. Kết quả là minh họa quan hệ vật lý, không phải xác định tiếp xúc hay khả năng giữ dầu của một cấu tạo thực.
Định nghĩa dấu trước khi viết công thức
Đặt z đo lên trên từ FWL, P_c=p_o−p_w, trong đó o là dầu không ưa bề mặt và w là nước ưa bề mặt. Ở trạng thái thủy tĩnh, dp_o/dz=−ρ_o g và dp_w/dz=−ρ_w g. Trừ hai phương trình:
dP_c/dz = (ρ_w−ρ_o)g.
Nếu khối lượng riêng không đổi và P_c=0 tại z=0, tích phân được P_c=Δρgh, hay h=P_c/(Δρg). P_c tính bằng Pa, Δρ=ρ_w−ρ_o bằng kg/m³, g bằng m/s², h bằng m. Giá trị h là cao độ theo phương trọng lực, không phải chiều dài dọc một vỉa nghiêng hoặc dọc lỗ khoan.
Điểm quan trọng của phép trừ là hai pha đều có gradient. Dầu càng gần khối lượng riêng nước thì cùng P_c cần một h càng lớn. Khi Δρ tiến về không, quan hệ h trở nên rất nhạy; đó là giới hạn của việc suy cao độ từ chênh áp, không phải bằng chứng có thể giữ một cột vô hạn trong địa chất thực.
Ví dụ với khối lượng riêng tại vỉa
Cho ρ_w=1.050 kg/m³, ρ_o=800 kg/m³, g=9,81 m/s². Gradient chênh áp là 2.452,5 Pa/m. Với P_c=60.000 Pa, ta được h=24,464832 m. Nếu nhầm dùng riêng ρ_w, kết quả chỉ còn 5,824960 m, nhỏ hơn đúng 4,2 lần.
Ngược lại, tại h=15 m, P_c=36.787,5 Pa. Hai phép tính thuận và nghịch là một kiểm tra đơn vị đơn giản. Nếu P_c được ghi bằng kPa nhưng chưa nhân 1.000, kết quả cao độ sẽ sai một nghìn lần. Nếu dữ liệu khối lượng riêng lấy ở bồn bề mặt trong khi áp suất, nhiệt độ vỉa khác, phép tính có thể nhất quán về đơn vị nhưng chưa đúng cơ sở vật lý.
Nối áp suất với hình học một ống mao dẫn
Với ống tròn lý tưởng bán kính trong r và góc tiếp xúc θ đo qua pha ưa bề mặt, cân bằng sức căng cho P_c=2σcosθ/r. σ là sức căng liên pha, N/m; r là m. OpenStax trình bày dạng chiều cao mao dẫn trong ống, cho thấy vai trò của góc tiếp xúc và bán kính. [2] Ở hệ hai chất lỏng, thay gradient trọng lượng bằng Δρg để nối với cao độ vừa suy ra.
Giả định σ=0,025 N/m, θ=30° và r=2 µm. Khi đổi 2 µm thành 2×10⁻⁶ m, P_c=21.650,635095 Pa và h=8,827986 m với hai khối lượng riêng ở trên. Đây là bán kính ống tương đương trong mô hình, không mặc nhiên là bán kính trung bình của mọi lỗ rỗng trong mẫu đá.
Nếu một báo cáo đưa đường kính cổ lỗ 4 µm, phải chia hai trước khi đưa vào r. Dùng 4 µm làm bán kính sẽ làm P_c và h giảm một nửa. Cũng cần phân biệt bán kính cổ lỗ với bán kính cong của mặt phân cách: quan hệ hình học chứa cosθ nối hai đại lượng, không nên thay chúng bằng nhau trong mọi góc tiếp xúc.
Độ nhạy đối với bề mặt và kích thước
Giữ r, σ và khối lượng riêng nhưng tăng θ từ 30° lên 60°, h giảm từ 8,827986 còn 5,096840 m. Tỷ số bằng cos60°/cos30°=0,577350. Với thay đổi nhỏ, d ln h=d ln σ−d ln r−d ln Δρ−tanθ dθ, trong đó dθ phải dùng radian.
Tăng θ thêm 5° từ điểm 30° làm h giảm chính xác 5,412469%; tuyến tính hóa cho khoảng 5,0383%. Khi góc tiến gần 90°, độ nhạy tương đối tăng mạnh. Qua 90°, dấu cos thay đổi; cần quay lại quy ước pha ưa bề mặt và chiều áp suất, thay vì lấy trị tuyệt đối để giữ một chiều cao dương thuận mắt.
Tăng r 20% làm h còn 1/1,2 giá trị cũ, giảm 16,6667%. Tăng sức căng liên pha 20% làm h tăng 20%. Các phép so này chỉ giữ những đại lượng khác cố định; trong thực tế, thay thành phần lưu chất có thể đồng thời đổi sức căng, góc tiếp xúc và khối lượng riêng.
Chuyển một phép đo phòng thí nghiệm sang hệ lưu chất khác
Giả định cùng hình học mao dẫn và cùng nhánh xâm nhập, hệ thí nghiệm có σ_lab=0,072 N/m, θ_lab=0° và đo P_c,lab=120 kPa. Hệ vỉa giả định có σ_res=0,025 N/m, θ_res=30°. Chia hai phương trình 2σcosθ/r, bán kính triệt tiêu:
P_c,res = P_c,lab × (σ_res cosθ_res)/(σ_lab cosθ_lab).
Kết quả là 36.084,391824 Pa, tương ứng h=14,713310 m. Không chuyển trực tiếp 120 kPa thành chiều cao vỉa trước khi xử lý khác biệt hệ lưu chất. Tuy nhiên, phép tỷ lệ này chỉ là suy luận từ hình học giữ nguyên; nếu chuẩn bị mẫu làm thay bề mặt, ứng suất làm đổi cổ lỗ hoặc lịch sử xâm nhập khác, cần dữ liệu thích hợp hơn để chuyển đường cong.
Tính thấm ướt mô tả ưu tiên tiếp xúc bề mặt, không phải chỉ báo rằng pha nào đang chiếm nhiều thể tích. Schlumberger cũng phân biệt bề mặt trung gian với hệ có nhiều trạng thái thấm ướt cùng tồn tại. [1] Vì vậy một góc đại diện duy nhất là giả thiết cần xem xét khi mẫu có cấu trúc phức tạp.
Từ một cổ lỗ sang đường bão hòa theo chiều cao
Để minh họa phân bố, chọn một đường Brooks–Corey dạng P_c=P_e S_e^(−1/λ), với S_e=(S_w−S_wr)/(1−S_wr). MOOSE trình bày quan hệ này và ý nghĩa áp suất vào P_e. [3] Chọn P_e=8.000 Pa, λ=2 và S_wr=0,20; mọi tham số là giả định, không lấy từ phép thử lõi.
Ở vùng P_c≥P_e, đảo hàm cho S_w=0,20+0,80(P_e/P_c)². Với P_c=Δρgh, cao độ bắt đầu xâm nhập theo đường giả định là h_e=3,261978 m. Dưới h_e, ta đặt S_w=1 cho nhánh minh họa này, không ngoại suy lũy thừa ra giá trị lớn hơn một.
Tại h=5; 15; 30 m, S_w lần lượt bằng 0,540496; 0,237833; 0,209458. Một mốc P_c=0, một ngưỡng dầu bắt đầu xâm nhập và một mốc bão hòa được chọn để gọi tiếp xúc là ba định nghĩa khác nhau. Không nên ép chúng trùng nhau chỉ để đơn giản hóa hình vẽ phân lớp.
Khi gradient thay đổi theo chiều cao
Nếu Δρ thay đổi, công thức đúng là P_c(h)=g∫₀ʰΔρ(z)dz. Chọn riêng một mô hình thử Δρ(z)=250+0,2z, với hệ số 0,2 mang đơn vị kg/m⁴. Khi đó P_c/g=250h+0,1h². Với 60 kPa, nghiệm dương h=24,229995 m, khác 24,464832 m của mật độ không đổi.
Ví dụ cho thấy phải tích phân gradient khi tính chất biến thiên đáng kể. Ngoài ra, nếu đang có dòng chảy, gradient nhớt và không cân bằng có thể làm chênh áp không còn chỉ do trọng lực. Trước khi chuyển đường mao dẫn thành phân bố dầu, cần xác định điều kiện tĩnh, nhánh lịch sử, nhóm đá, tính chất lưu chất và mốc cao độ. Những thông tin ấy quyết định ý nghĩa của đường cong nhiều hơn số chữ số thập phân được in ra.
Nguồn tham khảo
[1] Fundamentals of Wettability — Oilfield Review, Summer 2007, pp.44–61 — Wael Abdallah et al.; Schlumberger Oilfield Review. Truy cập ngày 09/11/2026.
[2] College Physics — 11.8 Cohesion and Adhesion in Liquids: Surface Tension and Capillary Action — Paul Peter Urone and Roger Hinrichs; OpenStax, Rice University. Truy cập ngày 09/11/2026.
[3] Capillary pressure — PorousFlow module — Idaho National Laboratory, MOOSE Framework. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
