Một mẫu lõi có thể cho nhiều giá trị độ thấm khí
Nếu hai phòng thí nghiệm dùng cùng mẫu nhưng khác áp suất đo, các số độ thấm khí có thể khác dù cấu trúc mẫu chưa thay đổi. Có ít nhất hai việc phải tách: lưu lượng thể tích thay đổi khi khí giãn nở dọc mẫu, và khả năng dòng khí trượt tại thành lỗ rỗng làm độ thấm biểu kiến phụ thuộc áp suất.
DTU mô tả phép đo ở nhiều áp suất trung bình, sau đó hồi quy độ thấm theo nghịch đảo áp suất để xác định giá trị hiệu chỉnh Klinkenberg. [1] Bài này tự dựng một bộ số giả định để giải thích phép biến đổi, không cung cấp quy trình vận hành thiết bị chịu áp.
Từ Darcy tại một tiết diện đến công thức khí nén được
Xét lõi trụ dài L, tiết diện A, dòng một chiều ổn định, khí lý tưởng đẳng nhiệt, độ nhớt μ không đổi. Bỏ trọng lực, quán tính và rò quanh áo mẫu. Ban đầu giả sử độ thấm k không phụ thuộc vị trí hay áp suất. Darcy cục bộ viết q(x)=−kA(dp/dx)/μ, với q là lưu lượng thể tích tại đúng vị trí x.
Bảo toàn khối lượng cùng nhiệt độ không đổi cho p(x)q(x)=p_out q_out. Thay vào Darcy rồi tích phân từ đầu vào tới đầu ra:
q_out = k A (p_in²−p_out²)/(2 μ L p_out).
k_app = 2 μ L p_out q_out/{A(p_in²−p_out²)}.
p_in và p_out phải là áp suất tuyệt đối, đơn vị Pa; μ dùng Pa·s, L dùng m, A dùng m², q_out dùng m³/s thì k có đơn vị m². Lưu lượng tại đầu ra không được thay bằng lưu lượng quy chuẩn khi chưa quy đổi về p_out và nhiệt độ mẫu.
Đặt p_m=(p_in+p_out)/2 và Δp=p_in−p_out. Vì p_in²−p_out²=2p_mΔp, có thể viết k_app=μLq_out p_out/(A p_mΔp). Dạng này làm rõ hệ số hiệu chỉnh khí nén được p_out/p_m so với phép tính chất lỏng.
Mô hình trượt khí và giá trị ngoại suy
Báo cáo Kansas Geological Survey dùng quan hệ Klinkenberg k_g=k_liquid(1+b/P), đồng thời nêu hệ số b thay đổi theo cấu trúc lỗ rỗng. [2] Với ký hiệu trong bài, k_∞ là giá trị ngoại suy khi nghịch đảo áp suất tiến về không; b có đơn vị áp suất.
k_app = k_∞(1+b/p_m).
Có thể kiểm tính nhất quán bằng giả định cục bộ k(p)=k_∞(1+b/p). Khi tích phân Darcy, tích p k(p)=k_∞(p+b) cho:
p_out q_out = k_∞ A { (p_in²−p_out²)/2 + bΔp }/(μL).
Chia kết quả theo định nghĩa k_app ở trên thu đúng k_∞(1+b/p_m). Vì vậy trong các giả thiết này, áp suất trung bình số học xuất hiện từ tích phân, không phải một lựa chọn tùy ý. Nếu khí thực, nhiệt độ hay độ nhớt thay đổi đáng kể, tích phân cần các hàm tính chất phù hợp.
Dựng ba điểm đo giả định
Cho đường kính mẫu D=0,025 m, L=0,050 m, μ=1,8×10⁻⁵ Pa·s. Diện tích A=πD²/4=0,000490874 m². Chọn k_∞=20 mD và b=0,080 MPa; quy đổi trong ví dụ dùng 1 mD=9,869233×10⁻¹⁶ m².
Giữ Δp=0,020 MPa. Ba lần đo có p_m lần lượt 0,20; 0,40; 0,80 MPa tuyệt đối. Do đó các cặp p_in/p_out là 0,21/0,19; 0,41/0,39; 0,81/0,79 MPa. Theo mô hình, k_app lần lượt 28; 24; 22 mD.
Thay các giá trị vào công thức khí nén được cho q_out lần lượt 19,038225; 15,900056; 14,390557 cm³/min. Đây là lưu lượng tại áp suất đầu ra từng lần đo, không phải ba lưu lượng ở cùng điều kiện tham chiếu. Bộ số được tạo từ mô hình để có thể kiểm ngược, không phải dữ liệu của một mỏ thực.
Hồi quy đúng trục và đơn vị
Nếu trục x là 1/p_m với p_m tính bằng MPa, ba hoành độ là 5; 2,5; 1,25 MPa⁻¹. Đường thẳng k_app=c+s x đi qua ba điểm với c=20 mD và s=1,6 mD·MPa. Từ đó k_∞=c và b=s/c=0,080 MPa.
Không lấy độ dốc 1,6 làm trực tiếp hệ số b; phải chia cho tung độ gốc. Nếu đổi trục áp suất sang Pa, độ dốc trở thành 1,6×10⁶ mD·Pa nhưng b vẫn mô tả cùng áp suất 80.000 Pa. Bảng tính cần giữ đơn vị trục hồi quy ngay bên cạnh hệ số.
Ngoại suy tới 1/p_m=0 là thao tác toán học. Nó không đòi hỏi tăng áp suất thí nghiệm tới vô hạn. Khoảng đo, độ phân giải và mức phù hợp đường thẳng quyết định mức tin cậy của tung độ gốc.
Hai kiểu hiệu chỉnh không thay thế nhau
Tại điểm đầu, nếu dùng công thức chất lỏng k=μLq_out/(AΔp), kết quả là 29,473684 mD, lớn hơn độ thấm khí biểu kiến đúng 28 mD khoảng 5,263158%. Hiệu chỉnh khí nén được đưa về 28 mD; hiệu chỉnh trượt khí mới đưa tiếp về 20 mD.
Nếu chỉ lấy 28 mD rồi gọi đó là độ thấm nội tại, sai khác còn 40% so với giá trị đã đặt. Ngược lại, áp dụng hệ số trượt hai lần sẽ làm giảm quá mức. Báo cáo cần ghi rõ số nào đã hiệu chỉnh lưu lượng khí nén được và số nào đã ngoại suy Klinkenberg.
Chênh áp không đổi khi chuyển từ áp suất tuyệt đối sang áp suất dư, nhưng p_out và tổng p_in+p_out có đổi. Vì vậy việc dùng Δp đúng chưa bảo đảm công thức bình phương áp suất đúng. Mọi cột áp suất phải ghi rõ mốc trước khi tính.
Độ nhạy của phép đo và phép ngoại suy
Do A tỷ lệ D², chỉ riêng việc nhập D lớn hơn đúng 1% sẽ làm độ thấm suy ra thấp hơn 1,970395%, khi các đầu vào khác giữ nguyên. Tăng sai số chiều dài 1% lại tăng k suy ra 1%; hai sai lệch hình học không có cùng dấu.
Xét bộ số mới chỉ thay điểm áp suất cao nhất từ 22 lên 22,4 mD. Hồi quy bình phương tối thiểu cả ba điểm cho c=20,4 mD, s=1,508571 mD·MPa và b=0,073950 MPa. Một thay đổi 0,4 mD ở một điểm đã làm b giảm khoảng 7,56%, trong khi k_∞ tăng 2%.
Điều này không chứng minh cần loại điểm cuối. Nó cho thấy hai hệ số tương quan và một đường thẳng đi qua ít điểm có thể che độ nhạy. Cần xem số đo lặp, phần dư, độ không đảm bảo áp suất/lưu lượng và độ trải của trục 1/p_m, thay vì chỉ báo một hệ số tương quan cao.
Giữ điều kiện cơ học của mẫu có ý nghĩa
Kansas Geological Survey tách ảnh hưởng ứng suất giam giữ khỏi trượt khí; đóng các cổ lỗ rỗng có thể làm thay đổi độ thấm thực. [2] Khi tăng áp suất lỗ rỗng mà giữ áp suất ngoài áo không đổi, chênh áp nén mẫu thay đổi. Khi ấy xu hướng k theo p_m có thể chứa cả cơ học mẫu lẫn trượt khí.
Bộ số minh họa giả định cấu trúc mẫu và trạng thái cơ học không đổi. Thí nghiệm thực cần nêu cách kiểm soát ứng suất, lịch sử tăng giảm tải và tình trạng mẫu. Không dùng một tương quan hiệu chỉnh riêng của tập đá khác để thay số đo trên loại đá đang xét.
Phạm vi kết luận từ một mẫu lõi
SLB phân biệt độ thấm một pha với độ thấm hiệu dụng và độ thấm tương đối khi nhiều pha cùng có mặt. [3] Giá trị k_∞ của mẫu khô không tự cho lưu lượng dầu hoặc khí trong vỉa có nước, cũng không tự đại diện các khe nứt trên quy mô mỏ.
DTU mô tả việc ghi lưu lượng và áp suất sau khi đạt trạng thái ổn định. [1] Với lõi rất kín, thời gian đạt ổn định, rò nền và dòng vòng qua áo có thể chi phối số đo. Ở tốc độ lớn, giả thiết Darcy tuyến tính cũng cần được kiểm tra bằng các mức chênh áp khác nhau.
Một báo cáo có thể kiểm lại phải giữ kích thước, hướng lấy lõi, trạng thái làm sạch, loại khí, nhiệt độ, áp suất tuyệt đối hai đầu, lưu lượng cùng điều kiện quy chiếu, ứng suất ngoài và dữ liệu từng điểm. Con số độ thấm cuối chỉ có ý nghĩa khi đi cùng các điều kiện đã tạo ra nó.
Nguồn tham khảo
[1] Gas permeability — Technical University of Denmark, Center for Energy Resources Engineering. Truy cập ngày 09/11/2026.
[2] Chapter 4. Reservoir Characterization — Influence of Confining Stress and Klinkenberg Effect — Alan P. Byrnes, Martin K. Dubois and John H. Doveton; Kansas Geological Survey, Open-file Report 2007-06. Truy cập ngày 09/11/2026.
[3] The Defining Series: Defining and Determining Permeability — SLB, Oilfield Review, Richard Nolen-Hoeksema, 9 September 2015. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Nếu hai phòng thí nghiệm dùng cùng mẫu nhưng khác áp suất đo, các số độ thấm khí có thể khác dù cấu trúc mẫu chưa thay đổi. Có ít nhất hai việc phải tách: lưu lượng thể tích thay đổi khi khí giãn nở dọc mẫu, và khả năng dòng khí trượt tại thành lỗ rỗng làm độ thấm biểu kiến phụ thuộc áp suất.
DTU mô tả phép đo ở nhiều áp suất trung bình, sau đó hồi quy độ thấm theo nghịch đảo áp suất để xác định giá trị hiệu chỉnh Klinkenberg. [1] Bài này tự dựng một bộ số giả định để giải thích phép biến đổi, không cung cấp quy trình vận hành thiết bị chịu áp.
Từ Darcy tại một tiết diện đến công thức khí nén được
Xét lõi trụ dài L, tiết diện A, dòng một chiều ổn định, khí lý tưởng đẳng nhiệt, độ nhớt μ không đổi. Bỏ trọng lực, quán tính và rò quanh áo mẫu. Ban đầu giả sử độ thấm k không phụ thuộc vị trí hay áp suất. Darcy cục bộ viết q(x)=−kA(dp/dx)/μ, với q là lưu lượng thể tích tại đúng vị trí x.
Bảo toàn khối lượng cùng nhiệt độ không đổi cho p(x)q(x)=p_out q_out. Thay vào Darcy rồi tích phân từ đầu vào tới đầu ra:
q_out = k A (p_in²−p_out²)/(2 μ L p_out).
k_app = 2 μ L p_out q_out/{A(p_in²−p_out²)}.
p_in và p_out phải là áp suất tuyệt đối, đơn vị Pa; μ dùng Pa·s, L dùng m, A dùng m², q_out dùng m³/s thì k có đơn vị m². Lưu lượng tại đầu ra không được thay bằng lưu lượng quy chuẩn khi chưa quy đổi về p_out và nhiệt độ mẫu.
Đặt p_m=(p_in+p_out)/2 và Δp=p_in−p_out. Vì p_in²−p_out²=2p_mΔp, có thể viết k_app=μLq_out p_out/(A p_mΔp). Dạng này làm rõ hệ số hiệu chỉnh khí nén được p_out/p_m so với phép tính chất lỏng.
Mô hình trượt khí và giá trị ngoại suy
Báo cáo Kansas Geological Survey dùng quan hệ Klinkenberg k_g=k_liquid(1+b/P), đồng thời nêu hệ số b thay đổi theo cấu trúc lỗ rỗng. [2] Với ký hiệu trong bài, k_∞ là giá trị ngoại suy khi nghịch đảo áp suất tiến về không; b có đơn vị áp suất.
k_app = k_∞(1+b/p_m).
Có thể kiểm tính nhất quán bằng giả định cục bộ k(p)=k_∞(1+b/p). Khi tích phân Darcy, tích p k(p)=k_∞(p+b) cho:
p_out q_out = k_∞ A { (p_in²−p_out²)/2 + bΔp }/(μL).
Chia kết quả theo định nghĩa k_app ở trên thu đúng k_∞(1+b/p_m). Vì vậy trong các giả thiết này, áp suất trung bình số học xuất hiện từ tích phân, không phải một lựa chọn tùy ý. Nếu khí thực, nhiệt độ hay độ nhớt thay đổi đáng kể, tích phân cần các hàm tính chất phù hợp.
Dựng ba điểm đo giả định
Cho đường kính mẫu D=0,025 m, L=0,050 m, μ=1,8×10⁻⁵ Pa·s. Diện tích A=πD²/4=0,000490874 m². Chọn k_∞=20 mD và b=0,080 MPa; quy đổi trong ví dụ dùng 1 mD=9,869233×10⁻¹⁶ m².
Giữ Δp=0,020 MPa. Ba lần đo có p_m lần lượt 0,20; 0,40; 0,80 MPa tuyệt đối. Do đó các cặp p_in/p_out là 0,21/0,19; 0,41/0,39; 0,81/0,79 MPa. Theo mô hình, k_app lần lượt 28; 24; 22 mD.
Thay các giá trị vào công thức khí nén được cho q_out lần lượt 19,038225; 15,900056; 14,390557 cm³/min. Đây là lưu lượng tại áp suất đầu ra từng lần đo, không phải ba lưu lượng ở cùng điều kiện tham chiếu. Bộ số được tạo từ mô hình để có thể kiểm ngược, không phải dữ liệu của một mỏ thực.
Hồi quy đúng trục và đơn vị
Nếu trục x là 1/p_m với p_m tính bằng MPa, ba hoành độ là 5; 2,5; 1,25 MPa⁻¹. Đường thẳng k_app=c+s x đi qua ba điểm với c=20 mD và s=1,6 mD·MPa. Từ đó k_∞=c và b=s/c=0,080 MPa.
Không lấy độ dốc 1,6 làm trực tiếp hệ số b; phải chia cho tung độ gốc. Nếu đổi trục áp suất sang Pa, độ dốc trở thành 1,6×10⁶ mD·Pa nhưng b vẫn mô tả cùng áp suất 80.000 Pa. Bảng tính cần giữ đơn vị trục hồi quy ngay bên cạnh hệ số.
Ngoại suy tới 1/p_m=0 là thao tác toán học. Nó không đòi hỏi tăng áp suất thí nghiệm tới vô hạn. Khoảng đo, độ phân giải và mức phù hợp đường thẳng quyết định mức tin cậy của tung độ gốc.
Hai kiểu hiệu chỉnh không thay thế nhau
Tại điểm đầu, nếu dùng công thức chất lỏng k=μLq_out/(AΔp), kết quả là 29,473684 mD, lớn hơn độ thấm khí biểu kiến đúng 28 mD khoảng 5,263158%. Hiệu chỉnh khí nén được đưa về 28 mD; hiệu chỉnh trượt khí mới đưa tiếp về 20 mD.
Nếu chỉ lấy 28 mD rồi gọi đó là độ thấm nội tại, sai khác còn 40% so với giá trị đã đặt. Ngược lại, áp dụng hệ số trượt hai lần sẽ làm giảm quá mức. Báo cáo cần ghi rõ số nào đã hiệu chỉnh lưu lượng khí nén được và số nào đã ngoại suy Klinkenberg.
Chênh áp không đổi khi chuyển từ áp suất tuyệt đối sang áp suất dư, nhưng p_out và tổng p_in+p_out có đổi. Vì vậy việc dùng Δp đúng chưa bảo đảm công thức bình phương áp suất đúng. Mọi cột áp suất phải ghi rõ mốc trước khi tính.
Độ nhạy của phép đo và phép ngoại suy
Do A tỷ lệ D², chỉ riêng việc nhập D lớn hơn đúng 1% sẽ làm độ thấm suy ra thấp hơn 1,970395%, khi các đầu vào khác giữ nguyên. Tăng sai số chiều dài 1% lại tăng k suy ra 1%; hai sai lệch hình học không có cùng dấu.
Xét bộ số mới chỉ thay điểm áp suất cao nhất từ 22 lên 22,4 mD. Hồi quy bình phương tối thiểu cả ba điểm cho c=20,4 mD, s=1,508571 mD·MPa và b=0,073950 MPa. Một thay đổi 0,4 mD ở một điểm đã làm b giảm khoảng 7,56%, trong khi k_∞ tăng 2%.
Điều này không chứng minh cần loại điểm cuối. Nó cho thấy hai hệ số tương quan và một đường thẳng đi qua ít điểm có thể che độ nhạy. Cần xem số đo lặp, phần dư, độ không đảm bảo áp suất/lưu lượng và độ trải của trục 1/p_m, thay vì chỉ báo một hệ số tương quan cao.
Giữ điều kiện cơ học của mẫu có ý nghĩa
Kansas Geological Survey tách ảnh hưởng ứng suất giam giữ khỏi trượt khí; đóng các cổ lỗ rỗng có thể làm thay đổi độ thấm thực. [2] Khi tăng áp suất lỗ rỗng mà giữ áp suất ngoài áo không đổi, chênh áp nén mẫu thay đổi. Khi ấy xu hướng k theo p_m có thể chứa cả cơ học mẫu lẫn trượt khí.
Bộ số minh họa giả định cấu trúc mẫu và trạng thái cơ học không đổi. Thí nghiệm thực cần nêu cách kiểm soát ứng suất, lịch sử tăng giảm tải và tình trạng mẫu. Không dùng một tương quan hiệu chỉnh riêng của tập đá khác để thay số đo trên loại đá đang xét.
Phạm vi kết luận từ một mẫu lõi
SLB phân biệt độ thấm một pha với độ thấm hiệu dụng và độ thấm tương đối khi nhiều pha cùng có mặt. [3] Giá trị k_∞ của mẫu khô không tự cho lưu lượng dầu hoặc khí trong vỉa có nước, cũng không tự đại diện các khe nứt trên quy mô mỏ.
DTU mô tả việc ghi lưu lượng và áp suất sau khi đạt trạng thái ổn định. [1] Với lõi rất kín, thời gian đạt ổn định, rò nền và dòng vòng qua áo có thể chi phối số đo. Ở tốc độ lớn, giả thiết Darcy tuyến tính cũng cần được kiểm tra bằng các mức chênh áp khác nhau.
Một báo cáo có thể kiểm lại phải giữ kích thước, hướng lấy lõi, trạng thái làm sạch, loại khí, nhiệt độ, áp suất tuyệt đối hai đầu, lưu lượng cùng điều kiện quy chiếu, ứng suất ngoài và dữ liệu từng điểm. Con số độ thấm cuối chỉ có ý nghĩa khi đi cùng các điều kiện đã tạo ra nó.
Nguồn tham khảo
[1] Gas permeability — Technical University of Denmark, Center for Energy Resources Engineering. Truy cập ngày 09/11/2026.
[2] Chapter 4. Reservoir Characterization — Influence of Confining Stress and Klinkenberg Effect — Alan P. Byrnes, Martin K. Dubois and John H. Doveton; Kansas Geological Survey, Open-file Report 2007-06. Truy cập ngày 09/11/2026.
[3] The Defining Series: Defining and Determining Permeability — SLB, Oilfield Review, Richard Nolen-Hoeksema, 9 September 2015. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
