Điện trở suất là đầu vào của một mô hình đá
Một đường log có điện trở suất cao không tự xác định được bao nhiêu dầu khí đang có trong lớp đá. Cùng điện trở suất ấy, thay độ rỗng, độ mặn nước hoặc cấu trúc liên thông có thể cho độ bão hòa nước khác nhau. Công thức Archie hữu ích khi các giả thiết phù hợp và các thông số được xác định đúng điều kiện.
Bài gốc của Archie trình bày các quan hệ thực nghiệm gần đúng và nhấn mạnh vai trò của dữ liệu phòng thí nghiệm về đá. [1] Phần dưới tự dựng các phép tính để kiểm đơn vị, độ nhạy và giới hạn diễn giải; các con số không phải kết quả đánh giá một giếng thật.
Ba điện trở suất cần phân biệt
Gọi R_w là điện trở suất của nước vỉa ở điều kiện xét; R₀ là điện trở suất của đá khi phần lỗ rỗng đang xét bão hòa hoàn toàn bởi nước đó; R_t là điện trở suất của đá ở trạng thái bão hòa nước cần tìm. Cả ba dùng Ω·m, không phải Ω.
Điện trở đo qua một mẫu trụ, ký hiệu R_e, dùng Ω. Với dòng một chiều đều qua tiết diện A và chiều dài L trong mô hình hình học lý tưởng, điện trở suất là R_e A/L. Đưa số điện trở chưa hiệu chỉnh hình học vào vị trí R_t sẽ làm kết quả phụ thuộc kích thước mẫu một cách sai lệch.
DTU mô tả phép đo xoay chiều và lưu ý tránh bay hơi làm mẫu khô trong quá trình đo. [2] Tình trạng bão hòa, tiếp xúc điện cực, nhiệt độ và ứng suất cần được giữ hoặc ghi lại. Độ chính xác của đồng hồ điện trở cũng không đồng nghĩa độ chính xác của độ bão hòa suy bằng mô hình đá.
Từ hệ số thành hệ đến độ bão hòa
Hệ số thành hệ F=R₀/R_w và chỉ số điện trở suất I=R_t/R₀ đều không có đơn vị. Trong dạng có hệ số hiệu chuẩn a dùng ở bài này:
F=aφ⁻ᵐ; I=S_w⁻ⁿ; R_t=aR_w/(φᵐS_wⁿ).
S_w={aR_w/(φᵐR_t)}^(1/n).
φ là độ rỗng dưới dạng phần số; S_w là tỷ lệ thể tích nước trong phần lỗ rỗng của mô hình; a, m, n là các hệ số không thứ nguyên. Dạng gốc liên hệ F với độ rỗng và độ tăng điện trở suất khi giảm nước; a được giữ ở đây như một hệ số hiệu chuẩn, không mặc định bằng nhau cho mọi tập đá. [1]
Cơ sở vật lý của cách dùng đơn giản là nước dẫn điện trong một khung đá và các pha còn lại có đóng góp dẫn điện có thể bỏ qua. Glover liên hệ số mũ bão hòa với mức liên thông của các pha dẫn. [3] Vì vậy n mô tả nhiều hơn một hệ số tùy chỉnh để làm kết quả trông hợp lý.
Ví dụ cơ sở và hai mẫu số khác nhau
Tự chọn φ=0,20; R_w=0,080 Ω·m; R_t=20 Ω·m; a=1; m=2; n=2. Khi đó F=25, R₀=2 Ω·m và I=10. Độ bão hòa nước là S_w=√(2/20)=0,316228, tương đương 31,6228% thể tích lỗ rỗng.
Nếu chỉ có nước và hydrocarbon trong lỗ rỗng của mô hình, phần hydrocarbon là 1−S_w=0,683772. Tuy nhiên tỷ lệ so với toàn bộ thể tích đá lại là φ(1−S_w)=0,136754. Không dùng 68,3772% thay cho 13,6754% khi đổi từ thể tích lỗ rỗng sang thể tích khối đá.
Cả hai con số đều chưa là trữ lượng thu hồi. Chúng chưa chứa diện tích, chiều dày hữu hiệu, biến đổi thể tích chất lưu từ vỉa lên mặt, khả năng dịch chuyển hay hiệu suất thu hồi. Phép tính điện trở suất chỉ giải một phần của chuỗi đánh giá.
Lỗi phần trăm có thể tạo kết quả rất thuyết phục
Nếu ô độ rỗng hiển thị 20% nhưng phần tính toán nhận số 20 thay vì 0,20, φ² tăng 10.000 lần. Với n=2, S_w sẽ nhỏ đi 100 lần, còn 0,00316228. Kết quả này vẫn nằm giữa không và một, nên kiểm giới hạn đơn thuần không phát hiện được lỗi.
Phải kiểm dữ liệu nội bộ, không chỉ hình thức ô. Tương tự, R_w và R_t cần cùng đơn vị. Dùng một số theo Ω·cm và số kia theo Ω·m tạo sai lệch tỷ số dù từng con số trông hợp lý trong tài liệu gốc.
Hiệu chuẩn m bằng dữ liệu bão hòa nước hoàn toàn
Xét bộ dữ liệu giả định độc lập: hai mẫu cùng tập đá có φ₁=0,15; φ₂=0,25. Đo được F₁=29,411461 và F₂=11,143047, là các giá trị được tạo từ a=0,8 và m=1,9.
Lấy logarit cho lnF=lna−m lnφ. Từ hai điểm:
m=ln(F₁/F₂)/ln(φ₂/φ₁); a=F₁φ₁ᵐ.
Tính ngược thu m=1,9 và a=0,8. Hai điểm chỉ đủ xác định hai hệ số cho dạng quan hệ đã chọn. Chúng chưa xác nhận mọi mẫu thuộc cùng một quan hệ, nên việc phân tập đá, đo lặp và kiểm phần dư ở các mẫu khác vẫn cần thiết.
Không lấy m từ phép thử thay đổi bão hòa trên một mẫu rồi gọi đó là số mũ độ rỗng. Hai phép thử làm thay đổi hai đại lượng khác nhau; ký hiệu m và n không thể hoán đổi chỉ vì đều xuất hiện dưới dạng lũy thừa.
Hiệu chuẩn n trên cùng mẫu và cùng nước
Một bộ số khác giữ mẫu và nước không đổi. Ở S_w=0,60 và 0,40, chọn I lần lượt 3,076573 và 7,507028, được tạo từ n=2,2. Khi đó n=ln(I₂/I₁)/ln(S_w₁/S_w₂)=2,2.
R₀ phải tương ứng cùng mẫu và cùng nước tham chiếu. Nếu dùng R₀ của mẫu khác hoặc nước có độ dẫn khác, chỉ số I bị thay đổi trước khi hồi quy. Lịch sử giảm hoặc tăng bão hòa và phân bố chất lưu cũng cần được nêu; một số mũ khớp một thí nghiệm không tự bảo đảm áp dụng cho mọi lịch sử của vỉa.
Phân tích độ nhạy bằng logarit
Đặt X=aR_w/(φᵐR_t). Vì lnS_w=lnX/n, các độ nhạy tương đối cục bộ đối với a, R_w, R_t, φ lần lượt là 1/n; 1/n; −1/n; −m/n. Ở trường hợp cơ sở, tăng R_w một tỷ lệ nhỏ làm S_w tăng xấp xỉ một nửa tỷ lệ đó; tăng φ một tỷ lệ nhỏ làm S_w giảm gần cùng tỷ lệ.
Với hai số mũ, đạo hàm theo thay đổi tuyệt đối là ∂lnS_w/∂m=−lnφ/n=0,804719 và ∂lnS_w/∂n=−lnS_w/n=0,575646. Không diễn giải các hệ số này như độ nhạy theo phần trăm của m hoặc n nếu chưa nhân thêm số mũ tương ứng.
Tính đầy đủ cho m=1,8 và 2,2, giữ các đầu vào khác như cơ sở, thu S_w=0,269217 và 0,371447. Nếu thay n bằng 1,8 và 2,2 thì thu 0,278256 và 0,351119. Cùng mức thay đổi 0,2 không tạo cùng tác động cho hai số mũ.
Đổi φ từ 0,20 thành 0,18 hoặc 0,22 cho S_w=0,351364 hoặc 0,287480. Hai thay đổi đối xứng của đầu vào tạo thay đổi đầu ra không hoàn toàn đối xứng. Đây là lý do cần tính lại ở các kịch bản khi khoảng biến động đủ lớn, thay vì kéo dài tuyến tính hóa quá xa.
Một đường dẫn điện bổ sung làm lệch kết quả ra sao
Tự dựng một mô hình phản ví dụ, không phải công thức hiệu chỉnh đá sét. Giữ a=1, m=n=2, φ=0,20, R_w=0,080 Ω·m nhưng đặt S_w thật của mô hình bằng 0,30. Độ dẫn điện phần Archie là σ_A=φ²S_w²/R_w=0,045 S/m.
Giả định một đường dẫn độc lập song song đóng góp thêm σ_x=0,025 S/m. Tổng σ=0,070 S/m, nên R_t=1/σ=14,285714 Ω·m. Nếu vẫn đưa R_t này vào công thức chỉ có nước dẫn, người tính suy ra S_w=0,374166, cao hơn giá trị đã đặt 0,30.
Ví dụ chứng minh chiều sai lệch trong mô hình song song đã chọn, không đưa ra một hệ số sửa chung. Glover thảo luận cách nhìn nhiều pha dẫn và tính liên thông; lượng vật liệu dẫn cùng cách chúng nối nhau đều ảnh hưởng phản ứng điện. [3] Cần một mô hình phù hợp cùng dữ liệu khoáng vật, nước và phép đo để xử lý trường hợp thực.
Kết quả ngoài miền cần được giải thích
Nếu aR_w/(φᵐR_t)>1 thì công thức cho S_w>1. Cắt tự động về một làm mất dấu hiệu mâu thuẫn đầu vào hoặc mô hình. Nên giữ giá trị tính thô và cờ kiểm tra, rồi xem đơn vị, R_w, φ, hệ số đá, ảnh hưởng lớp lân cận và trạng thái log.
Archie đã lưu ý ảnh hưởng hình học giếng, lớp mỏng, xâm nhập và nhiệt độ đối với việc đối chiếu phép đo. [1] Báo cáo cuối nên ghi cả dữ liệu đã hiệu chỉnh, cách chọn thông số, tập đá áp dụng và khoảng độ nhạy. Khi đó người đọc biết con số bão hòa đến từ giả thiết nào và điều gì cần kiểm chứng thêm.
Nguồn tham khảo
[1] The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics — G. E. Archie, Shell Oil Co.; Transactions of the AIME, volume 146, 1942, pp. 54–62; bản lưu tại Pennsylvania State University. Truy cập ngày 09/11/2026.
[2] Electrical resistivity — Technical University of Denmark, Center for Energy Resources Engineering. Truy cập ngày 09/11/2026.
[3] A new theoretical interpretation of Archie's saturation exponent — Paul W. J. Glover, University of Leeds; Solid Earth 8, 805–816, 2017, DOI 10.5194/se-8-805-2017. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một đường log có điện trở suất cao không tự xác định được bao nhiêu dầu khí đang có trong lớp đá. Cùng điện trở suất ấy, thay độ rỗng, độ mặn nước hoặc cấu trúc liên thông có thể cho độ bão hòa nước khác nhau. Công thức Archie hữu ích khi các giả thiết phù hợp và các thông số được xác định đúng điều kiện.
Bài gốc của Archie trình bày các quan hệ thực nghiệm gần đúng và nhấn mạnh vai trò của dữ liệu phòng thí nghiệm về đá. [1] Phần dưới tự dựng các phép tính để kiểm đơn vị, độ nhạy và giới hạn diễn giải; các con số không phải kết quả đánh giá một giếng thật.
Ba điện trở suất cần phân biệt
Gọi R_w là điện trở suất của nước vỉa ở điều kiện xét; R₀ là điện trở suất của đá khi phần lỗ rỗng đang xét bão hòa hoàn toàn bởi nước đó; R_t là điện trở suất của đá ở trạng thái bão hòa nước cần tìm. Cả ba dùng Ω·m, không phải Ω.
Điện trở đo qua một mẫu trụ, ký hiệu R_e, dùng Ω. Với dòng một chiều đều qua tiết diện A và chiều dài L trong mô hình hình học lý tưởng, điện trở suất là R_e A/L. Đưa số điện trở chưa hiệu chỉnh hình học vào vị trí R_t sẽ làm kết quả phụ thuộc kích thước mẫu một cách sai lệch.
DTU mô tả phép đo xoay chiều và lưu ý tránh bay hơi làm mẫu khô trong quá trình đo. [2] Tình trạng bão hòa, tiếp xúc điện cực, nhiệt độ và ứng suất cần được giữ hoặc ghi lại. Độ chính xác của đồng hồ điện trở cũng không đồng nghĩa độ chính xác của độ bão hòa suy bằng mô hình đá.
Từ hệ số thành hệ đến độ bão hòa
Hệ số thành hệ F=R₀/R_w và chỉ số điện trở suất I=R_t/R₀ đều không có đơn vị. Trong dạng có hệ số hiệu chuẩn a dùng ở bài này:
F=aφ⁻ᵐ; I=S_w⁻ⁿ; R_t=aR_w/(φᵐS_wⁿ).
S_w={aR_w/(φᵐR_t)}^(1/n).
φ là độ rỗng dưới dạng phần số; S_w là tỷ lệ thể tích nước trong phần lỗ rỗng của mô hình; a, m, n là các hệ số không thứ nguyên. Dạng gốc liên hệ F với độ rỗng và độ tăng điện trở suất khi giảm nước; a được giữ ở đây như một hệ số hiệu chuẩn, không mặc định bằng nhau cho mọi tập đá. [1]
Cơ sở vật lý của cách dùng đơn giản là nước dẫn điện trong một khung đá và các pha còn lại có đóng góp dẫn điện có thể bỏ qua. Glover liên hệ số mũ bão hòa với mức liên thông của các pha dẫn. [3] Vì vậy n mô tả nhiều hơn một hệ số tùy chỉnh để làm kết quả trông hợp lý.
Ví dụ cơ sở và hai mẫu số khác nhau
Tự chọn φ=0,20; R_w=0,080 Ω·m; R_t=20 Ω·m; a=1; m=2; n=2. Khi đó F=25, R₀=2 Ω·m và I=10. Độ bão hòa nước là S_w=√(2/20)=0,316228, tương đương 31,6228% thể tích lỗ rỗng.
Nếu chỉ có nước và hydrocarbon trong lỗ rỗng của mô hình, phần hydrocarbon là 1−S_w=0,683772. Tuy nhiên tỷ lệ so với toàn bộ thể tích đá lại là φ(1−S_w)=0,136754. Không dùng 68,3772% thay cho 13,6754% khi đổi từ thể tích lỗ rỗng sang thể tích khối đá.
Cả hai con số đều chưa là trữ lượng thu hồi. Chúng chưa chứa diện tích, chiều dày hữu hiệu, biến đổi thể tích chất lưu từ vỉa lên mặt, khả năng dịch chuyển hay hiệu suất thu hồi. Phép tính điện trở suất chỉ giải một phần của chuỗi đánh giá.
Lỗi phần trăm có thể tạo kết quả rất thuyết phục
Nếu ô độ rỗng hiển thị 20% nhưng phần tính toán nhận số 20 thay vì 0,20, φ² tăng 10.000 lần. Với n=2, S_w sẽ nhỏ đi 100 lần, còn 0,00316228. Kết quả này vẫn nằm giữa không và một, nên kiểm giới hạn đơn thuần không phát hiện được lỗi.
Phải kiểm dữ liệu nội bộ, không chỉ hình thức ô. Tương tự, R_w và R_t cần cùng đơn vị. Dùng một số theo Ω·cm và số kia theo Ω·m tạo sai lệch tỷ số dù từng con số trông hợp lý trong tài liệu gốc.
Hiệu chuẩn m bằng dữ liệu bão hòa nước hoàn toàn
Xét bộ dữ liệu giả định độc lập: hai mẫu cùng tập đá có φ₁=0,15; φ₂=0,25. Đo được F₁=29,411461 và F₂=11,143047, là các giá trị được tạo từ a=0,8 và m=1,9.
Lấy logarit cho lnF=lna−m lnφ. Từ hai điểm:
m=ln(F₁/F₂)/ln(φ₂/φ₁); a=F₁φ₁ᵐ.
Tính ngược thu m=1,9 và a=0,8. Hai điểm chỉ đủ xác định hai hệ số cho dạng quan hệ đã chọn. Chúng chưa xác nhận mọi mẫu thuộc cùng một quan hệ, nên việc phân tập đá, đo lặp và kiểm phần dư ở các mẫu khác vẫn cần thiết.
Không lấy m từ phép thử thay đổi bão hòa trên một mẫu rồi gọi đó là số mũ độ rỗng. Hai phép thử làm thay đổi hai đại lượng khác nhau; ký hiệu m và n không thể hoán đổi chỉ vì đều xuất hiện dưới dạng lũy thừa.
Hiệu chuẩn n trên cùng mẫu và cùng nước
Một bộ số khác giữ mẫu và nước không đổi. Ở S_w=0,60 và 0,40, chọn I lần lượt 3,076573 và 7,507028, được tạo từ n=2,2. Khi đó n=ln(I₂/I₁)/ln(S_w₁/S_w₂)=2,2.
R₀ phải tương ứng cùng mẫu và cùng nước tham chiếu. Nếu dùng R₀ của mẫu khác hoặc nước có độ dẫn khác, chỉ số I bị thay đổi trước khi hồi quy. Lịch sử giảm hoặc tăng bão hòa và phân bố chất lưu cũng cần được nêu; một số mũ khớp một thí nghiệm không tự bảo đảm áp dụng cho mọi lịch sử của vỉa.
Phân tích độ nhạy bằng logarit
Đặt X=aR_w/(φᵐR_t). Vì lnS_w=lnX/n, các độ nhạy tương đối cục bộ đối với a, R_w, R_t, φ lần lượt là 1/n; 1/n; −1/n; −m/n. Ở trường hợp cơ sở, tăng R_w một tỷ lệ nhỏ làm S_w tăng xấp xỉ một nửa tỷ lệ đó; tăng φ một tỷ lệ nhỏ làm S_w giảm gần cùng tỷ lệ.
Với hai số mũ, đạo hàm theo thay đổi tuyệt đối là ∂lnS_w/∂m=−lnφ/n=0,804719 và ∂lnS_w/∂n=−lnS_w/n=0,575646. Không diễn giải các hệ số này như độ nhạy theo phần trăm của m hoặc n nếu chưa nhân thêm số mũ tương ứng.
Tính đầy đủ cho m=1,8 và 2,2, giữ các đầu vào khác như cơ sở, thu S_w=0,269217 và 0,371447. Nếu thay n bằng 1,8 và 2,2 thì thu 0,278256 và 0,351119. Cùng mức thay đổi 0,2 không tạo cùng tác động cho hai số mũ.
Đổi φ từ 0,20 thành 0,18 hoặc 0,22 cho S_w=0,351364 hoặc 0,287480. Hai thay đổi đối xứng của đầu vào tạo thay đổi đầu ra không hoàn toàn đối xứng. Đây là lý do cần tính lại ở các kịch bản khi khoảng biến động đủ lớn, thay vì kéo dài tuyến tính hóa quá xa.
Một đường dẫn điện bổ sung làm lệch kết quả ra sao
Tự dựng một mô hình phản ví dụ, không phải công thức hiệu chỉnh đá sét. Giữ a=1, m=n=2, φ=0,20, R_w=0,080 Ω·m nhưng đặt S_w thật của mô hình bằng 0,30. Độ dẫn điện phần Archie là σ_A=φ²S_w²/R_w=0,045 S/m.
Giả định một đường dẫn độc lập song song đóng góp thêm σ_x=0,025 S/m. Tổng σ=0,070 S/m, nên R_t=1/σ=14,285714 Ω·m. Nếu vẫn đưa R_t này vào công thức chỉ có nước dẫn, người tính suy ra S_w=0,374166, cao hơn giá trị đã đặt 0,30.
Ví dụ chứng minh chiều sai lệch trong mô hình song song đã chọn, không đưa ra một hệ số sửa chung. Glover thảo luận cách nhìn nhiều pha dẫn và tính liên thông; lượng vật liệu dẫn cùng cách chúng nối nhau đều ảnh hưởng phản ứng điện. [3] Cần một mô hình phù hợp cùng dữ liệu khoáng vật, nước và phép đo để xử lý trường hợp thực.
Kết quả ngoài miền cần được giải thích
Nếu aR_w/(φᵐR_t)>1 thì công thức cho S_w>1. Cắt tự động về một làm mất dấu hiệu mâu thuẫn đầu vào hoặc mô hình. Nên giữ giá trị tính thô và cờ kiểm tra, rồi xem đơn vị, R_w, φ, hệ số đá, ảnh hưởng lớp lân cận và trạng thái log.
Archie đã lưu ý ảnh hưởng hình học giếng, lớp mỏng, xâm nhập và nhiệt độ đối với việc đối chiếu phép đo. [1] Báo cáo cuối nên ghi cả dữ liệu đã hiệu chỉnh, cách chọn thông số, tập đá áp dụng và khoảng độ nhạy. Khi đó người đọc biết con số bão hòa đến từ giả thiết nào và điều gì cần kiểm chứng thêm.
Nguồn tham khảo
[1] The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics — G. E. Archie, Shell Oil Co.; Transactions of the AIME, volume 146, 1942, pp. 54–62; bản lưu tại Pennsylvania State University. Truy cập ngày 09/11/2026.
[2] Electrical resistivity — Technical University of Denmark, Center for Energy Resources Engineering. Truy cập ngày 09/11/2026.
[3] A new theoretical interpretation of Archie's saturation exponent — Paul W. J. Glover, University of Leeds; Solid Earth 8, 805–816, 2017, DOI 10.5194/se-8-805-2017. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
