Thạch địa hóa là gì? Bản chất, phương pháp và ứng dụng

Thạch địa hóa (tiếng Anh: petrogeochemistry hoặc lithogeochemistry) là phân ngành khoa học liên ngành giữa thạch học và địa hóa học, chuyên nghiên cứu thành phần hóa học toàn đá, hàm lượng và quy luật phân bố của các nguyên tố chính, nguyên tố vết cũng như các hệ đồng vị trong các loại đá tự nhiên. Mục tiêu trọng tâm của thạch địa hóa là xác định nguồn gốc thạch luận, điều kiện hóa lý hình thành (như nhiệt độ, áp suất, độ fugacity oxy) và tái tạo bối cảnh kiến tạo của các thành tạo đá magma, đá biến chất và đá trầm tích trong lịch sử tiến hóa của vỏ Trái Đất.

Bản chất và đối tượng nghiên cứu của thạch địa hóa

Trong địa chất học hiện đại, thạch địa hóa đóng vai trò cầu nối định lượng giữa quan sát thạch học định tính dưới kính hiển vi và các mô hình động lực học manti - vỏ Trái Đất. Bằng cách phân tích hóa học mẫu đá tổng thể (whole-rock analysis), các nhà địa chất có thể nhận diện những dấu vết địa hóa còn lưu giữ lại từ quá trình nóng chảy từng phần của manti hoặc vỏ, quá trình kết tinh phân đoạn của magma, cũng như các tương tác đồng hóa magma và nhiễm bẩn vỏ.

Đối tượng khảo sát của thạch địa hóa bao gồm toàn bộ các kiểu đá trong thạch quyển:

• Đá magma: Nghiên cứu thành phần hóa học để xác định nguồn gốc magma nguyên sinh, mức độ tiến hóa của chuỗi magma, điều kiện tích tụ khoáng sản quặng liên quan và tái lập môi trường kiến tạo cổ (như cung núi lửa, sống núi giữa đại dương, hay đới nội mảng).
• Đá biến chất: Sử dụng địa hóa nguyên tố trơ và đồng vị để xác định đá gốc nguyên thủy (protolith), đồng thời ứng dụng các cặp khoáng vật cân bằng để tính toán áp suất và nhiệt độ biến chất (geothermobarometry).
• Đá trầm tích: Đánh giá mức độ phong hóa của đá mẹ, điều kiện oxy hóa khử của môi trường lắng đọng trầm tích, bối cảnh kiến tạo của lưu vực và nguồn cung cấp vật liệu trầm tích (provenance analysis).

Các nhóm nguyên tố trong phân tích thạch địa hóa

Dữ liệu thạch địa hóa học được phân tách thành ba nhóm thông tin cơ bản: oxit nguyên tố chính, nguyên tố vết và tỷ số đồng vị.

1. Nhóm oxit nguyên tố chính

Theo tài liệu chuyên khảo của Rollinson và Pease (2021), nhóm oxit nguyên tố chính (major elements) gồm 10 oxit cơ bản: $\mathrm{SiO_2}$ (silic dioxit), $\mathrm{TiO_2}$ (titan dioxit), $\mathrm{Al_2O_3}$ (nhôm oxit), $\mathrm{Fe_2O_3}$ và $\mathrm{FeO}$ (sắt oxit), $\mathrm{MnO}$ (mangan oxit), $\mathrm{MgO}$ (magie oxit), $\mathrm{CaO}$ (canxi oxit), $\mathrm{Na_2O}$ (natri oxit), $\mathrm{K_2O}$ (kali oxit) và $\mathrm{P_2O_5}$ (photpho pentoxit). Nhóm oxit này thường chiếm trên 99% khối lượng của hầu hết các loại đá magma không chứa nước.

Hàm lượng của các oxit chính phản ánh trực tiếp tổ hợp khoáng vật tạo đá chủ yếu (như thạch anh, felspat, pyroxen, olivin, amphibol, mica). Chỉ số $\mathrm{SiO_2}$ là thước đo phân loại cơ bản mức độ axit - bazơ của đá magma, chia thành đá siêu mafic ($\mathrm{SiO_2} < 45\%$), mafic ($45\% \le \mathrm{SiO_2} < 52\%$), trung tính ($52\% \le \mathrm{SiO_2} < 63\%$) và axit ($\mathrm{SiO_2} \ge 63\%$).

2. Nhóm nguyên tố vết

Nguyên tố vết (trace elements) là những nguyên tố có hàm lượng dưới 0,1% khối lượng (thường đo bằng ppm hoặc ppb). Dựa vào bán kính ion, điện tích và đặc tính địa hóa, nguyên tố vết được chia thành các phân nhóm chủ đạo:

