- Xác ñịnh thông số vỉa dầu theo lịch sử khai thác trên cơ sở Phương pháp Gauss-Newton sửa ñổi. - Bài báo trình bày việc áp dụng phuơng pháp giải lặp Gauss-Newton sửa ñổi cho hệ phương trình phi tuyến ñể xác ñịnh thông số mô hình mô phỏng của các vỉa dầu theo lịch sử khai thác (phục hồi lịch sử). - Các thông số vỉa như ñộ rỗng, ñộ thấm hay bất kỳ thông số mô hình mô phỏng vỉa dầu cần phải xác ñịnh ñều có thể coi là các tham biến ñể cực tiểu hoá hàm tổng bình phương sai số và bất kỳ số liệu thực tế nào cũng có thể sử dụng ñể ñưa vào ñánh giá bình phương sai số. - Một số kỹ thuật cũng ñã ñược sử dụng ñể ñảm bảo tính hội tụ và giảm thời gian tính khi bài toán lớn và phi tuyến.. - Trong các tính toán, chương trình mô phỏng vỉa là chương trình IMEX của CMG, Canada.. - Các kết quả cũng ñã ñược so sánh, ñánh giá và cho thấy quá trình hội tụ là tương ñối tốt và ñảm bảo ñộ chính xác. - Xác ñịnh các thông số cho các mô hình mô phỏng là bài toán cần thiết và có ý nghĩa quan trọng không những trong công nghệ mô phỏng các vỉa dầu mà còn có thể áp dụng cho nhiều vấn ñề khác – vấn ñề xác ñịnh thông số cho các mô hình mô phỏng nói chung.. - Việc xác ñịnh thông số mô hình vỉa theo lịch sử khai thác là công ñoạn quan trọng ñể xây dựng mô hình số phù hợp mô phỏng hoạt ñộng, tính toán khai thác hiệu quả các vỉa dầu.. - Công việc này yêu cầu một quá trình thay ñổi, tính thử theo các thông số vỉa giả ñịnh và ñánh giá sai số giữa kết quả tính toán bằng mô hình và số liệu quan sát thực tế. - Quá trình này ñược thực hiện liên tục cho ñến khi nhận ñược sự phù hợp cần thiết giữa kết quả tính và số liệu khai thác thực tế và hiện nay chủ yếu làm bằng thủ công, tốn nhiều thời gian, ñộ tin cậy chưa cao vì vậy gây không ít khó khăn trong việc ñề xuất kế hoạch khai thác các vỉa dầu, ñặc biệt là các viả có cấu trúc phức tạp, các vỉa lớn, hoặc mỏ có nhi ều vỉa phân tán [1,2].. - Nhiều phương pháp xác ñịnh thông số mô hình vỉa bằng các chương trình - Tự ñộng phục hồi lịch sử ñã ñược phát triển [3,4,5]. - Tuy nhiên một số phương pháp ñòi hỏi nhiều kỹ thuật phức tạp, và thường yêu cầu phát triển riêng chương trình tính toán mô phỏng vỉa phù hợp (hoặc phải có chương trình nguồn ñể thay ñổi). - Áp dụng phương pháp lặp Gauss-Newton cho hệ phương trình phi tuyến có ưu ñiểm là ñơn giản, cho phép sử dụng các phần mềm mô phỏng vỉa có sẵn (như là một hàm ẩn), vấn ñề là ñiều khiển ñược các file dữ liệu vào và file kết quả ra của chương trình mô phỏng [8]. - Trên cơ sở phương pháp và thuật toán ñã phát triển, chương trình tính toán xác ñịnh một số loại thông số vỉa cơ bản như ñộ thấm, ñộ rỗng...cũng ñã ñược tạo lập. - Chương trình có thể sử dụng trực tiếp khi dùng chương trình mô phỏng vỉa RESSIM (Viện Cơ học) hoặc IMEX (CMG, Canada, ñang sử dụng rộng rãi ở các Công ty dầu khí trên Thế giới và Việt nam) hoặc với những sửa ñổi nhỏ khi dùng với một phần mềm mô phỏng vỉa khác. - Chương trình ñã sử dụng ñể xác ñịnh thông số cho một số mô hình vỉa ví dụ. - Các kết quả cho thấy tất các bài toán ñều hội tụ và cho kết quả nghiệm với ñộ chính xác khá cao và hoàn toàn có thể ñưa vào. - áp dụng cho tính toán thực tế vỉa hiện có của Việt Nam.. - Xác ñịnh thông số mô hình vỉa dầu-khí qua số liệu khai thác. - Mô hình mô phỏng vỉa dầu-khí. - Trong tính toán khai thác các vỉa dầu khí hay nước ngầm, công cụ hữu hiệu ñã ñược phát triển và sử dụng rộng rãi ñể tính toán các phương án khai thác là các chương trình, phần mềm tính toán mô phỏng vỉ a. - Trong các chương trình, phần mềm này, mô hình sử dụng hiệu quả nhất hiện nay là mô hình dòng chảy ba pha (dầu-khí-nước. - Trong ñó các thông số của mô hình gồm các thông số ñặc trưng cho tính chất hình học và vật lý của môi trường cũng như của các chất lỏng như: ñộ rỗng ( φ. - Các tham biến của mô hình là áp suất (P), ñộ bão hoà của các pha (S).. - Các phương trình (1)-(3) cùng với các biểu thức ñóng kín ñược giải bằng phương pháp số ñể thu nhận các phân bố theo thời gian của các tham biến trạng thái chính của vỉa như áp suất, ñộ bão hoà pha…trong quá trình khai thác.. - Trong bài báo này chương trình mô phỏng sử dụng là chương trình imex (cmg-canada), phương pháp giải là phương pháp sai phân hữu hạn. - Mặc dù vẫn tồn tại những khó khăn nhất ñịnh, các chương trình mô phỏng vỉa ñã ñược phát triển và sử dụng rất thành công trong công tác tính toán các phương án khai thác các vỉa dầu khí, nước ngầm trên thế giới cũng như ở Việt Nam. - Tuy nhiên ñộ chính xác của các kết quả tính bằng mô hình phụ thuộc vào tính chính xác của các thông số mô hình vỉa và hiện nay, khó khăn lớn nhất trong kỹ thuật vỉa vẫn là xác ñịnh các thông số mô hình, công việc tốn rất nhiều thời gian, công sức.. - Xác ñịnh thông số vỉa qua số liệu khai thác Khi một mô hình mô phỏng vỉa ñã ñược xây dựng, bài toán tính toán dự báo các phương án khai thác là có thể ñược giải quyết, vấn ñề còn lại trong quản lý khai thác là chọn phương án trên cơ sở các tính toán dự báo. - Bài toán quản lý này có thể giải quyết trên cơ sở ghép nối bài toán mô phỏng vỉa với bài toán ñiều khiển tối ưu.. - Khó khăn ñầu tiên khi áp dụng phương pháp sử dụng công cụ mô phỏng vào thực tế là các thông số cho mô hình mô phỏng phải ñược xác ñịnh phù hợp, sao cho các kết quả mô phỏng mô tả ñúng các quan sát thực tế. - Một mô hình mô phỏng vỉa ñúng phải thoả mãn hai ñiều kiện sau:. - Công cụ mô phỏng phải cho phép lập mô hình mô tả ñúng với thực tế quan sát vỉa.. - Các thông số của mô hình phải ñược xác ñịnh ñúng với các thông số vỉa thực tế.. - Trong thực tế, việc xác ñịnh chính xác các thông số cho các mô hình cho toàn bộ vỉa là hết sức khó khăn. - có thể xác ñịnh một cách dễ dàng tại các giếng (khai thác hoặc bơm ép) trong quá trình khai thác vỉa (lịch sử khai thác). - Từ ñó, với yêu cầu một mô hình ñúng phải cho kết quả ñúng với các giá trị của các biến trạng thái quan sát khi các thông số vỉa ñược chọn ñúng và ngược lại, chúng ta có thể sử dụng các giá trị quan sát của các biến trạng thái ñể xác ñịnh các thông số vỉa thông qua các mô hình mô phỏng vỉa (phục hồi lịch sử).. - Phương pháp lặp Newton cho hệ phương trình phi tuyến áp dụng cho bài toán xác ñịnh thông số vỉa. - Áp dụng phương pháp ñ iều khiển tối ưu cho bài toán xác ñịnh thông số vỉa. - Trong thực tế mô phỏng vỉa, các tham biến như áp suất tại giếng, tỷ số dầu-nước, tỷ số khí- dầu thường ñược sử dụng ñể xác ñịnh, hiệu chỉnh các thông số mô hình vỉa. - Trong quá trình xác ñịnh, hiệu chỉnh các thông số vỉa, một số thông số và tính chất chất lỏng coi như ñã biết, một số thông số khác cần phải xác ñịnh. - Tập các thông số cần phải xác ñịnh, hiệu chỉnh ký hiệu là X(X 1 , X 2 ,...X N. - Mục ñích ñặt ra trong bài báo này là xác ñịnh tập các thông số X sao cho sai số bình phương giữa kết quả tính toán từ mô hình và các giá trị quan sát thực tế của các tham biến là bé nhất. - P k obs là giá trị áp suất quan sát, P k cal là giá trị áp suất tính toán bằng các chương trình mô phỏng vỉa, m là số ñiểm quan sát.. - Các thông số cần xác ñịnh X(X 1 ,X 2 ,...X N ) thông thường là ñộ rỗng, ñộ thấm của các vùng khác nhau trong vỉa, mỗi vùng có thể một hay nhiều ô lưới. - Dựa trên các hiểu biết về vỉa thông qua các kỹ thuật vỉa, chúng ta có thể có ñược các ràng buộc cho miền giá trị của các thông số. - Như vậy, bài toán xác ñịnh thông số vỉa sử dụng các chương trình mô phỏng vỉa ñược phát biẻu dưới dạng bài toán tối ưu như sau: Xác ñịnh tập các giá trị X(X 1 ,X 2 ,..X N ) với các ràng buộc (5) ñể hàm sai số (4), E(X 1 ,X 2 ,..X N ) ñạt giá trị cực tiểu.. - i=1,2,...N (8) ðây là hệ phương trình dạng cổ ñiển của phương pháp Gauss-Newton cho bài toán bình phương tối thiểu. - Trong ñó, các ñạo hàm riêng của hàm ẩn F(X) tính ñược bằng cách tính toán mô phỏng với một thay ñổi của một thông số X i. - trong khi giữ nguyên các thông số khác. - ðối với bài toán xác ñịnh thông số vỉa, phương trình (8) là phi tuyến vì vậy các phương pháp lặp phải ñược áp dụng và luôn gặp phải vấn ñề hội tụ. - Tuy nhiên khi áp dụng cho các bài toán thực tế, các kỹ thuật này cũng không mấy hiệu quả vì yêu cầu chạy nhiều mô phỏng xuôi, quá trình lặp hay bị dừng ở cực tiểu ñịa phương hay các ñiểm yên ngựa.. - Trong chương trình tính của bài báo này chúng tôi ñã áp dụng, phát triển và kiểm tra một kỹ thuật ñơn giản và hiệu quả hơn, tóm tắt như sau:. - ðể thoát khỏi các ñiểm dừng hoặc cực tiểu ñịa phương, cách gây nhiễu ñược thực hiện bằng cách tăng lần lượt số gia của các thông số khi tính ñạo hàm riêng số của hàm mục tiêu trong (10).. - Trên cơ sở phương pháp và các kỹ thuật ñã trình bày ở trên, chương trình tính ñã ñược viết bằng ngôn ngữ FORTRAN. - Chương trình sau ñó ñã chạy thử, kiểm tra và tính toán cho một số ví dụ áp dụng.. - Sơ ñồ khối chương trình xác ñịnh thông số vỉa.. - Giả sử bộ thông số X o =(X o 1 ,X o 2 ,...X o n. - Gọi chương trình mô phỏng - Tính hàm sai số E o =E(X o. - Gọi chương trình mô phỏng - Tính hàm sai số E=E(X). - Các ví dụ tính toán áp dụng. - Xác ñịnh thông số mô hình vỉa của bài toán SPE1. - Trong bài toán này, mô hình vỉa ñề xuất bởi Aziz, S.O. - [14] và ñược SPE (Society of Petroleum Engineers) sử dụng như một trong những bài toán mẫu (SPE1) ñể kiểm tra và chuẩn hoá các chương trình mô phỏng vỉa trước khi công bố và ñưa vào áp dụng thực tế. - Bài toán này thuộc loại bài toán dòng chảy 3 pha, 3 chiều, bơm ép khí 5 ñiểm. - Quá trình hoạt ñộng của vỉa ñã ñược mô phỏng ñầy ñủ bằng IMEX với mô hình vỉa gồm 10x10x3 ô lưới. - Trong vỉa có 2 giếng, 1 giếng khai thác ñặt tại ô lưới giếng bơm ép ñặt tại ô lưới (1,1,1) (lớp trên cùng là 1). - Thời gian mô phỏng là 10 năm.. - Kết quả mô phỏng với tập thông số gốc của bài toán này ñược giữ lại một phần như là số liệu quan sát (59 giá trị áp suất tại ô giếng khai thác theo thời gian). - Sau ñó các thông số vỉa ñược thay ñổi ñể chương trình xác ñịnh thông số tìm lại các thông số sao cho giá trị tính toán mô phỏng phù hợp với số liệu quan sát.. - Các thông số cần xác ñịnh là ñộ thấm theo hướng I của các lớp (3 thông số, giá trị ñộ thấm là bằng nhau trong mỗi lớp). - Chương trình ñã sử dụng cho hai trường hợp là dòng chảy 1 pha và hai pha (tương ứng là có và không có giếng bơm ép). - Trên Bảng 1 là kết quả của quá trình hội tụ và trên Hình 2, Hình 3 và Hình 4 là kết quả phục hồi lịch sử sủa áp suất tại ô lưới chứa giếng khai thác, lưu lượng dầu và tỷ số khí dầu (GOR) cho bài toán một pha. - Tương tự Bảng 2 và Hình 5, Hình 6 và Hình 7 là kết quả tính toán cho bài toán SPE1 hai pha (có giếng bơm ép).. - Kết quả cho thấy nghiệm tìm ñược ở bài toán này gần như chính xác vì số thông số ít, vỉa khá ñồng nhất và kích thước vỉa không lớn. - Hội tụ (1P). - Phục hồi áp suất ô giếng khai thác (SPE1, 1 pha).. - Hội tụ (2P).. - Phục hồi áp suất ô giếng khai thác (SPE1, 2 pha).. - Lưu lượng dầu khai thác (2 pha).. - Xác ñịnh thông số cho mô hình vỉa của bài toán SPE9. - Bài toán kiểm tra mô hình vỉa thứ 9 của hiệp hội dầu khí (SPE9) [15] là bài toán 3 pha, 3 chiều, vỉa có kích thước 2200x2300x110m và có tính bất ñồng nhất cao (ñộ thấm và ñộ rỗng thay ñổi mạnh theo không gian). - Vỉa ñược khai thác bởi 25 giếng khai thác và 1 giếng bơm ép, thời gian hoạt ñộng là 2.5 năm. - Trong bài toán này, mô hình vỉa chia thành 24x25x15 ô lưới có ñộ thấm khác nhau, ñộ rỗng của các ô khác nhau theo lớp. - ðể kiểm tra chương trình xác ñịnh thông số mô hình vỉa cho bài toán này, cách làm cũng ñược tiến hành tương tự như bài. - Kết quả mô phỏng của bài toán SPE9 là áp suất tại các ô giếng khai thác cũng ñược giữ lại như là số liệu quan sát về áp suất.. - Sau ñó thông số vỉa là ñộ rỗng cũng bị thay ñổi ñể chương trình tự xác ñịnh các giá trị ñộ thấm cho từng lớp (15 lớp, lớp trên cùng là 1) sao cho kết quả mô phỏng phù hợp với kết quả quan sát. - Quá trình hội tụ của bài toán SPE9 Số lần. - Thời gian tính của các bài toán. - ô lưới thông số Sô lần lặp Sai số Thời gian. - Bài toán SPE1-1P phút. - Bài toán SPE1-2P phút. - Bài toán SPE phút. - Trên cơ sở sai số bình phương tối thiểu giữa số liệu quan sát và kết quả tính toán từ mô hình, thuật toán tối ưu ñã ñược áp dụng ñể xây dựng và giải bài toán xác ñịnh, hiệu chỉnh thông số mô hình cho các mô hình mô phỏng các vỉa dầu-khí theo số liệu quan sát từ thực tế khai thác. - Một chương trình tính ñã ñược lập và ñã tiến hành tính toán với một số vỉa ví dụ. - Các kết quả tính kiểm tra cho thấy quá trình hội tụ là tương ñối tốt. - Chương trình hoàn toàn có khả năng phát triển, áp dụng cho thực tế khai thác các vỉa dầu-khí cũng như các vỉa nước ngầm hiện ñang khai thác của Việt nam.. - [1] Phạm Thu Thuỷ, Xây dựng mô hình khai thác cho tầng móng nứt nẻ, ví dụ áp dụng cho mỏ Ruby, bể Cửu long, Tuyển tập báo cáo Hội nghị KHCN “30 năm Dầu khí Việt nam: Cơ hội mới, thách thức mới”, Hà nội, 7/2005, quyển 1, tr. - [3] Trương Công Tài, Nguyễn Văn Út, Một số vấn ñề về công nghệ thiết kế khai thác và kết quả khai thác tầng móng mỏ Bạch Hổ, Tuyển tập báo cáo Hội nghị KHCN “30 năm Dầu khí Việt nam:. - [7] ðặng Thế Ba, Áp dụng lý thuyết ñiều khiển tối ưu cho bài toán xác ñịnh thông số mô phỏng của các mô hình mô phỏng vỉa dầu nhiều pha, Tuyển tập Báo cáo hội nghị khoa học Thuỷ khí Toàn quốc, ðà nẵng, 7/2003, trang 1.