BÁO CÁO CÁ NHÂN “BÁO CÁO XÂY DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT BẰNG CÁC PHƯƠNG PHÁP VIỄN THÁM, ĐỊA CHẤT, ĐỊA VẬT LÝ”

NHIỆM VỤ THƯỜNG XUYÊN THEO CHỨC NĂNGNĂM 2018TÊN NHIỆM VỤ: NGHIÊN CỨU X Y DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA, NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆTMà SỐ: VĐCKS.2018.02 (032018PGV NVTXTCN) BÁO CÁO CÁ NH N“BÁO CÁO X Y DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT BẰNG CÁC PHƯƠNG PHÁP VIỄN THÁM, ĐỊA CHẤT, ĐỊA VẬT LÝ”MỤC LỤCMỤC LỤC I. TỔNG QUAN 11. Tổng quan về địa nhiệt 12. Nguồn Năng Lượng Địa Nhiệt 1II. CÁC PHƯƠNG PHÁP ĐIỀU TRA ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT 71. Tổng quan các phương pháp điều tra đánh giá nguồn địa nhiệt 72. Phương pháp viễn thám 73 Các phương pháp Địa chất 94. Phuơng pháp Địa vật lý 13III. SƠ LƯỢC HIỆN TRẠNG ĐIỀU TRA ĐÁNH GIÁ TIỀM NĂNG ĐỊA NHIỆT VIỆT NAM.1 24IV. VÍ DỤ THỰC TẾ Đà THỰC HIỆN ĐỂ ĐIỀU TRA ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT TẠI VIỆT NAM. 271. Các phương pháp địa vật lý đã thi công và khối lượng 282. Kỹ thuật thi công công tác địa vật lý 282.1. Mạng lưới khảo sát 282.2. Kỹ thuật thi công các phương pháp địa vật lý 283. Công tác xử lý và tổng hợp tài liệu 303.1. Tách các trường địa phương và khu vực 303.2 Phân tích định lượng tài liệu điện theo mô hình hai chiều (2D) 304. Kết quả công tác địa vật lý 334.1. Phương pháp từ 334.2. Kết quả của phương pháp mặt cắt phân cực 334.3 Phương pháp đo sâu phân cực 334.4 Kết quả tính lượng sinh nhiệt do phóng xạ 34KẾT LUẬN 39TÀI LIỆU DẪN 40I. TỔNG QUAN1. Tổng quan về địa nhiệt Địa Nhiệt là nguồn nhiệt năng có sẵn trong lòng đất. Cụ thể hơn, nguồn năng lượng nhiệt này tập trung ở khoảng vài km dưới bề mặt Trái Đất, phần trên cùng của vỏ Trái Đất. Cùng với sự tăng nhiệt độ khi đi sâu vào vỏ Trái Đất, nguồn nhiệt lượng liên tục từ lòng đất này được ước đoán tương đương với với một khoảng năng lượng cỡ 42 triệu MW. Lòng đất thì vẫn tiếp tục nóng hằng tỷ năm nữa, đảm bảo một nguồn nhiệt năng gần như vô tận. Chính vì vậy Địa Nhiệt được liệt vào dạng năng lượng tái tạo. Nguồn nhiệt lượng này được chuyển lên mặt đất qua dạng hơi hoặc nước nóng khi nước chảy qua đất đá nóng. Nhiệt lượng thường được sử dụng trực tiếp, ví dụ như hệ thống điều hòa nhiệt độ (bơm địa nhiệt), hoặc chuyển thành điện năng (nhà máy nhiệt điện). Địa nhiệt là dạng năng lượng sạch và bền vững. So với các dạng năng lượng tái tạo khác như gió, thủy điện hay điện mặt trời, địa nhiệt không phụ thuộc vào các yếu tố thời tiết và khí hậu. Do đó địa nhiệt cũng có hệ số công suất rất cao, nguồn địa nhiệt luôn sẵn sàng 24hngày,7 ngày trong tuần.Cho đến nay, hơn 30 quốc gia trên thế giới đã khai thác tổng cộng 12.000 MW địa nhiệt cho các ứng dụng trực tiếp và sản xuất hơn 8.000 MW điện. Tại một vài quốc gia đang phát triển, địa nhiệt điện chiếm một vai trò đáng kể trong việc đáp ứng nhu cầu điện.Các nhà máy địa nhiệt có giới hạn công suất từ 100 kW cho đến 100 MW, phụ thuộc vào nguồn năng lượng vào nhu cầu điện năng. Kỹ thuật này rất thích hợp cho điện khí hóa nông thôn và các ứng dụng mạng lưới mini (minigrid), bên cạnh ứng dụng trong việc hòa mạng quốc gia. Tại các quốc gia có nguồn tài nguyên eo hẹp hoặc có điều kiện khí hậu khắc nghiệt, địa nhiệt điện có thể đóng một vai trò rất hữu dụng. Các ứng dụng trực tiếp của địa nhiệt có thể góp phần tăng đáng kể sản lượng nông nghiệp và ngư nghiệp (nuôi trồng thủy hải sản) và cung cấp nhiệt cho các quá trình xử lý công nghiệp phụ trợ. Nguồn địa nhiệt được xem là đặc biệt quan trọng đối với các quốc gia đang phát triển mà lại không có các nguồn tài nguyên năng lượng như than, dầu và khí tự nhiên.2. Nguồn Năng Lượng Địa Nhiệt Nguồn địa nhiệt liên quan mật thiết đến cấu trúc nhiệt độ của Trái Đất và chu trình đối lưu nhiệt trong lòng Trái Đất. Nhiệt độ của Trái Đất tăng dần theo độ sâu và đạt đến 4.200oC tại tâm (Hình1.a). Hình 1. a) Cấu trúc nhiệt độ sơ giản và b) các dạng biểu hiện của nguồn địa nhiệt thường quan sát được trên mặt đất như lỗkhenúi lửa (fumarole), suối nước nóng(hot spring), bọt bùn (mud pot) và một số dạng khác (Nguồn: Geothermal education Ofice) Hình 2. a) Vỏ trái đất không phải là một khối vật chất liền mạch mà bị phân chia ra thành các mảng lớn. Do động lực của đối lưu nhiệt phía dưới, các mảng này có thể dịch chuyển ra xa nhau (đới phân kỳ) hoặc hút vào nhau (đới hội tụ), hoặc trượt ngang qua nhau ở một vận tốc rất chậm (vài cmnăm).b) tại đới phân kỳ ở giữa các đại dương (spreading center) và tại nft lục địa , vỏ Trái đất được “tái sinh” . Tại các đới hội tụ các mảnh có thể trượt chồm lên nhau (subduction). ở rìa của mảng trượt chìm, magma có thể dâng lên về phía bề mặt Trái đất. Tại những nơi vỏ Trái đất “yếu”, magma có thể phụ trào lên mặt đất trở thành lava. Tuy nhiên đa phần magma không phun trào lên đến mặt đất mà chỉ thường đun nóng một diện tích khổng lồ các lớp đất đá bên dưới mặt đất, đồng thời đun nóng các dòng nước ngầm (bắt nguồn từ nước mưa thấm qua bề mặt đất hoặc len theo các đứt gãy bề mặt xuống sâu trong lòng đất). Nước ngầm bị nung nóng có thể dâng lên mặt đất dưới dạng hơi nước hoặc dưới dạng nước nóng. Đây chính là nguồn địa nhiệt mà con người đang khai thác.Một phần trong tổng khối nhiệt lượng khổng lồ trong lòng Trái Đất này bắt nguồn từ quá trình hình thành hành tinh trong khoảng 4,5 tỷ năm trước (Trái Đất hình thành từ một khối cầu vật chất cực nóng, nguội dần từ trong ra ngoài qua quá trình quay quanh trục), và phần còn lại là kết quả của quá trình phân rã của các nguyên tố phóng xạ tồn tại trong lõi Trái Đất. Theo nguyên lý tuần hoàn nhiệt lượng từ nơi nhiệt độ cao xuống nhiệt độ thấp, dòng nhiệt của Trái Đất di chuyển từ trong lõi ra ngoài vỏ. Dưới tác động của một quá trình địa chất gọi là kiến tạo mảng (xem tài liệu VnGG), vỏ Trái Đất được phân ra thành 12 mảng lớn (mảng kiến tạo, Hình 2.a) và được tái tạo (tái sinh) một cách chậm chạp qua hàng triệu năm. Các mảng này di chuyển tương đối với nhau (phân tách hoặc hội tụ) với tốc độ vài cmnăm. Khi hai mảng kiến tạo va chạm vào nhau, 1 mảng có thể hút chìm xuống mảng còn lại, tạo nên các trũng đại dương và gây ra động đất (Hình 2.b). Đây chính là nơi vỏ Trái Đất trở nên yếu hơn bình thường, cho phép vật chất nóng từ trong lòng đất dịch chuyển lên mặt. Ở độ sâu lớn tại đới hội tụ, ngay bên dưới mảng sụp chìm, nhiệt độ tăng lên đủ cao đến nung chảy đất đá và tạo ra magma (nham thạch). Do có mật độ thấp hơn khối đất đá xung quanh, magma di chuyển lên phía trên vỏ Trái Đất và mang theo nhiệt lượng cùng với nó. Đôi khi magma di chuyển lên tới bề mặt Trái Đất thông qua các điểm yếu của vỏ Trái Đất và phun trào lava tại các miệng núi lửa. Tuy nhiên, đa phần magma được giữ lại trong vỏ Trái Đất và nung nóng đất đá và các khối nước ngầm (subterranean water). Một phần khối nước nóng này có thể di chuyển lên mặt đất thông qua các đới đứt gãy hoặc khe đá rạn (cracks), hình thành suối nước nóng (hay là geysers, mạch nước nóng). Một khi khối nước nóng và hơi nước này bị “bẫy” do khối đất đá không thấm (impermeable) ở phía bên trên và được giữ lại trong khối đất đá thấm (permeable), bồn trũng địa nhiệt được hình thành (Hình 3). Các bồn trũng này chính là nguồn địa nhiệt có thể được dùng trực tiếp hoặc để sản xuất điện qua hệ thống turbine hơi nước (steam turbine). Hình 3. a) Bể địa nhiệt được hình thành khi nước nóng và hơi nước dâng lên bề mặt và bị “bẫy“ bởi tầng đá “nóc” không thấm và được giữa lại bởi đá rỗng; b) mô hình một hệ địa nhiệt lý tưởng. Một hệ địa nhiệt bao gồm 3 thành phần: nguồn nhiệt (heat source), bể địa nhiệt (reservoir) và lưu chất địa nhiệt (heat trasntsfer fluid). Nguồn nhiệt có thể là magma intrution có nhiệt độ cao (>6000) dâng lên phần trên cùng của vỏ trái đất (510km), hoặc chỉ đơn giản là nhiệt độ chuẩn của trái đất tăng theo độ sâu trong một số trường hợp hệ địa nhiệt nhiệt độ thấp. Bể địa nhiệt là một khối đất thấm (permeable) cho phép chất lỏng hồi chuyển nhiệt có thể chảy qua thu nhiệt từ nguồn. Bể địa nhiệt thường được “phủ” bởi đá không thấm ((impermeable rock) và thông với một khu vực “tái sinh” trên bề mặt thông qua đó lưu chất địa nhiệt có thể được lưu thông và tái sử dụng qua hệ thống các bơm lỗ khoan. Chất lỏng địa nhiệt có thể là nước ngầm đa phần là meteoric water ở pha lỏng hoặc khí tùy thuộc vào nhiệt độ và áp suất bồn. lưu chất này thường mang các hợp chất hóa học hoặc khí như CO2, H2S… (Nguồn: Geothermal education Ofice IGA) Có 5 dạng nguồn địa nhiệt khác nhau, trong đó chỉ có bồn trũng thủy địa nhiệt (hydrothermal reservoirs) và năng lượng trái đất (earth energy) là đã được đưa vào khai thác thương mại. 3 dạng còn lại, nước muối địa áp (geopressureed brine), đá khô nóng (dry hot rock) và magma, vẫn còn yêu cầu phát triển các kỹ thuật caotân tiến. Bể thủy nhiệt là các bể chứa hơi hoặc nước nóng bị bẫy trong đá porous (Hình 3). Để sản xuất điện, hơi hoặc nước nóng được bơm từ các bể lên mặt đất để vận hành các turbin phát điện. Do nguồn hơi nước tương đối hiếm, nên hầu hết các nhà máy địa nhiệt sử dụng nguồn nước nóng. Chi tiết về kỹ thuật sẽ được giới thiệu ở phần sau. Đá khô nóng: địa nhiệt có thể được khai thác từ một số các nguồn đá khô, không thấm ở độ sâu khoảng 510 m dưới mặt đất, hoặc thậm chí nông hơn ở một số khu vực. Ý tưởng chủ đạo là bơm nước lạnh xuống nguồn đá khô này tại một giếng khoan, cho khối nước này chảy qua nguồn đá khô và được nung nóng, sau đó dẫn khối nước được nung nóng ra một giếng khoan khác và trữ trong bể địa nhiệt. Tuy nhiên hiện nay vẫn chưa có ứng dụng thương mại nào cho kỹ thuật này (xem Hình 4). Hình 4. Sơ đồ công nghệ khai thácđịa nhiệt dạng khô nóng (hot dry rock. HDR). Dự án HDR đầu tiên được thử nghiệm vào năm 1970 tại Los Alamos, Neww Mexico, Hoa Kỳ. Thông qua các giếng khoan đặc biệt, nước áp suất cao được bơm vào khu đá nhiệt độ cao, nén cặt nhằm gây ra hiện tượng “rạn” thủy lực (hydraulic fracturing). Nước sẽ thấm vào các khe rạn “nhân tạo”này, từ đó thu nhiệt của lớp địa nhietj cao ở xung quanh, tương tự như dạng bồn địa nhiệt tự nhiên . Bồn địa nhiệt này sau đó sẽ được khoan đến qu một giếng thứ haiddeer thu nhiệt của lưu chất được đun nóng. Do đó, hệ HDR bao gồm (i) lỗ khoan dùng để tạo “rạn” thủy lực qua đó được bơm vào (ii) bồn địa nhiệt nhân tạo và (iii) giếng khoanđể thu lưu chất nhiệt độ cao. Toàn bộ hệ thống này có thể tạo thành một chu kỳ hoàn toàn khép kín(garnish, 987) ( Nguồn: IGA và Deep Heat Mining project (DHM), Thụy Sỹ)Magma: tất cả các kỹ thuật địa nhiệt hiện nay đều chỉ khai thác “gián tiếp” nhiệt năng từ lòng đất do magma chuyển lên. Hiện tại vẫn chưa có kỹ thuật này cho phép khai thác trực tiếp nhiệt lượng từ magma, mặc dù magma là nguồn nhiệt lượng cực kỳ dồi dào trong vỏ Trái Đất. Nước muối địa áp là dạng nước nóng, áp suất cao và chứa methane hòa tan. Cả nhiệt và methane đều có thể được sử dụng để sản xuất điện thông qua turbine.II. CÁC PHƯƠNG PHÁP ĐIỀU TRA ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT 1. Tổng quan các phương pháp điều tra đánh giá nguồn địa nhiệtMặc dù nguồn địa nhiệt không thể quan sát trực tiếp trên bề mặt đất nhưng lại có rất nhiều phương pháp có thể cung cấp các thông tin về sự tồn tại của nó trong những điều kiện nhất định, thậm chí có thể cung cấp cả những thông tin về chất lượng của nguồn địa nhiệt từ những vị trí trên mặt đất hay từ không gian. Các phương pháp điều tra, khảo sát trên mặt nguồn địa nhiệt rất ít khi thành công bởi các kết quả thu được thường đưa ra một bức tranh không hoàn chỉnh về địa chất địa nhiệt. Tuy nhiên, các phương pháp này lại thường kinh tế hơn nhiều so với các phương pháp khảo sát bên dưới mặt đất. Các phương pháp địa chất đưa ra những nhận định quan trọng ban đầu cho mọi phương pháp khảo sát nguồn địa nhiệt khác, bao gồm phân tích các tài liệu địa chất và khảo sát ngoài thực địa. Phương pháp viễn thám, nghĩa là phương pháp điều tra từ máy bay hay và vệ tinh đã đang trở thành một công cụ hữu hiệu cho nghiên cứu nguồn địa nhiệt. Cuối cùng, các phương pháp địa vật lý mà đặc biệt là phương pháp điện và phương pháp địa chấn đã cung cấp cho ta những thông tin gián tiếp về nguồn địa nhiệt. Do vậy, các thông số địa chất địa nhiệt phải được suy luận, giải đoán ra từ các tài liệu khảo sát trên bề mặt. Tuy nhiên, trong khi nghiên cứu cần phải kết họp với tài liệu khảo sát dưới mặt đất đế đánh giá độ họp lý của các tài liệu thu được từ các phương pháp bề mặt.2. Phương pháp viễn thámCác hình ảnh về trái đất thu được từ máy bay hay vệ tinh ở những bước sóng điện từ khác nhau có thể cung cấp cho ta những thông tin hữu ích liên quan đến điều kiện tồn tại của nguồn địa nhiệt. Công nghệ viễn thám đã và đang phát triển nhanh chóng những năm gần đây trong khi việc ứng dụng của nó trong nghiên cứu tài nguyên địa nhiệt vẫn còn đang được nghiên cứu và khám phá. Hơn nữa, việc sẵn có các bức ảnh này ở các hãng tư nhân và các cơ quan đã làm thúc đẩy việc ứng dụng chúng.Nghiên cứu các bức ảnh máy bay đen trắng ba chiều có thể thu được những thông tin quan trọng. Các cấu trúc địa chất, màu sắc và địa hình quan sát được có thể phân biệt được sự khác nhau về điều kiện địa chất, các loại đất, độ ẩm của đất, thảm thực vật và hiện trạng sử dụng đất. Do đó các bức ảnh về địa chất có thể phân biệt được giữa các loại đá với các loại đất cho ta biết được tính thấm và diện phân bố của chúng, từ đó xác định được diện tích của các nguồn cấp và nguồn thoát của nguồn địa nhiệt. Các bản đồ phân vùng khả năng cung cấp nguồn địa nhiệt thành khu vực cung cấp địa nhiệt tốt, trung bình và kém có thể được thành lập. Bảng 1 tóm tắt vai trò của ảnh viễn thám trong việc trợ giúp công tác minh giải điều kiện địa chất – địa nhiệt.Các bức ảnh từ không gian cũng có thể chỉ ra các khe nứt. Các khe nứt này có thể liên quan đến độ rỗng, tính thấm và cuối cùng là lưu lượng của lỗ khoan. Vị trí các con suối và các khu vực đầm lầy chỉ ra rằng ở những nơi đó mực nguồn địa nhiệt nằm tương đối nông. Việc nghiên cứu về thuỷ thực vật từ ảnh viễn thám có thể rất hữu ích. Các thực vật nhạy cảm địa nhiệt sẽ giúp ích cho việc xác định độ sâu của mực nguồn địa nhiệt do chúng có thể hút địa nhiệt từ mực nguồn địa nhiệt ở độ sâu không lớn. Ngoài ra, nghiên cứu các phổ điện từ cũng giữ một vai trò quan trọng trong hệ thống phương pháp ảnh ứng dụng trong khảo sát và nghiên cứu địa chất – địa nhiệt. Phương pháp ảnh hồng ngoại là phương pháp ghi lại sự khác biệt về nhiệt độ tại các vị trí khác nhau trên bề mặt trái đất. Phương pháp này có thể cho ta thông tin về độ ẩm của đất, sự vận động của nguồn địa nhiệt, và các đứt gãy đóng vai trò như là những tầng dẫn địa nhiệt. Một trong những kết quả thú vị nhất của phương pháp ảnh hồng ngoại là đã thành lập được bản đồ của các dòng chảy địa nhiệt nóng và địa nhiệt lạnh dưới đất ở khu vực ven biển, ở những diện tích phân bố đá bazan hoặc đá vôi. Hình 2 cho thấy những dòng thoát của nguồn địa nhiệt có nhiệt độ thấp hơn địa nhiệt biển đã lộ ra xung quanh vùng đảo Hawaii. Ảnh rada có thể cung cấp cho ta thông tin về độ ẩm trên bề mặt hoặc ở những vị trí nông dưới mặt đất. Cuối cùng, những khảo sát điện từ tần thấp đã cho ta thấy vị trí các dòng mặt bị chôn vùi . Địa hình + Đánh giá điều kiện địa hình toàn vùng + Đánh giá điều kiện địa hình khu vực Thực vật không ưa địa nhiệt Các thành tạo địa chất thấm dẫn nhiệt từ dưới lên Các hồ, suối Các hồ suối có hiện tượng bốc hơi nóng hoặc nước sôi. Các điểm địa nhiệt nhân tạo + Lỗ khoan + Giếng đào + Các bồn chứa + Các kênh đào....3 Các phương pháp Địa chấtCác nghiên cứu địa chất có khả năng đánh giá tiềm năng nguồn địa nhiệt cho một vùng rộng lớn một cách nhanh chóng và kinh tế. Một khảo sát địa chất được bắt đầu bằng việc thu thập, phân tích và minh giải về điều kiện địa chất – địa nhiệt từ các bản đồ địa hình, ảnh viễn thám, bản đồ địa chất, tài liệu địa vật lý lỗ khoan và các tài liệu liên quan khác. Ngoài các tài liệu thu thập được, nếu có thể cũng cần tiến hành khảo sát thực địa để đánh giá mức độ tin cậy của các tài liệu thu thập được thông qua việc nghiên cứu các hiện tượng địa nhiệt trên mặt, lỗ khoan, nguồn cấp và thoát của nguồn địa nhiệt, đặc điểm địa nhiệt.Những hiểu biết về sự hoạt động của kiến tạo trẻ cho biết được quy mô của các thành tạo nguồn địa nhiệt. Thông tin về các loại đá sẽ cho ta thông tin về mức độ chứa địa nhiệt có thể phù hợp cho việc cung cấp địa nhiệt cho hộ gia đình chứ không chỉ phù hợp cho quy mô công nghiệp hoặc cho một thành phố. Địa tầng và lịch sử địa chất của một vùng có thể cho ta biết về các tầng chứa địa nhiệt ở bên dưới, tính liên tục và mối liên hệ giữa các tầng chứa địa nhiệt với nhau. Thành phần vật chất, chiều dày các tầng chắn bên trên, thế nằm của các thành tạo chứa địa nhiệt sẽ giúp ta tính toán được độ sâu lỗ khoan. Tương tự, từ dấu hiệu các tầng chứa địa nhiệt có thể phán đoán được chiều sâu đặt máy khai thác. Các đứt gãy đã hình thành nên các đới dập vỡ là miềm thoát của nguồn địa nhiệt từ dưới sâu.3. 1.Vị trí của lãnh thổ so với những vành đai động của hành tinh1Những trường địa nhiệt (geothermal field) lớn của Trái Đất chủ yếu nằm trùng với các vành đai động (mobile belt) của hành tinh, thường là ranh giới giữa các mảng kiến tạo, dọc theo đó hay xảy ra sự va chạm, hút chìm, trượt ngang, tách giãn, gây nên những trận động đất, phun trào núi lửa và những biểu hiện địa động lực khác, như sự xuất lộ những hố phun bùn sôi, lỗ phun khí núi lửa (fumarol, mofet, solfata), geyser và những mạch nước nóng,... Những hiện tượng trên cũng thường thấy ở các sống núi giữa đại dương (midocean ridge), đới rift, các điểm nóng (hot spot) dưới biển và trên đất liền.Nhìn lên bản đồ kiến tạo toàn cầu ta thấy rõ lãnh thổ Việt Nam ở vào một vị trí khá đặc biệt: nằm “lọt thỏm” giữa đôi “gọng kìm” của hai vành đai động: trước mặt là nhánh tây của “Vòng lửa quanh Thái Bình Dương” (Circum Pacific Ring of Fire), sau lưng và dưới chân là nhánh đông của “đai AlpesHimalaya” (Alpine Himalayan Belt hay Alpides) (hình 1), nhưng (may thay) khoảng cách từ dải đất hình chữ S của chúng ta đến cái “vòng kim cô oan nghiệt” đó luôn duy trì ở mức vừa đủ để tránh được những thảm họa khủng khiếp như đã từng xảy ra ở các nước láng g

BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG VIỆN KHOA HỌC ĐỊA CHẤT VÀ KHOÁNG SẢN NHIỆM VỤ THƯỜNG XUYÊN THEO CHỨC NĂNG NĂM 2018 TÊN NHIỆM VỤ: MÃ SỐ: NGHIÊN CỨU XÂY DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA, NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT VĐCKS.2018.02 (03/2018/PGV- NVTXTCN) BÁO CÁO CÁ NHÂN “BÁO CÁO XÂY DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT BẰNG CÁC PHƯƠNG PHÁP VIỄN THÁM, ĐỊA CHẤT, ĐỊA VẬT LÝ” HÀ NỘI, 11-2018 BỘ TÀI NGUYÊN VÀ MÔI TRƯỜNG 1 VIỆN KHOA HỌC ĐỊA CHẤT VÀ KHOÁNG SẢN NHIỆM VỤ THƯỜNG XUYÊN THEO CHỨC NĂNG NĂM 2018 TÊN NHIỆM VỤ: MÃ SỐ: NGHIÊN CỨU XÂY DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA, NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT VĐCKS.2018.02 (03/2018/PGV- NVTXTCN) BÁO CÁO CÁ NHÂN “BÁO CÁO XÂY DỰNG TỔ HỢP CÁC PHƯƠNG PHÁP ĐIỀU TRA NGHIÊN CỨU ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT BẰNG CÁC PHƯƠNG PHÁP VIỄN THÁM, ĐỊA CHẤT, ĐỊA VẬT LÝ” CHỦ TRÌ THỰC HIỆN NGƯỜI THỰC HIỆN Nguyễn Hồng Quang Bùi Xuân Khải CƠ QUAN CHỦ TRÌ VIỆN KHOA HỌC ĐỊA CHẤT VÀ KHOÁNG SẢN HÀ NỘI, 11-2018 2 MỤC LỤC 3 I TỔNG QUAN Tổng quan địa nhiệt Địa Nhiệt nguồn nhiệt có sẵn lòng đất Cụ thể hơn, nguồn lượng nhiệt tập trung khoảng vài km bề mặt Trái Đất, phần vỏ Trái Đất Cùng với tăng nhiệt độ sâu vào vỏ Trái Đất, nguồn nhiệt lượng liên tục từ lòng đất ước đoán tương đương với với khoảng lượng cỡ 42 triệu MW Lịng đất tiếp tục nóng tỷ năm nữa, đảm bảo nguồn nhiệt gần vơ tận Chính Địa Nhiệt liệt vào dạng lượng tái tạo Nguồn nhiệt lượng chuyển lên mặt đất qua dạng nước nóng nước chảy qua đất đá nóng Nhiệt lượng thường sử dụng trực tiếp, ví dụ hệ thống điều hịa nhiệt độ (bơm địa nhiệt), chuyển thành điện (nhà máy nhiệt điện) Địa nhiệt dạng lượng bền vững So với dạng lượng tái tạo khác gió, thủy điện hay điện mặt trời, địa nhiệt không phụ thuộc vào yếu tố thời tiết khí hậu Do địa nhiệt có hệ số công suất cao, nguồn địa nhiệt sẵn sàng 24h/ngày,7 ngày tuần Cho đến nay, 30 quốc gia giới khai thác tổng cộng 12.000 MW địa nhiệt cho ứng dụng trực tiếp sản xuất 8.000 MW điện Tại vài quốc gia phát triển, địa nhiệt điện chiếm vai trò đáng kể việc đáp ứng nhu cầu điện Các nhà máy địa nhiệt có giới hạn cơng suất từ 100 kW 100 MW, phụ thuộc vào nguồn lượng vào nhu cầu điện Kỹ thuật thích hợp cho điện khí hóa nơng thơn ứng dụng mạng lưới mini (mini-grid), bên cạnh ứng dụng việc hòa mạng quốc gia Tại quốc gia có nguồn tài nguyên eo hẹp có điều kiện khí hậu khắc nghiệt, địa nhiệt điện đóng vai trị hữu dụng Các ứng dụng trực tiếp địa nhiệt góp phần tăng đáng kể sản lượng nông nghiệp ngư nghiệp (nuôi trồng thủy hải sản) cung cấp nhiệt cho q trình xử lý cơng nghiệp phụ trợ Nguồn địa nhiệt xem đặc biệt quan trọng quốc gia phát triển mà lại nguồn tài nguyên lượng than, dầu khí tự nhiên Nguồn Năng Lượng Địa Nhiệt Nguồn địa nhiệt liên quan mật thiết đến cấu trúc nhiệt độ Trái Đất chu trình đối lưu nhiệt lòng Trái Đất Nhiệt độ Trái Đất tăng dần theo độ sâu đạt đến 4.200oC tâm (Hình1.a) 4 Hình a) Cấu trúc nhiệt độ sơ giản b) dạng biểu nguồn địa nhiệt thường quan sát mặt đất lỗ/khenúi lửa (fumarole), suối nước nóng(hot spring), bọt bùn (mud pot) số dạng khác (Nguồn: Geothermal education Ofice) 5 Hình a) Vỏ trái đất khơng phải khối vật chất liền mạch mà bị phân chia thành mảng lớn Do động lực đối lưu nhiệt phía dưới, mảng dịch chuyển xa (đới phân kỳ) hút vào (đới hội tụ), trượt ngang qua vận tốc chậm (vài cm/năm).b) đới phân kỳ đại dương (spreading center) nft lục địa , vỏ Trái đất “tái sinh” Tại đới hội tụ mảnh trượt chồm lên (subduction) rìa mảng trượt chìm, magma dâng lên phía bề mặt Trái đất Tại nơi vỏ Trái đất “yếu”, magma phụ trào lên mặt đất trở thành lava Tuy nhiên đa phần magma không phun trào lên đến mặt đất mà thường đun nóng diện tích khổng lồ lớp đất đá bên mặt đất, đồng thời đun nóng dịng nước ngầm (bắt nguồn từ nước mưa thấm qua bề mặt đất len theo đứt gãy bề mặt xuống sâu lịng đất) Nước ngầm bị nung nóng dâng lên mặt đất dạng nước dạng nước nóng Đây nguồn địa nhiệt mà người khai thác 6 Một phần tổng khối nhiệt lượng khổng lồ lòng Trái Đất bắt nguồn từ trình hình thành hành tinh khoảng 4,5 tỷ năm trước (Trái Đất hình thành từ khối cầu vật chất cực nóng, nguội dần từ ngồi qua q trình quay quanh trục), phần cịn lại kết q trình phân rã nguyên tố phóng xạ tồn lõi Trái Đất Theo nguyên lý tuần hoàn nhiệt lượng từ nơi nhiệt độ cao xuống nhiệt độ thấp, dòng nhiệt Trái Đất di chuyển từ lõi ngồi vỏ Dưới tác động q trình địa chất gọi kiến tạo mảng (xem tài liệu VnGG), vỏ Trái Đất phân thành 12 mảng lớn (mảng kiến tạo, Hình 2.a) tái tạo (tái sinh) cách chậm chạp qua hàng triệu năm Các mảng di chuyển tương (phân tách hội tụ) với tốc độ vài cm/năm Khi hai mảng kiến tạo va chạm vào nhau, mảng hút chìm xuống mảng cịn lại, tạo nên trũng đại dương gây động đất (Hình 2.b) Đây nơi vỏ Trái Đất trở nên yếu bình thường, cho phép vật chất nóng từ lòng đất dịch chuyển lên mặt Ở độ sâu lớn đới hội tụ, bên mảng sụp chìm, nhiệt độ tăng lên đủ cao đến nung chảy đất đá tạo magma (nham thạch) Do có mật độ thấp khối đất đá xung quanh, magma di chuyển lên phía vỏ Trái Đất mang theo nhiệt lượng với Đơi magma di chuyển lên tới bề mặt Trái Đất thông qua điểm yếu vỏ Trái Đất phun trào lava miệng núi lửa Tuy nhiên, đa phần magma giữ lại vỏ Trái Đất nung nóng đất đá khối nước ngầm (subterranean water) Một phần khối nước nóng di chuyển lên mặt đất thông qua đới đứt gãy khe đá rạn (cracks), hình thành suối nước nóng (hay geysers, mạch nước nóng) Một khối nước nóng nước bị “bẫy” khối đất đá không thấm (impermeable) phía bên giữ lại khối đất đá thấm (permeable), bồn trũng địa nhiệt hình thành (Hình 3) Các bồn trũng nguồn địa nhiệt dùng trực tiếp để sản xuất điện qua hệ thống turbine nước (steam turbine) 7 Hình a) Bể địa nhiệt hình thành nước nóng nước dâng lên bề mặt bị “bẫy“ tầng đá “nóc” khơng thấm lại đá rỗng; b) mơ hình hệ địa nhiệt lý tưởng Một hệ địa nhiệt bao gồm thành phần: nguồn nhiệt (heat source), bể địa nhiệt (reservoir) lưu chất địa nhiệt (heat trasntsfer fluid) Nguồn nhiệt magma intrution có nhiệt độ cao (>6000) dâng lên phần vỏ trái đất (5-10km), đơn giản nhiệt độ chuẩn trái đất tăng theo độ sâu số trường hợp hệ địa nhiệt nhiệt độ thấp Bể địa nhiệt khối đất thấm (permeable) cho phép chất lỏng hồi chuyển nhiệt chảy qua thu nhiệt từ nguồn Bể địa nhiệt thường “phủ” đá không thấm ((impermeable rock) thông với khu vực “tái sinh” bề mặt thông qua lưu chất địa nhiệt lưu thông tái sử dụng qua hệ thống bơm lỗ khoan Chất lỏng địa nhiệt nước ngầm đa phần meteoric water pha lỏng khí tùy thuộc vào nhiệt độ áp suất bồn lưu chất thường mang hợp chất hóa học khí CO 2, H2S… (Nguồn: Geothermal education Ofice &IGA) 8 Có dạng nguồn địa nhiệt khác nhau, có bồn trũng thủy địa nhiệt (hydrothermal reservoirs) lượng trái đất (earth energy) đưa vào khai thác thương mại dạng cịn lại, nước muối địa áp (geopressureed brine), đá khơ nóng (dry hot rock) magma, cịn u cầu phát triển kỹ thuật cao/tân tiến • Bể thủy nhiệt bể chứa nước nóng bị bẫy đá porous (Hình 3) Để sản xuất điện, nước nóng bơm từ bể lên mặt đất để vận hành turbin phát điện Do nguồn nước tương đối hiếm, nên hầu hết nhà máy địa nhiệt sử dụng nguồn nước nóng Chi tiết kỹ thuật giới thiệu phần sau • Đá khơ nóng: địa nhiệt khai thác từ số nguồn đá khô, không thấm độ sâu khoảng 5-10 m mặt đất, chí nơng số khu vực Ý tưởng chủ đạo bơm nước lạnh xuống nguồn đá khô giếng khoan, cho khối nước chảy qua nguồn đá khơ nung nóng, sau dẫn khối nước nung nóng giếng khoan khác trữ bể địa nhiệt Tuy nhiên chưa có ứng dụng thương mại cho kỹ thuật (xem Hình 4) Hình Sơ đồ cơng nghệ khai thácđịa nhiệt dạng khơ nóng (hot dry rock HDR) Dự án HDR thử nghiệm vào năm 1970 Los Alamos, Neww Mexico, Hoa Kỳ Thông qua giếng khoan đặc biệt, nước áp suất cao bơm vào khu đá nhiệt độ cao, nén cặt nhằm gây tượng “rạn” thủy lực (hydraulic 9 fracturing) Nước thấm vào khe rạn “nhân tạo”này, từ thu nhiệt lớp địa nhietj cao xung quanh, tương tự dạng bồn địa nhiệt tự nhiên Bồn địa nhiệt sau khoan đến qu giếng thứ haiddeer thu nhiệt lưu chất đun nóng Do đó, hệ HDR bao gồm (i) lỗ khoan dùng để tạo “rạn” thủy lực qua bơm vào (ii) bồn địa nhiệt nhân tạo (iii) giếng khoanđể thu lưu chất nhiệt độ cao Toàn hệ thống tạo thành chu kỳ hồn tồn khép kín(garnish, 987) ( Nguồn: IGA Deep Heat Mining project (DHM), Thụy Sỹ) • • Magma: tất kỹ thuật địa nhiệt khai thác “gián tiếp” nhiệt từ lòng đất magma chuyển lên Hiện chưa có kỹ thuật cho phép khai thác trực tiếp nhiệt lượng từ magma, magma nguồn nhiệt lượng dồi vỏ Trái Đất Nước muối địa áp dạng nước nóng, áp suất cao chứa methane hịa tan Cả nhiệt methane sử dụng để sản xuất điện thông qua turbine II CÁC PHƯƠNG PHÁP ĐIỀU TRA ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT Tổng quan phương pháp điều tra đánh giá nguồn địa nhiệt Mặc dù nguồn địa nhiệt quan sát trực tiếp bề mặt đất lại có nhiều phương pháp cung cấp thơng tin tồn điều kiện định, chí cung cấp thơng tin chất lượng nguồn địa nhiệt từ vị trí mặt đất hay từ khơng gian Các phương pháp điều tra, khảo sát mặt nguồn địa nhiệt thành cơng kết thu thường đưa tranh khơng hồn chỉnh địa chất - địa nhiệt Tuy nhiên, phương pháp lại thường kinh tế nhiều so với phương pháp khảo sát bên mặt đất Các phương pháp địa chất đưa nhận định quan trọng ban đầu cho phương pháp khảo sát nguồn địa nhiệt khác, bao gồm phân tích tài liệu địa chất khảo sát thực địa Phương pháp viễn thám, nghĩa phương pháp điều tra từ máy bay hay vệ tinh trở thành công cụ hữu hiệu cho nghiên cứu nguồn địa nhiệt Cuối cùng, phương pháp địa vật lý mà đặc biệt phương pháp điện phương pháp địa chấn cung cấp cho ta thông tin gián tiếp nguồn địa nhiệt Do vậy, thông số địa chất - địa nhiệt phải suy luận, giải đoán từ tài liệu khảo sát bề mặt Tuy nhiên, nghiên cứu cần phải kết họp với tài liệu khảo sát mặt đất đế đánh giá độ họp lý tài liệu thu từ phương pháp bề mặt Phương pháp viễn thám Các hình ảnh trái đất thu từ máy bay hay vệ tinh bước 10 10 khoan dầu khí có chiều sâu 3-4 nghìn mét không nghiên cứu đầy đủ địa nhiệt) - Việc nghiên cứu chất lượng địa nhiệt nóng chủ yếu nhằm vào nhiệt độ thành phần ion Chưa có tài liệu phân tích khống vật đá chứa tích tụ travertin, đá tufa vơi (sinter), nên chưa đủ giải thích nguồn gốc, điều kiện hình thành nhiệt dịch - Việc tính tốn trữ địa nhiệt thực số điểm lộ riêng lẻ sử dụng phương pháp đơn giản, nên số liệu kết mức dự báo, độ tin cậy thấp Vì vậy, để có hiểu biết đầy đủ, toàn diện tài nguyên địa nhiệt đất địa nhiệt , thời gian tới cần tiếp tục điều tra nghiên cứu chi tiết hơn, đồng thời tính tốn kỹ mặt kinh tế Dẫu sao, dựa lý thuyết quy luật thành tạo, phân bố trường nhiệt Trái Đất kết nghiên cứu bước đầu thu được, khẳng định đánh giá tổng quan tiềm chất lượng địa nhiệt lãnh thổ Việt Nam nêu có sở khoa học thực tế, làm định hướng nghiên cứu, đề xuất giải pháp khai thác sử dụng địa nhiệt hợp lý, góp phần phát triển lượng địa nhiệt ta Khu vực Trung du miền núi phía Bắc Việt Nam với đặc trưng địa chất kiến tạo đại nhiều cơng trình nghiên cứu đề cập đến nhiên đại đa số khơng đề cập tới tiềm khống sản nhiên liệu tiềm địa nhiệt khu vực [4] Gần số cơng trình nghiên cứu tác giả công bố báo cáo từ dự án điều tra địa nhiệt khống nóng, kết nghiên cứu địa nhiệt chưa đề cập sâu [5, 6, 7], Các nghiên cứu văn liệu quốc tế bước đầu cho thấy nguồn địa nhiệt khu vực trung du miền núi phía Bắc có tiềm cho ứng dụng khai thác lượng Theo Muraokal, (2008) [8] với nhiệt độ nguồn địa nhiệt > 53oC sử dụng cơng nghệ Kalina để phát điện Trên sở đó, chúng tơi lựa chọn nguồn địa nhiệt xuất lộ có nhiệt độ > 53 oC để tính tốn lượng tự nhiên mà khơng khai thác nhiệt lượng nguồn địa nhiệt tỏa khơng khí xung quanh lãng phí Trên thực tế, thơng thường nguồn địa nhiệt nguồn có nhiệt độ bồn chứa cao (> 100oC), đủ phát điện nhờ công nghệ Chu kỳ Nhị nguyên (Binary Cycle), chúng tơi sử dụng phương pháp ước tính cơng suất phát điện Muffler, P Cataldi, R (1978) [9] để tính cho nguồn địa nhiệt 29 29 Trong tổng số 164 nguồn địa nhiệt nóng xuất lộ khu vực nghiên cứu có tới 18 nguồn có nhiệt độ mặt >53oC, thiếu số thông số cần thiết để tính tốn lượng ước tính cơng suất phát điện nên 18 nguồn lựa chọn để tính tốn, 18 nguồn địa nhiệt nằm địa phận tỉnh Sơn La, Điện Biên, Lai Châu, Hà Giang, Tuyên Quang, Yên Bái, Lào Cai Nghệ An (Hình 1) Hình 14 Sơ đồ vị trí nguồn địa nhiệt triển vọng cấu trúc địa chất gắn với cấu trúc địa nhiệt tiềm đơn giản hóa từ đồ Địa chất Khống sản Việt Nam tỷ lệ 1/200.000 thuộc vùng trung du miền núi phía Bắc [20] Các số vịng trọn biểu thị nguồn: 1- Pe Luông, 2-Na Hai, 3- Pom Lót/Uva, 4-Pa Thơm, 5-Pa Bát, 6-Pác Ma, 7-La Si, 8-Sin Chải, 9Nậm Cải, 10-Làng Sang, 11-Nậm Păm, 12-Lũng Pơ, 13-Bó Đướt, 14-Quảng Ngần, 15-Quảng Ngun, 16-Mỹ Lâm, 17-Nam Ron, 18-Kim Đa IV VÍ DỤ THỰC TẾ ĐÃ THỰC HIỆN ĐỂ ĐIỀU TRA ĐÁNH GIÁ NGUỒN ĐỊA NHIỆT TẠI VIỆT NAM Kết công tác địa vật lý: Nghiên cứu đánh giá tiềm địa nhiệt 30 30 vùng Pom Lót thuộc đề tài “Nghiên cứu đánh giá tiềm địa nhiệt vùng Tây Bắc Bắc Bộ triển vọng sử dụng chúng phục vụ phát triển kinh tế - xã hội” Viện Khoa học Địa chất khoáng sản, 2002 Các phương pháp địa vật lý thi công khối lượng Căn vào khối lượng kinh phí cơng tác địa vật lý cấp năm 2002 đề tài “Nghiên cứu đánh giá tiềm địa nhiệt vùng Tây Bắc Bắc Bộ triển vọng sử dụng chúng phục vụ phát triển kinh tế - xã hội”, tiến hành thi công tổ hợp phương pháp địa vật lý khối lượng vùng Pom Lót cụ thể sau: ST T Đo từ mặt đất Khối lượng kế hoạch 3000 điểm Khối lượng thực 3091 điểm Đo mặt cắt phân cực 800 điểm 1092 điểm Đo sâu phân cực 20 điểm 31 điểm Phương pháp Kỹ thuật thi công công tác địa vật lý 2.1 Mạng lưới khảo sát Tại điểm nước nóng Pom Lót chúng tơi tiến hành đo phương pháp địa vật lý 14 tuyến, từ tuyến T17 đến tuyến T16 Khoảng cách tuyến đo 100 - 200m khoảng cách điểm đo tuyến 10m Phương vị tuyến đo 80 Chiều dài trung bình tuyến khảo sát 2000m 2.2 Kỹ thuật thi công phương pháp địa vật lý 2.2.3 Phương pháp từ Phương pháp từ thực tất tuyến với khoảng cách điểm đo tuyến 10m Máy đo từ kế proton ENVIMAG Canada, có độ xác cao với sai số ± 1nT, máy đo trường từ toàn phần T Các điểm đo tuyến đo lần tính giá trị trung bình Việc hiệu chỉnh biến thiên thực cách hàng ngày đo biến thiên vị trí nơi đóng quân, phút đọc giá trị Việc đo kiểm tra tiến hành đo độc lập tuyến T13, T3, T0 T6 với tỷ lệ 10 %, kết tính sai số 1,3nT, sai só nhỏ sai số cho phép 31 31 2.2.2 Phng phỏp mt ct phõn cc kớch thớch Phơng pháp mặt cắt phân cực kích thích xoay chiều (đo mặt cắt phân cực miền tần số) với tần số 0.5Hz đợc tiến hành 12 đoạn tuyến (từ tuyến T17 đến tuyến T8), khoảng cách điểm đo tun lµ 10m - Máy đo: Chúng tơi sử dụng trạm phân cực xoay chiều (đầu thu V5 đầu phát T-3) Canađa thiết kế dựa tiến kỹ thuật vi xử lý Các số liệu đo đạc, tính tốn đổ trực tiếp vào máy tính - Thiết bị đo Qua kết đo thí nghiệm, sử dụng phương pháp mặt cắt gradient trung gian với AB = 500m, MN = 60m để thi cơng cho tồn vùng - Kỹ thuật đo đạc thực địa: Kỹ thuật đo đạc thực địa tuân thủ theo quy định phương pháp lý lịch máy Các giá trị đo đạc tính tốn tự động máy đổ trực tiếp vào máy tính để xử lý - Kiểm tra đánh giá chất lượng tài liệu: Việc đánh giá chất lượng tài liệu phương pháp đo thực theo quy phạm công tác địa vật lý điện Chúng đo kiểm tra độc lập đoạn tuyến T13 với tỷ lệ đo kiểm tra 5% Sai số điểm đo tính theo công thức sau: δ i = 2( ρ ρ k1i k1i −ρ ) +ρ k2i × 100% k2i n Cơng thức tính sai số cho tồn vùng là: δ = ∑δ I n Trong : - ρk1i, ρk2i giá trị đo lần thứ lần thứ hai điểm đo thứ i - n tổng số điểm đo kiểm tra Cơng thức tính sai số cho giá trị pha tương tự Kết qủa tính tốn sai số cho hệ số pha 3.7% điện trở 2,4% Các sai số nhỏ sai số cho phép, chất lượng tài liệu đủ độ tin cậy để xử lý luận giải địa chất 32 32 2.2.3 Phương pháp đo sâu phân cực kích thích Phương pháp đo sâu phân cực kích thích tiến hành đoạn tuyến T0P (tuyến T0P qua điểm nước nóng với phương vị 46 o) đoạn tuyến T6 Máy móc sử dụng chế độ đo đạc thu thập số liệu tiến hành tương tự phương pháp mặt cắt phân cực kích thích Chúng tơi đo sâu phân cực hệ thiết bị cực đối xứng Wenner, AB max = 900m, MNmax = 300m khoảng cách điểm đo tuyến 20m Công tác xử lý tổng hợp tài liệu 3.1 Tách trường địa phương khu vực Các trường địa vật lý đo phản ánh yếu tố địa chất trường tổng cộng Để giải nhiệm vụ cụ thể từ trường tổng cần phải tách thành phần riêng biệt liên quan đến đối tượng nghiên cứu Chúng tơi sử dụng phương pháp trung bình trượt cải tiến để tách trường khu vực địa phương" Thơng qua việc xác định bán kính tương quan dị thường chúng tơi chọn kích thước cửa sổ để tách trường khu vực địa phương 150m để tách trường khu vực địa phương cho giá trị trường từ (Ta) Kết tách trường xây dựng thành sơ đồ đồ thị tỷ lệ 1: 10.000 1: 5.000 3.2 Phân tích định lượng tài liệu điện theo mơ hình hai chiều (2D) 3.2.1 Giải toán thuận 2D Lý thuyết trường điện dòng chiều phần lý thuyết điện động lực, thoả mãn hệ phương trình Maxwell tính chất định luật ơm Kyrchop Phương trình tốn phân bố điện U với nguồn điện chiều I có dạng: Div ( σgradU ) = − Iδ q (3.2-1) Với điều kiện biên tính liên tục hàm ranh giới phẳng độ dẫn điện σ suy giảm hàm theo khoảng cách từ nguồn đến điểm quan sát: ∂U ∂z = z=0 ∂U ∂n + U cos θ / r = (3.2-2) Trong đó: σ độ dẫn điện mơi trường, độ nguồn I), θ góc pháp tuyến δq n xung Dirac (đặc trưng cho đơn vị biên r , r khoảng cách từ nguồn đến biên 33 33 Để giải toán (3.2-1) (3.2-2) phân bố chiều (2D), từ luận điểm vật lý khác người ta dẫn đến hai hướng giải phương pháp phương trình tích phân phương pháp phương trình vi phân Tuy nhiên, phương pháp sử dụng rộng rãi hiệu phương pháp sai phân hữu hạn Phương pháp sai phân hữu hạn dựa diễn giải trực tiếp phép tính vi phân phương trình (3.2-1) cho phần tử mơi trường có diện tích S ij độ dẫn σ ij với đặc trưng điện nguồn điểm dịng I phân bố mơi trường gây (hình 3.2-1) I=1 N-1 N J= U ij ∆z1 ∆z 4 ρ ij = σ ij ρ ij ∆z M −1 ∆z M −1 ∆x1 ∆x ∆x N −1∆x N Hình 3.2-1 Lưới rời rạc theo lược đồ sai ph©n Viết phương trình (3.2-1) theo sơ đồ sai phân cho phần tử S ij phần tử lân cận đảm bảo điều kiện biên (3.2-2) ta hệ phương M-1 ~ trình đại số tuyến tính với hàm phổ điểm nút: ~ ~ U ij ~ ẩn số đặc trưng cho trường ~ ~ C e ( i, j ) U i +1, j + C s ( i, j ) U i, j +1 + C w ( i, j ) U i −1, j + C u ( i, j ) U i , j −1 − C p ( i, j ) U i , j = − Iε i , j M ε (3.2-3) ; i, j , x s ∈ Pij , , x s ∉ Pij Phương trình (3.2-3) viết dạng ma trận: i = 1, L , j − 1, M =1 ~ CU = b C- , b- (3.2-4) ma trận hệ số phụ thuộc vào kích thước tính chất dẫn điện véc tơ số hạng Iε ij 34 34 S ij Giá trị hàm U ( x, z ) tính thơng qua hàm phổ việc giải phương trình (3.2-4), phép biến đổi fourier ngược: ( ) ( ~ sin K y ∂ y U x, K y , z ∞  U ( x, y , z ) = ∫ πy ∂K y Ky ~ U , có từ ) dK y  (3.2-5) - hàm Macdonal bậc bậc theo công thức chuỗi (A.Ango, 1965) Từ giá trị hàm U ij xác định nút lưới rời rạc nguồn I gây ra, ta tính giá trị điện trở suất phân cực 3.2.2 Giải tốn ngược Mục đích phân tích định lượng mặt cắt sâu điện hai chiều (2D) xác định tham số phần tử toàn mặt cắt cho trường quan sát mặt cắt trùng với trường tính tốn cách tốt nhất, thực chất cực tiểu hóa phiếm hàm độ lệch bình phương trung bình giá trị điện trở suất, phân cực fi đo toàn tuyến đo  n [ f i − g ( x + δx ) ]  ∇F ( x + δ x ) = ∇  ∑  g i2 ( x)  n i =1  g i (x) với giá trị tính mơ hình lý thuyết : 0,5 →0 (3.2-6) Trong : x - Tham số mơ hình, tọa độ điện trở suất hay giá trị phân cực phần tử điện tích Khi giá trị S ij , δx - δx đủ bước thay đổi tham số x nhỏ, hàm g ( x + δx ) viết lại dạng đa thức: g ( x + δx) = g ( x ) + Jδx Trong J ma trận Jacobi (3.2-4), xác định cho lần thay đổi mơ hình Độ lệch số liệu thực tế mơ hình ban đầu xác định biểu thức : ε = f − g ( x) Phiếm hàm (3.2-6) đạt cực tiểu * * ε ≈ Jδx , từ ta xác định bước thay đổi δx = J ε , với J ma trận nghịch đảo J Các tham số mơ hình điều chỉnh sau lần tính lặp dừng kết hai lần tính liên tiếp nhỏ giá trị sai số cho trước Việc cực tiểu hoá phiếm hàm sai số (3.2-6) có nhiều thuật tốn giải khác nhau, với tốc độ hội tụ nhanh đạt đến độ sai số 1% số liệu tính lý thuyết sau 5- 10 lần tính lặp 35 35 Với thuật tốn nêu trên, cho phép lập chương trình tự động hố q trình tính, giới có số phần mềm phân tích 2D tương đối phổ biến hiệu phần mềm RES2DINV, RESIX IP2DI Chúng tơi sử dụng chương trình RESIX IP2DI để phân tích tài liệu đo sâu phân cực Kết công tác địa vật lý 4.1 Phương pháp từ Kết đo từ theo điện tích 14 tuyến xử lý để tách trường khu vực với bán kính 150m trường dư thể sơ đồ đồ thị giá trị trường từ (Hình vẽ số 1) Kết đo từ diện tích nghiên cứu cho thấy trường từ thay đổi phạm vi nhỏ, chủ yếu biến đổi khoảng biên độ từ 10 đến 30nT chúng phản ánh diện phân bố đá cuội kết, cát kết, bột kết, cát kết quắc zít hệ tầng Suối Bàng hệ tầng Tây Chang, không khoanh định dải dị thường Theo kết đo từ diện tích nghiên cứu khơng có hoạt động magma nguồn nước nóng khơng liên quan với magma 4.2 Kết phương pháp mặt cắt phân cực Các kết đo mặt cắt phân cực hệ thiết bị gradient trung gian với kích thước AB=500m, MN=60m, tần số 0,5Hz đoạn tuyến (T2, T0B, T0 T0A) thể sơ đồ đồ thị giá trị điện trở pha tỷ lệ 1: 5.000 1: 2.500 (Hình vẽ số 2) Kết đo mặt cắt cho thấy giá trị điện trở suất biểu kiến (ρk) thay đổi phạm vi lớn từ đến 2100Ωm, giá trị pha thay đổi từ 15 đến 80 mrad Kết đo mặt cắt phân cực khoanh định hai dải dị thường có giá trị điện trở suất biểu kiến nhỏ 30Ωm giá trị pha lớn 40 mrad Dải dị thường số có quy mơ lớn, chiều rộng 200m, chiều dài khống chế 1,5km phản ánh hệ thống đứt gãy sâu theo phương tây bắc - đông nam đứt gãy liên quan trực tiếp với nguồn nước nóng Pom Lót Dải dị thường thứ có bề rộng gần 100m, phản ánh đứt gãy theo phương gần bắc - nam khơng liên quan với nguồn nước nóng Pom Lót 4.3 Phương pháp đo sâu phân cực 4.3.1 Kết đo sâu phân cực tuyến Top Kết đo sâu phân cực xoay chiều hệ thiết bị cực đối xứng Wenner với khoảng điểm đo tuyến 20m, đoạn tuyến T0p từ cọc -64 đến đến cọc 36 36 -32 (điểm lộ nước nóng khe suối cọc -50 (-500m hình vẽ)) phân tích định lượng chương trình 2D (RESIX IP2DI) thể hình 3, hình vẽ số Kết đo sâu phân cực cho phép ta theo dõi khoanh định đới dị thường điện trở suất thấp nhỏ 30Ωm có giá trị pha lớn 50mrad, với kích thước rộng gần 200m tồn đến độ sâu 100m (vì vướng sơng, khơng có đường rải dây nên khống chế đến độ sâu 100m) Dị thường Dị thường phản ánh đới đứt gãy tây bắc đơng nam, liên quan với tầng nước nóng Pom Lót 4.3.2 Kết đo sâu phân cực tuyến T6 Kết đo sâu phân cực xoay chiều đoạn tuyến T6 từ cọc -42 đến cọc -12 phân tích định lượng chương trình 2D (RES2DINV) thể hình 5, hình vẽ số Kết đo sâu phân cực cho phép ta theo dõi khoanh định đới dị thường điện trở suất thấp nhỏ 25Ωm có giá trị pha lớn 80mrad, với kích thước gần 200m, tồn đến độ sâu 150m (vì vướng sơng, khơng có đường rải dây nên khống chế đến độ sâu 150m) Dị thường phản ánh đới đứt gãy tây bắc đông nam, liên quan với nguồn nước nóng Pom Lót Đới đứt gãy cắm dốc nghiêng đơng bắc 4.4 Kết tính lượng sinh nhiệt phóng xạ Một số nguyên nhân sinh nguồn địa nhiệt phân hủy phóng xạ đất đá lớp vỏ đất Để xác định nguyên nhân gây nguồn địa nhiệt vùng Pom Lót xem chúng có phải q trình phân rã phóng xạ lớp vỏ đất khoảng chiều dày 10 km? chúng tơi tính lượng sinh nhiệt phóng xạ gây dựa vào kết đo mật độ, hàm lượng U, Th K mẫu đá diện tích nghiên cứu Các ngun tố phóng xạ tự nhiên đất đá có thời gian bán phân rã từ 109 đên 1010 năm nguồn cung cấp nhiệt vô tận, phân hủy khối lượng hạt nhân nguyên tử nguyên tố phóng xạ giải phóng lượng lớn mà chủ yếu nhiệt Tất đồng vị tự nhiên trình phân rã sinh nhiệt theo phương trình: Qsn-px = σ ( α Cu + β CTh + γ Ck ) Kw/ km3 Trong : σ mật độ đá chứa U,Th, K tính kg/ km3 37 37 Cu hàm lượng U mẫu đá tính theo ppm CTh hàm lượng Th mẫu đá tính theo ppm Ck hàm lượng K mẫu đá tính theo % Khi nhân phóng xạ U, Th, K phân hủy giải phóng lượng định để đạt tới cân phóng xạ ứng với nguyên tố có số sinh nhiệt định gọi số sinh nhiệt phóng xạ Qua thực nghiệm người ta xác định số sau: α - số sinh nhiệt U 9,525 x 10- kw/kg β - số sinh nhiệt Th 2,561 x 10- kw/kg γ - số sinh nhiệt K 3,477 x 10- kw/kg Thay giá trị trên,Với đơn vị đo mật độ (σ) g/cm3 ta có cơng thức tính sau: Qsn-px = 0,01 σ (9,525.Cu + 2,561.CTh + 3,477.Ck ) Kw/ km3 Kết tính lượng sinh nhiệt mẫu diện tích nghiên cứu đưa bảng sau: Kết tính lượng sinh nhiệt phóng xạ Bảng1 Hàm lượng phóng xạ Lượng Số Mật độ sinh U TT hiệu Th (ppm) nhiệt (g/cm3) K (%) (ppm) mẫu phóng xạ 500 2,18 2,10 4,80 7,80 1.6 556 2,50 1,70 4,50 7,80 1.7 558 2,46 1,70 4,40 13,50 2.0 553 2,41 2,10 4,90 17,80 2.4 554 2,43 2,30 5,20 22,00 2.8 559 2,41 1,80 5,50 9,30 2.0 550 2,24 1,70 4,60 7,80 1.6 534 2,42 1,20 3,40 9,30 1.5 535 2,66 1,60 5,50 10,70 2.3 10 542 2,44 1,90 5,00 9,30 1.9 11 522 2,44 2,00 3,70 9,30 1.6 12 523 2,38 1,80 3,80 10,70 1.7 13 532 2,29 1,70 5,40 12,10 2.0 14 576 1,93 1,80 5,20 9,30 1.5 15 580 1,98 1,70 3,50 10,70 1.3 16 586 2,07 2,00 3,40 13,50 1.5 38 38 Tên đá cát kết xericit cát kết cát bột kết cát kết đa khoáng cát kết đa khoáng cát kết bột kết bột kết cát kết bột kết Cát bột kết cát kết bột kết bột kết bột kết phiến sét TT 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 Số Hàm lượng phóng xạ Lượng Mật độ U hiệu Th (ppm) sinh (g/cm3) K (%) (ppm) mẫu nhiệt 528 2,25 2,10 4,00 10,70 1.6 545 2,34 1,70 5,30 14,90 2.2 547 2,30 1,80 3,40 13,50 1.7 527 2,22 1,70 6,00 10,70 2.0 531 2,22 1,70 4,10 13,50 1.8 621 2,33 1,70 4,40 14,90 2.0 619 1,98 1,70 2,80 13,50 1.3 602 2,19 1,90 5,20 9,30 1.8 603 2,28 1,70 5,50 7,80 1.8 606 2,35 1,70 5,30 7,80 1.8 607 2,15 2,10 5,20 9,30 1.7 609 2,42 1,70 4,30 12,10 1.9 610 2,38 1,80 5,50 10,70 2.0 634 2,10 1,80 4,70 7,80 1.5 633 2,32 1,70 4,10 13,50 1.8 631 2,13 1,70 5,30 13,50 1.9 624 2,06 1,80 3,40 12,10 1.4 612 2,26 1,50 4,90 10,70 1.8 680 2,40 1,70 5,80 9,30 2.0 683 2,01 1,80 5,00 13,50 1.8 684 2,12 1,50 5,00 13,50 1.9 658 2,20 1,30 4,10 13,50 1.7 659 2,29 1,30 4,30 10,70 1.7 695 2,16 2,00 4,30 12,10 1.7 690 2,04 1,90 5,50 10,70 1.8 691 2,26 1,70 5,70 9,30 1.9 638 2,58 1,60 5,40 12,10 2.3 639 2,57 1,70 5,70 13,50 2.4 700 2,41 1,70 5,40 12,10 2.1 702 2,45 2,00 5,60 19,20 2.7 32/2 2,71 1,80 5,50 32,00 3.8 23/2 2,44 2,10 4,90 22,00 2.7 21/2 2,34 1,90 4,60 20,60 2.4 27/2 2,68 2,10 6,10 23,50 3.4 36/2 2,53 1,30 3,50 32,00 3.0 71/2 2,56 1,40 3,60 6,40 1.4 37/2 2,80 1,40 3,70 7,80 1.7 25/2 2,64 1,30 3,60 5,00 1.4 34/2 2,43 1,30 3,70 6,40 1.4 29/2 2,20 1,70 5,00 16,40 2.1 63/2 2,71 1,40 3,70 6,40 1.5 39 39 Tên đá bột kết cát kết bột kết cát kết Cát bột kết bột kết bột kết bột kết bột kết cát kết cát sạn kết dăm kết cát kết bột kết bột kết bột kết bột kết cát kết cát kết bột kết cát kết bột kết cát kết bột kết bột kết cát kết cát kết bột kết cát kết bột kết dăm kết bột kết cát kết, bột kết dăm kết đá vôi đá hoa? đá phiến đá vôi TT 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 Số Hàm lượng phóng xạ Lượng Mật độ U hiệu Th (ppm) sinh (g/cm3) K (%) (ppm) mẫu nhiệt 400/2 2,03 1,90 4,10 17,80 1.9 30/2 2,74 1,30 3,60 7,80 1.6 46/2 2,67 2,00 4,40 17,80 2.5 703 2,14 1,80 4,60 9,30 1.6 704 2,49 1,80 5,20 10,70 2.1 705 2,23 1,80 3,50 12,10 1.6 706 2,17 1,60 3,70 13,50 1.6 707 2,28 1,70 5,00 12,10 1.9 708 2,16 1,50 4,70 9,30 1.6 709 2,36 1,80 4,30 12,10 1.8 710 2,18 2,20 5,50 9,30 1.8 711 2,51 1,60 5,30 10,70 2.1 675 2,55 2,20 5,30 12,10 2.3 692 2,30 2,00 5,80 10,70 2.1 511 2,49 1,60 8,10 9,30 2.7 512 2,34 1,80 8,70 9,30 2.6 513 2,50 1,90 4,10 13,50 2.0 Tên đá bột kết đá vôi Cát bột kết bột kết cát kết bột kết bột kết cát kết bột kết cát kết bột kết cát kết cát kết bột kết cát kết cát kết cát kết dạng 75 514 2,52 2,00 4,70 20,60 2.6 quăczit 76 515 2,14 2,00 6,10 17,80 2.4 Bột kết 77 516 2,18 2,00 6,10 16,40 2.3 Bột kết 78 517 2,35 2,00 6,00 14,90 2.4 Bột kết cát kết dạng 79 519 2,52 2,00 5,80 23,50 3.1 quăczit 80 536 2,34 1,60 5,50 10,70 2.0 cát kết Từ kết tính tốn bảng cho thấy loại đất đá diện tích nghiên cứu có lượng nhiệt sinh phóng xạ thay đổi từ 1,3 đến 3,8 kw/km3, nhiệt lượng nhỏ thơng lượng nhiệt trung bình vỏ đất (thơng lượng trung bình đất 10Kw/ km 3), khơng có khả làm sơi nước mà góp phần bổ sung cho lượng nhiệt mà Như nguyên nhân gây nguồn nước nóng diện tích nghiên cứu khơng phải phóng xạ từ vỏ mà nhiệt từ sâu lên theo đứt gãy sâu trẻ tái hoạt động chế tuần hoàn đối lưu 40 40 KẾT LUẬN Từ kết thu thập tài liệu, việc xây dựng tổ hợp phương pháp điều tra nghiên cứu đánh giá nguồn địa nhiệt phương pháp Viễn Thám, Địa chất, Địa Vật lý” Thăm dị tìm kiếm địa nhiệt cần phải áp dụng tổ hợp phương pháp: Viễn Thám, Địa Chất Địa Vật lí, kết phương pháp cho ta tham số ứng với đối tượng địa chất khác Việc lựa chọn tổ hợp phương pháp Viễn Thám, Địa Chất Địa Vật lí thăm dị, tìm kiếm địa nhiệt phải hợp lí với tình hình địa chất, đối tượng cần nghiên cứu nhằm nâng cao hiệu công tác tiết kiệm chi phí Ứng dụng nguồn địa nhiệt thiết thực, hiệu kinh tế an tồn với mơi trường, cụ thể Viện Khoa học Địa chất Khoáng sản sử dụng nguồn địa nhiệt đời sống sinh hoạt hàng ngày cần áp dụng, phổ biến sâu rộng toàn Quốc Trên kết thu tập tài liệu tổng hợp cá nhân, kính trình Viện nghiệm thu Trân trọng cảm ơn Hà nội, ngày 20 tháng 10 năm 2018 Tác giả Bùi Xuân Khải 41 41 TÀI LIỆU DẪN [1*] Võ cơng nghiệp, 2011, cần có nhìn mực tiềm địa nhiệt Việt Nam, Tạp chí khoa hoc Trái đất [1] Baskov E.A., Surikov S N., 1975: Thủy nhiệt chỏm cầu Thái Bình Dương Trái Đất, Nxb “Nedra”, Moskva (Tiếng Nga), 200tr [2] Baskov E.A., Surikov S N., 1989: Thủy nhiệt Trái Đất, Nxb “Nedra”, Moskva (Tiếng Nga), 176tr [3] Hoàng Văn Chước (chủ biên), 1989: Báo cáo kết đề tài 52C-05-01 “Nghiên cứu sử dụng lượng địa nhiệt Hội Vân để sấy” Lưu trữ Đại học Bách khoa Hà Nội [4] Gadalia A., 1982 : Projet géothermique de moyenne énergie dans la province de Nghĩa Bình BRGM Orléans, France [5] Cao Duy Giang (chủ biên), 1998: Báo cáo tổng kết đề tài nghiên cứu “Đánh giá tiềm địa nhiệt vùng Bắc Trung Bộ Việt Nam” Lưu trữ Viện Khoa học Địa chất Khoáng sản, Hà Nội [6] Cao Duy Giang (chủ biên), 2004: Báo cáo tổng kết đề tài nghiên cứu “Đánh giá tiềm địa nhiệt vùng Tây Bắc Bắc Bộ triển vọng sử dụng chúng phục vụ phát triển kinh tế - xã hội” Lưu trữ Viện Khoa học Địa chất Khoáng sản, Hà Nội [7] Lại Thế Huyến (chủ biên), 1986: Báo cáo kết nghiên cứu đề tài “Sử dụng lượng địa nhiệt nhiệt thấp vào việc sấy sản phẩm” Viện Năng lượng điện khí hóa, Hà Nội [8] Koenig J., 1981: Geothermal resources exploration and exploitation (in Vietnam) Project document VIE/80/025, Hanoi [9] Thomas Mathew, 2008: Study on the Socio - economic framework for the use of Geothermal Energy in Vietnam FUGRO CONSULT GMBH Germany Procedings of the 8th Asian Geothermal symposium Ha Noi, - 10 December, 1619 [10] Maunder B.R., 1993: Geoscientific reconnaissance of selected Geothermal Areas South Central coastal region of Vietnam Wellington, KRTA Ltd New Zealand [11] Muraoka H., 2008: Development of a small and low - temperature geothermal power generation system and its marketability in Asia Procedings of the 8th Asian Geothermal symposium Ha Noi, - 10 December, 13 -15 [12] Trương Minh, 1999 Về chế độ địa nhiệt bể trầm tích thềm lục địa Việt Nam Tạp chí Kinh tế Địa chất Ngun liệu khống, 18, (2), 24-29 [13] Võ Công Nghiệp (chủ biên), 1988: Báo cáo tổng kết đề tài nghiên cứu 44-0404 “Đánh giá tài nguyên địa nhiệt làm sở thiết kế khai thác sử dụng thử nghiệm vào mục đích lượng số vùng triển vọng” Lưu trữ Địa chất, Hà Nội [14] Võ Công Nghiệp, 1997: Một vài kết nghiên cứu thử nghiệm sử dụng lượng địa nhiệt nguồn địa nhiệt nóng Hội Vân Tuyển tập “Địa chất - địa nhiệt Địa chất công trình miền Trung” Số Nha Trang [15] Võ Cơng Nghiệp (chủ biên), 1998: Danh bạ nguồn địa nhiệt khống 42 42 địa nhiệt nóng Việt Nam Cục Địa chất Khoáng sản Việt Nam Hà Nội, 200tr [16] Võ Công Nghiệp, 2003: Lại bàn phương hướng sử dụng lượng địa nhiệt Việt Nam Tuyển tập “Địa chất Khoáng sản”, số Viện Khoa học Địa chất Khoáng sản Hà Nội,346 -349 [17] ORMAT, 1995: Pre-feasibility study report for the development of the geothermal resources of the SR of Vietnam Prepared by ORMAT International, INC - USA [18] Hoàng Hữu Quý (chủ biên), 1995: Báo cáo “nghiên cứu đánh giá tiềm địa nhiệt lãnh thổ từ Quảng Nam - Đà Nẵng đến Bà Rịa - Vũng Tàu” Lưu trữ Viện nghiên cứu Địa chất Khoáng sản, Hà Nội [19] Renewable Energy Essentials 2010: geothermal WWW.IEA.org 5p [20] Schochet D.N., 1997: Vietnam geothermal project overview, description of Energy conversion and project detail Scientific seminar on the Potential for Development of Geothermal Energy in the Central regions of the SR of Vietnam” Hanoi, October, [21] Đinh Văn Toàn, 1994: Kết bước đầu nghiên cứu phân bố nhiệt độ thạch vùng trũng Hà Nội Tạp chí Các khoa học Trái Đất T.16, 2, Hà Nội, 12-21 [22] Đinh Văn Toàn, Trịnh Việt Bắc, Nguyễn Trọng Yêm, 1996 Kết bước đầu xác định giá trị dòng nhiệt Nam Việt Nam Tạp chí “Các khoa học Trái Đất” T.18, 2, Hà Nội, 59-65 [23] Toshihiro Uchida, 2008: The International Geothermal Association and the Western Pacific Regional Branch; collaboration among geothermists in Asian countries The 8th Asian Geothermal Symposium Ha Noi, Viet Nam, December 10/2008, 45-49 [24] Tournaye D., 1990 : Proposition pour l’étude et la valorisation des ressources gộothermales de la RS du Vietnam Compagnie franỗaise pour le Développement de la Géothermie et des Energies Nouvelles Orléans, France [25] Cao Đình Triều, 2009: Mơ hình vận tốc sóng dọc P thạch manti Đơng Nam Á Tạp chí Địa chất, A/314, 36 -42 [26] Đồn Văn Tuyến, Đinh Văn Toàn, 2010 : Đánh giá tiềm giải pháp sử dụng địa nhiệt để phát triển lượng tái tạo Việt Nam Tuyển tập Hội nghị khoa học kỷ niệm 35 năm Viện KH&CN Việt Nam Hà Nội, 10/2010, 168 -172 [27] Nguyễn Trọng Yêm, Duchkov A D, 1991: Preliminary results of heat flow studies in North Vietnam, 2nd Conference on Geology of Indochina Hanoi, Vietnam Nov 11 - 13/1991 43 43 ... nguyên tố có số sinh nhiệt định gọi số sinh nhiệt phóng xạ Qua thực nghiệm người ta xác định số sau: α - số sinh nhiệt U 9,525 x 10- kw/kg β - số sinh nhiệt Th 2,561 x 10- kw/kg γ - số sinh nhiệt K... -172 [27] Nguyễn Trọng Yêm, Duchkov A D, 1991: Preliminary results of heat flow studies in North Vietnam, 2nd Conference on Geology of Indochina Hanoi, Vietnam Nov 11 - 13/1991 43 43 ... (in Vietnam) Project document VIE/80/025, Hanoi [9] Thomas Mathew, 2008: Study on the Socio - economic framework for the use of Geothermal Energy in Vietnam FUGRO CONSULT GMBH Germany Procedings