• Nguyên tố thạch đới có bán kính ion lớn (LILE): Bao gồm $\mathrm{Rb}$ (rubidi), $\mathrm{Ba}$ (bari), $\mathrm{Sr}$ (stronti), $\mathrm{Cs}$ (cesi). Nhóm này có độ linh động cao trong dung dịch nhiệt dịch và quá trình phong hóa, thường phản ánh mức độ tham gia của chất lưu trong đới hút chìm.
• Nguyên tố có cường độ trường điện cao (HFSE): Bao gồm $\mathrm{Zr}$ (zirconi), $\mathrm{Hf}$ (hafni), $\mathrm{Nb}$ (niobi), $\mathrm{Ta}$ (tantal), $\mathrm{Th}$ (thori), $\mathrm{U}$ (urani), $\mathrm{Ti}$ (titan). Nhóm HFSE có tính trơ địa hóa, ít bị di động bởi dung dịch nước nhiệt dịch nên giữ vai trò tối quan trọng để nhận diện đá bị biến đổi.
• Nguyên tố đất hiếm (REE): Gồm chuỗi 15 nguyên tố nhóm lantan từ $\mathrm{La}$ (lantan) đến $\mathrm{Lu}$ (luteti). REE được chia thành đất hiếm nhẹ (LREE) và đất hiếm nặng (HREE). Đặc tính phân đoạn REE (ví dụ tỷ số $\mathrm{(La/Yb)_N}$) và dị thường Europi ($\mathrm{Eu/Eu^*}$) phản ánh trực tiếp sự kết tinh phân đoạn của plagioclas, granat hoặc amphibol trong buồng magma.
• Nguyên tố nhóm chuyển tiếp (TTE): Gồm $\mathrm{Ni}$ (niken), $\mathrm{Cr}$ (crom), $\mathrm{Co}$ (coban), $\mathrm{V}$ (vanadi), $\mathrm{Sc}$ (scandi). Các nguyên tố này tương thích cao với olivin, pyroxen và spinen trong các pha kết tinh sớm từ magma manti.

3. Hệ đồng vị phóng xạ và đồng vị bền

Các hệ đồng vị phóng xạ như $\mathrm{^{87}Sr/^{86}Sr}$, $\mathrm{^{143}Nd/^{144}Nd}$, $\mathrm{^{176}Hf/^{177}Hf}$ và các đồng vị chì ($\mathrm{Pb}$) cung cấp dấu vân tay đồng vị đặc trưng. Vì tỷ số đồng vị của các nguyên tố nặng không bị phân đoạn đáng kể bởi quá trình kết tinh phân đoạn, chúng là công cụ độc lập giúp xác định thành phần nguồn phát sinh magma (manti cạn kiệt, manti giàu hóa, hay vỏ lục địa cổ).

Phương pháp luận và các biểu đồ phân loại thạch địa hóa

Để trích xuất thông tin định lượng từ hàng chục tham số hóa học của mỗi mẫu đá, thạch địa hóa sử dụng hệ thống các biểu đồ nhị biến, tam giác và biểu đồ đa nguyên tố chuẩn hóa.

Biểu đồ chuỗi tiến hóa oxit chính (Biểu đồ Harker)

Biểu đồ Harker biểu diễn tương quan giữa hàm lượng oxit $\mathrm{SiO_2}$ trên trục hoành với các oxit khác ($\mathrm{TiO_2, Al_2O_3, FeO^*, MgO, CaO, Na_2O, K_2O}$) trên trục tung. Xu hướng biến thiên tuyến tính hoặc phi tuyến trên biểu đồ Harker cho phép nhận biết quy luật kết tinh phân đoạn của các khoáng vật cụ thể: ví dụ sự giảm đồng thời của $\mathrm{MgO}$ và $\mathrm{FeO}$ theo chiều tăng $\mathrm{SiO_2}$ phản ánh sự tách sớm của olivin và pyroxen.

Biểu đồ phân loại đá núi lửa TAS và sơ đồ nguyên tố trơ

Hệ thống phân loại đá núi lửa quốc tế dựa trên biểu đồ TAS (Total Alkali - Silica) biểu diễn tổng kiềm $\mathrm{(Na_2O + K_2O)}$ theo $\mathrm{SiO_2}$. Tuy nhiên, đối với các đá cổ bị biến chất hoặc biến đổi nhiệt dịch, các oxit kiềm thường bị rửa trôi hoặc làm giàu thứ sinh. Trong công trình kinh điển năm 1977, Winchester và Floyd đã đề xuất sơ đồ phân loại dựa trên các nguyên tố trơ không di động, sử dụng tỷ số $\mathrm{Zr/TiO_2}$ đại diện cho mức độ axit hóa và tỷ số $\mathrm{Nb/Y}$ đại diện cho tính kiềm của magma. Sơ đồ này cho phép định danh chính xác bản chất đá gốc ban đầu (basalt, andesit, dacit hay rhyolit) ngay cả khi đá đã trải qua biến chất mức độ cao.

Biểu đồ phân định bối cảnh kiến tạo (Pearce et al., 1984)

Nghiên cứu của Pearce, Harris và Tindle (1984) đã xây dựng các biểu đồ phân định môi trường kiến tạo cho đá granit dựa trên các cặp nguyên tố vết trơ và bán trơ như $\mathrm{Rb}$ tương quan với $\mathrm{(Y + Nb)}$, hoặc $\mathrm{Nb}$ tương quan với $\mathrm{Y}$. Sơ đồ phân chia các khối granit thành 4 bối cảnh kiến tạo rõ rệt:

1. Granit cung núi lửa (VAG - Volcanic Arc Granites): Hình thành tại các đới hút chìm mảng đại dương dưới mảng lục địa hoặc cung đảo, đặc trưng bởi hàm lượng $\mathrm{Nb}$ và $\mathrm{Y}$ thấp, giàu nguyên tố $\mathrm{Rb}$ và $\mathrm{Ba}$.
2. Granit nội mảng (WPG - Within-Plate Granites): Liên quan đến hiện tượng tách giãn nội mảng, điểm nóng manti (hotspot) hoặc rifting lục địa, đặc trưng bởi hàm lượng $\mathrm{Nb, Y, Ta, Ga}$ rất cao và tính chất loại A (A-type granite).
3. Granit va chạm (COLG - Collision Granites): Sinh ra trong các giai đoạn va mảng lục địa - lục địa tạo núi, đặc trưng bởi hàm lượng $\mathrm{Rb}$ rất cao do nóng chảy vật liệu vỏ trầm tích giàu sét.
4. Granit sống núi giữa đại dương (ORG - Ocean Ridge Granites): Tồn tại cục bộ trong các phức hợp ophiolit tại trung tâm tách giãn đáy đại dương, có hàm lượng các nguyên tố không tương thích rất thấp.

Biểu đồ nhện đa nguyên tố chuẩn hóa (Spidergram)

Để so sánh hàm lượng nguyên tố vết và REE giữa các mẫu đá khác nhau mà không bị ảnh hưởng bởi hiệu ứng độ phổ biến tự nhiên lẻ - chẵn của nguyên tố (hiệu ứng Oddo - Harkins), dữ liệu phân tích được chuẩn hóa theo thành phần manti nguyên thủy (primitive mantle) hoặc chondrit theo giá trị chuẩn của McDonough và Sun (1995). Đồ thị chuẩn hóa cho thấy các hình thế dị thường đặc trưng: ví dụ dị thường âm $\mathrm{Nb-Ta}$ và $\mathrm{Ti}$ là dấu hiệu đặc trưng không thể nhầm lẫn của magma cung hút chìm (subduction-related magmas).

Ứng dụng thạch địa hóa trong nghiên cứu địa chất Việt Nam

Tại Việt Nam, phương pháp thạch địa hóa đóng vai trò then chốt trong việc giải mã lịch sử tiến hóa địa chất phức tạp của khối Đông Dương, đới khâu Sông Mã, đới xiết trượt Sông Hồng và đới uốn nếp uốn nếp Trường Sơn.

Một ví dụ điển hình là công trình nghiên cứu thạch địa hóa của Trần Trọng Hòa và cs. (2015) trên các phức hệ đá magma axit thuộc bồn trũng Tú Lệ và vòm nâng Phan Si Pan ở khu vực Tây Bắc Việt Nam. Dữ liệu địa hóa nguyên tố chính, nguyên tố vết và đồng vị đã chứng minh các thành tạo magma peralkaline và granitoid loại A tại đây có tuổi Permi muộn (khoảng 255 đến 260 triệu năm), mang đặc trưng thạch địa hóa của bối cảnh nội mảng liên quan trực tiếp đến hoạt động của chùm manti thuộc Tỉnh magma lớn Emeishan (Emeishan Large Igneous Province). Phát hiện thạch địa hóa này đã góp phần định vị chính xác ranh giới phía tây nam của mảng Nam Trung Hoa và cơ chế dịch trượt dọc theo đới đứt gãy Sông Hồng trong kỷ Kainozoi.

So sánh các kỹ thuật phân tích thạch địa hóa

Dưới đây là bảng so sánh các kỹ thuật phân tích dụng cụ hiện đại dùng để thu thập số liệu thạch địa hóa trong phòng thí nghiệm:

Kỹ thuật phân tích	Đối tượng phân tích	Giới hạn phát hiện	Ưu điểm nổi bật	Hạn chế chính
Quang phổ huỳnh quang tia X (XRF)	Oxit nguyên tố chính và một số nguyên tố vết ($\mathrm{Ba, Sr, Rb, Zr}$)	Khoảng vài ppm đến mức phần trăm khối lượng	Độ chính xác và độ lặp lại rất cao đối với oxit chính; chuẩn bị mẫu bằng ép viên hoặc nung chảy đĩa thủy tinh	Giới hạn phát hiện kém với các nguyên tố đất hiếm nặng và nguyên tố vết hàm lượng siêu thấp
Khối phổ plasma cảm ứng cao tần (ICP-MS)	Toàn bộ nguyên tố vết, đất hiếm (REE) và nguyên tố nhóm platin (PGE)	Mức ppb hoặc ppt	Độ nhạy cực cao; phân tích đồng thời hầu hết các nguyên tố vết trong bảng tuần hoàn từ một dung dịch mẫu hòa tan	Quy trình hòa tan mẫu bằng axit mạnh ($\mathrm{HF, HNO_3}$) phức tạp; khó hòa tan triệt để các khoáng vật bền nhiệt như zircon
Vi phân tích đầu dò điện tử (EPMA)	Thành phần hóa học tại điểm vi mô trên từng khoáng vật riêng lẻ	Khoảng 0,01% khối lượng (100 ppm)	Độ phân giải không gian cao (vùng bắn phá tia điện tử khoảng 1 đến 2 micromet); không phá hủy mẫu lát mỏng	Chỉ đo được nguyên tố chính và nguyên tố phụ trong khoáng vật, không đo được nguyên tố vết hàm lượng thấp
Laser Ablation ICP-MS (LA-ICP-MS)	Nguyên tố vết và tỷ số đồng vị tại điểm vi mô trên từng hạt khoáng vật (như zircon, monazit, apatit)	Mức ppm đến sub-ppm trên diện tích điểm bắn 10 đến 50 micromet	Cho phép định tuổi U-Pb đồng thời với đo đạc hàm lượng nguyên tố vết và đồng vị Hf trên cùng một hạt khoáng vật	Yêu cầu mẫu chuẩn đối chiếu phù hợp với nền khoáng vật đo; hiệu ứng phân đoạn chùm laser

Hạn chế, điều kiện áp dụng và ranh giới khoa học

Mặc dù là công cụ định lượng mạnh mẽ, việc phiên giải dữ liệu thạch địa hóa đòi hỏi sự cẩn trọng nghiêm ngặt và phải luôn gắn liền với nghiên cứu thạch học thực địa:

• Hiện tượng biến đổi nhiệt dịch và phong hóa: Mẫu đá thu thập ngoài thực địa hiếm khi hoàn toàn tươi nguyên. Sự biến đổi thứ sinh làm di động các nguyên tố kiềm ($\mathrm{Na, K, Rb, Ba, Sr}$), dẫn đến việc xác định sai vị trí trên các biểu đồ TAS hoặc Harker. Giá trị tổn thất khi nung (Loss on Ignition - LOI) cao hơn 2% đến 3% là dấu hiệu cảnh báo cần loại bỏ các nguyên tố linh động và chỉ dựa vào nhóm nguyên tố trơ (HFSE, REE).
• Hiệu ứng tổng hằng số (Closure Problem / Constant Sum Effect): Dữ liệu thành phần oxit chính được biểu diễn dưới dạng tỷ lệ phần trăm khối lượng với tổng chuẩn hóa bằng 100%. Hiệu ứng toán học này khiến cho sự biến thiên của một oxit chiếm hàm lượng lớn (như $\mathrm{SiO_2}$) sẽ tự động ép hàm lượng phần trăm của các oxit còn lại giảm xuống một cách biểu kiến, tạo ra các tương quan giả tạo không phản ánh bản chất hóa lý thực tế.
• Tính đa nguyên sinh của các biểu đồ phân định kiến tạo: Các ranh giới trường kiến tạo trên biểu đồ thạch địa hóa (như biểu đồ Pearce 1984) được xây dựng dựa trên các mẫu đá magma hiện đại có bối cảnh địa động lực đã biết rõ. Khi áp dụng cho các thành tạo đá cổ Tiền Cambri hoặc Paleozoi, các quá trình phức tạp như đồng hóa vỏ kèm kết tinh phân đoạn (AFC), sự tái nóng chảy của vỏ lục địa dị chất, hoặc sự biến dạng sau magma có thể đưa thành phần đá rơi vào trường kiến tạo sai lệch.

Do đó, nguyên tắc tối thượng trong thạch địa hóa học là không bao giờ kết luận nguồn gốc hay bối cảnh kiến tạo chỉ dựa trên một biểu đồ đơn lẻ, mà phải luôn đối chiếu tổng thể giữa quan sát thạch học lát mỏng, địa chất cấu tạo thực địa, địa hóa đồng vị và địa thời học đồng vị phóng xạ.