Tổng quan luận án
Luận án "Bài toán ngược trong Lý thuyết thế vị" được tác giả Chu Đức Khánh thực hiện dưới sự hướng dẫn khoa học của Giáo sư Tiến sĩ Đặng Đình Áng và Phó Tiến sĩ Nguyễn Bích Huy, bảo vệ năm 1996 tại Trường Đại học Sư phạm Thành phố Hồ Chí Minh, thuộc chuyên ngành Toán Giải tích (Mã số: 1.01).
Lý thuyết thế vị phát sinh từ thế kỷ XIX như một bộ phận trọng yếu của Vật lý Toán, nghiên cứu trường thế của các lực tự nhiên như trường hấp dẫn và điện từ trường. Trong thực tế, nhiều bài toán kỹ thuật thuộc Vật lý Địa cầu, Trắc địa, Điện học và Cơ học môi trường liên tục được mô hình hóa dưới dạng các bài toán ngược. Tuy nhiên, phần lớn các bài toán này không chỉnh theo nghĩa của J. Hadamard: chúng có thể không tồn tại nghiệm, có vô số nghiệm, hoặc nghiệm không phụ thuộc liên tục vào dữ kiện đầu vào (một nhiễu nhỏ trong dữ kiện đo đạc có thể dẫn đến sai số rất lớn ở nghiệm, thậm chí làm bài toán trở nên vô nghiệm).
Khoảng trống nghiên cứu then chốt được luận án xác định là: trong các bài toán ngược không chỉnh, việc đánh giá sai số của các phương pháp chỉnh hóa thường gặp nhiều khó khăn giải tích và chỉ được giải quyết riêng lẻ cho từng dạng bài toán cụ thể. Đồng thời, việc khảo sát tính duy nhất nghiệm đóng vai trò tiên quyết nhằm kiểm tra xem dữ kiện đo đạc có đủ để xác định nghiệm duy nhất hay không, hoặc xác định nghiệm chính xác riêng biệt (nghiệm có chuẩn nhỏ nhất khi bài toán có vô số nghiệm) để làm cơ sở xây dựng nghiệm xấp xỉ ổn định.
Mục tiêu nghiên cứu của luận án bao gồm:
- Thiết lập mô hình toán học giải tích chặt chẽ dưới dạng các phương trình và hệ phương trình tích phân cho các bài toán ngược trong lý thuyết thế vị.
- Khảo sát điều kiện duy nhất nghiệm cho các bài toán ngược tuyến tính và phi tuyến, chứng minh sự tồn tại và duy nhất của nghiệm có chuẩn nhỏ nhất trong không gian Hilbert đối với các bài toán có vô số nghiệm.
- Xây dựng phương pháp chỉnh hóa và thiết lập đánh giá sai số cho các bài toán ngược tuyến tính thông qua bài toán moment tổng quát trên không gian Hilbert vô hạn chiều cho cả trường hợp dữ kiện chính xác và dữ kiện có nhiễu.
- Xây dựng phương pháp chỉnh hóa cho các bài toán ngược phi tuyến xác định miền nguồn thông qua lý thuyết bậc tôpô và phương pháp xấp xỉ không gian Banach bằng các không gian con hữu hạn chiều.
Đối tượng nghiên cứu của luận án là các bài toán ngược đối với phương trình Laplace và hàm thế vị trong không gian Euclid hai chiều ($\mathbb{R}^2$) và ba chiều ($\mathbb{R}^3$). Phạm vi khảo sát bao gồm: bài toán thác triển số liệu đo trường điều hòa, bài toán xác định phân bố mật độ nguồn, bài toán xác định lỗ hổng (vết nứt) bên trong vật thể, và bài toán xác định hình dạng miền nguồn sinh trường thế.
Tổng quan tài liệu và vị trí của luận án
Lý thuyết trường thế vị được xây dựng từ các định luật vật lý cổ điển như Định luật vạn vật hấp dẫn của Newton và Định luật tĩnh điện của Coulomb, dẫn đến phương trình đạo hàm riêng Laplace và lớp hàm điều hòa. Khái niệm về tính chỉnh (well-posedness) và không chỉnh (ill-posedness) của các bài toán vật lý toán được định hình từ các công trình kinh điển của J. Hadamard.
Đối với các bài toán không chỉnh, các công trình của A. N. Tikhonov và nhiều nhà toán học giải tích trên thế giới đã đặt nền móng cho lý thuyết chỉnh hóa (regularization), đưa ra các nguyên lý tìm nghiệm xấp xỉ ổn định thay thế cho nghiệm không ổn định ban đầu. Trong lĩnh vực giải tích phi tuyến và tôpô vi phân, định lý của J. Cronin về mối liên hệ giữa bậc tôpô (topological degree) và số nghiệm của phương trình phi tuyến mở ra hướng tiếp cận định lượng cho việc nghiên cứu sự tồn tại và ổn định nghiệm của các phương trình toán tử phi tuyến.
Vị trí của luận án nằm ở điểm giao giữa giải tích hàm hiện đại, lý thuyết phương trình đạo hàm riêng và phương pháp tính cho bài toán ngược. Luận án giải quyết khoảng trống học thuật bằng cách tiếp cận có hệ thống từ khâu mô hình hóa, khảo sát tính duy nhất nghiệm (hoặc nghiệm có chuẩn cực tiểu), cho đến xây dựng thuật toán chỉnh hóa với đánh giá sai số tường minh cho cả hai trường hợp dữ liệu: dữ liệu lý tưởng (chính xác) và dữ liệu thực tế (bị nhiễu).
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng các hệ thống lý thuyết và công cụ giải tích toán học:
- Lý thuyết hàm thế vị và giải tích điều hòa: Khảo sát thế vị Newton $U(x,y,z) = \int_D \frac{\rho(\xi,\eta,\zeta)}{r} d\xi d\eta d\zeta$ trong $\mathbb{R}^3$ và thế vị logarit $U(x,y) = \int_D \rho(\xi,\eta) \ln\frac{1}{r} d\xi d\eta$ trong $\mathbb{R}^2$. Sử dụng các tính chất vi phân của hàm thế, đẳng thức Poisson $\Delta U = -4\pi\rho$ (hoặc $-2\pi\rho$), tính trơn trong các không gian Sobolev $W^{2,p}(G)$ và $H^2(G)$, nguyên lý cực đại của hàm điều hòa, công thức tích phân Green và định lý duy nhất nghiệm của bài toán Cauchy cho phương trình Laplace.
- Giải tích hàm và hình học không gian Hilbert: Không gian hàm khả tích bình phương $L^2$, cấu trúc không gian afin lồi đóng $p + K$, định lý hình chiếu trên tập lồi đóng để xác định nghiệm có chuẩn nhỏ nhất. Biến đổi Fourier trong $L^2(\mathbb{R})$ và $L^2(\mathbb{R}^2)$ kết hợp đại số tích chập.
- Lý thuyết bài toán moment: Thiết lập và giải bài toán moment tổng quát trên không gian Hilbert vô hạn chiều làm công cụ chỉnh hóa các toán tử tích phân tuyến tính.
- Lý thuyết bậc tôpô và giải tích phi tuyến: Định lý J. Cronin về bậc tôpô để kiểm soát số nghiệm và tính ổn định của phương trình toán tử phi tuyến trên không gian Banach, kết hợp kỹ thuật xấp xỉ Galerkin/chiếu trên các không gian con hữu hạn chiều.
Nguồn tư liệu và dữ liệu nghiên cứu là các mô hình toán học giải tích thuần túy, xuất phát từ các bài toán đo đạc dị thường trọng lực, từ trường địa cầu, và phương pháp thăm dò điện từ trong kỹ thuật kiểm tra không phá hủy.
Nội dung chính theo từng chương
Chương 1: Tổng quan về bài toán ngược trong Lý thuyết thế vị
Chương này trình bày cơ sở giải tích của lý thuyết thế vị và thiết lập mô hình toán học cho các bài toán ngược. Tác giả xuất phát từ định luật Newton biểu diễn lực tương tác giữa hai chất điểm có khối lượng $m_1, m_2$ cách nhau khoảng $r$:
$$F = k \frac{m_1 m_2}{r^2}$$
và định luật Coulomb cho lực tĩnh điện giữa hai điện tích. Trường lực sinh bởi nguồn $D \subset \mathbb{R}^3$ có phân bố khối lượng $\rho$ được biểu diễn bởi:
$$\vec{F}(x,y,z) = \int_D \rho(\xi,\eta,\zeta) \frac{\vec{r}}{r^3} d\xi d\eta d\zeta$$
với hàm thế vị Newton tương ứng là:
$$U(x,y,z) = \int_D \frac{\rho(\xi,\eta,\zeta)}{\sqrt{(x-\xi)^2 + (y-\eta)^2 + (z-\zeta)^2}} d\xi d\eta d\zeta$$
Đối với mặt phẳng $\mathbb{R}^2$, trường lực và hàm thế vị logarit được xây dựng tương ứng:
$$U(x,y) = \int_D \rho(\xi,\eta) \ln \frac{1}{\sqrt{(x-\xi)^2 + (y-\eta)^2}} d\xi d\eta$$
Chương 1 thiết lập 5 nhóm bài toán ngược cụ thể (gồm 7 bài toán từ $P_1$ đến $P_7$):
- Bài toán 1 & 2 (Thác triển số liệu đo trường điều hòa): Xác định hàm điều hòa $u$ trong nửa mặt phẳng trên $\mathbb{R}^2_+$ (Bài toán 1) hoặc nửa không gian trên $\mathbb{R}^3_+$ (Bài toán 2) khi chỉ biết giá trị của $u$ và đạo hàm pháp tuyến của nó trên một tập mở bị chặn $I \times {k}$ hoặc $B \times {k}$ ($k > 0$). Bài toán được chuyển về hệ phương trình tích phân Poisson theo ẩn hàm $v$ trên biên đáy.
- Bài toán 3 & 4 (Tìm phân bố nguồn): Xác định mật độ nguồn $\rho$ trên miền bị chặn $\Omega \subset \mathbb{R}^3$ khi biết thế vị $V_\Omega$ trên một tập mở bị chặn $D \subset \mathbb{R}^3$ (Bài toán 3) hoặc biết $V_\Omega$ và $\frac{\partial V_\Omega}{\partial z}$ trên một mặt mở $\omega \subset \mathbb{R}^2 \times {0}$ (Bài toán 4). Cả hai bài toán dẫn đến phương trình tích phân Fredholm loại 1.
- Bài toán 5 (Xác định lỗ hổng bên trong vật thể): Xác định cặp $(\Omega, u)$ trong đó $\Omega$ có biên ngoài $\Gamma$ đã biết và biên trong $\gamma$ (biên các lỗ hổng/vết nứt) chưa biết, thỏa mãn hệ phương trình:
$$\begin{cases} \Delta u = 0 & \text{trong } \Omega \ u|{\Gamma_0} = f, \quad \frac{\partial u}{\partial n}\Big|{\Gamma_0} = g \ u|_\gamma = 0 \end{cases}$$
- Bài toán 6 & 7 (Xác định miền nguồn): Xác định biên trên của miền nguồn $\Omega$ có phân bố mật độ thuần nhất $\rho = 1$ trong $\mathbb{R}^2$ (Bài toán 6) và $\mathbb{R}^3$ (Bài toán 7) khi biết biên dưới và biết đạo hàm pháp tuyến của thế vị $\frac{\partial V_\Omega}{\partial y}$ hoặc $\frac{\partial V_\Omega}{\partial z}$ trên mặt phẳng chuẩn ($y=0$ hoặc $z=0$). Đây là các bài toán biên tự do phi tuyến.
Chương 2: Sự duy nhất nghiệm
Chương 2 khảo sát tính duy nhất nghiệm cho các bài toán đã thiết lập ở Chương 1, phân tách thành hai nhóm tuyến tính và phi tuyến:
| Mã bài toán |
Dạng toán tử |
Bản chất hình học / Dữ kiện |
Kết quả về tính duy nhất nghiệm |
| Bài toán 1 |
Tuyến tính |
Thác triển điều hòa Cauchy trong $\mathbb{R}^2_+$ |
Nghiệm duy nhất trong $L^2(\mathbb{R})$ |
| Bài toán 2 |
Tuyến tính |
Thác triển điều hòa Cauchy trong $\mathbb{R}^3_+$ |
Nghiệm duy nhất trong $L^2(\mathbb{R}^2)$ |
| Bài toán 3 |
Tuyến tính |
Tìm phân bố nguồn $\rho$ từ dữ kiện trên $D \subset \mathbb{R}^3$ |
Vô số nghiệm; duy nhất nghiệm chuẩn cực tiểu trong $L^2(\Omega)$ |
| Bài toán 4 |
Tuyến tính |
Tìm phân bố nguồn $\rho$ từ dữ kiện trên mặt $z=0$ |
Vô số nghiệm; duy nhất nghiệm chuẩn cực tiểu trong $L^2(\Omega)$ |
| Bài toán 5 |
Phi tuyến |
Xác định biên trong $\gamma$ của lỗ hổng/vết nứt |
Nhiều nhất một nghiệm $(\Omega, u)$ khi $f \ne 0$ hoặc $g \ne 0$ |
| Bài toán 6 |
Phi tuyến |
Xác định biên trên miền nguồn 2 chiều ($n=2$) |
Nghiệm duy nhất $\phi \in C(\bar{D})$ |
| Bài toán 7 |
Phi tuyến |
Xác định biên trên miền nguồn 3 chiều ($n=3$) |
Nghiệm duy nhất $\phi \in C(\bar{D})$ |
Chương 3: Chỉnh hóa một số bài toán ngược tuyến tính
Chương này giải quyết trọn vẹn việc xây dựng thuật toán chỉnh hóa cho các bài toán ngược tuyến tính thông qua công cụ giải tích hàm:
- Bài toán moment tổng quát trên không gian Hilbert (Tiết 3.1): Xét không gian Hilbert vô hạn chiều $H$. Bài toán moment yêu cầu tìm phần tử $u \in H$ thỏa mãn hệ phương trình vô hạn $\langle u, \phi_n \rangle = \mu_n$ ($n = 1, 2, \dots$). Tác giả thiết lập thuật toán chỉnh hóa, chứng minh tính hội tụ và đưa ra các đánh giá sai số chính xác cho hai trường hợp:
- Trường hợp dữ kiện chính xác ${\mu_n}$.
- Trường hợp dữ kiện bị nhiễu ${\mu_n^\epsilon}$ với mức độ sai số $|\mu - \mu^\epsilon| \le \epsilon$.
- Chỉnh hóa bài toán thác triển số liệu đo (Tiết 3.2): Hệ phương trình tích phân từ Bài toán 1 và Bài toán 2 được chuyển đổi về dạng bài toán moment trong không gian Hilbert. Bằng cách áp dụng kết quả của Tiết 3.1, tác giả thu được nghiệm chỉnh hóa xấp xỉ ổn định của trường điều hòa trên toàn miền biên, kèm theo chặn sai số giải tích theo tham số chỉnh hóa và mức nhiễu dữ kiện.
- Chỉnh hóa bài toán tìm phân bố nguồn (Tiết 3.3): Phương trình tích phân Fredholm loại 1 của Bài toán 3 và Bài toán 4 được chuyển về bài toán moment trên không gian Hilbert $L^2(\Omega)$. Thuật toán chỉnh hóa hướng trực tiếp đến việc xấp xỉ nghiệm duy nhất có chuẩn nhỏ nhất $\rho_0$ đã được xác lập ở Chương 2, đảm bảo tính ổn định số trị khi dữ liệu đo đạc bị nhiễu.
Chương 4: Chỉnh hóa một số bài toán ngược phi tuyến
Chương 4 tập trung xây dựng nghiệm chỉnh hóa cho bài toán phi tuyến xác định miền nguồn (Bài toán 6 trong $\mathbb{R}^2$ và Bài toán 7 trong $\mathbb{R}^3$):
- Sử dụng Định lý của J. Cronin về mối liên hệ giữa bậc tôpô của ánh xạ và số nghiệm của phương trình toán tử phi tuyến trên không gian Banach.
- Áp dụng phương pháp xấp xỉ không gian Banach vô hạn chiều bằng dãy các không gian con hữu hạn chiều (phương pháp xấp xỉ hữu hạn chiều kiểu Galerkin/chiếu tôpô).
- Xây dựng nghiệm xấp xỉ ổn định trong cả hai tình huống:
- Dữ kiện biên của trường thế vị chính xác tuyệt đối.
- Dữ kiện biên bị nhiễu sai số thực nghiệm.
- Thiết lập các định lý về sự hội tụ của nghiệm xấp xỉ về nghiệm chính xác của bài toán biên tự do khi mức độ nhiễu tiến về 0 và số chiều xấp xỉ tăng dần.
Kết quả và những đóng góp mới
Luận án đạt được các kết quả giải tích mới:
- Thiết lập mô hình hóa toán học hoàn chỉnh: Chuyển đổi thành công các bài toán vật lý thế vị thực tế (thác triển trường điều hòa, tìm phân bố nguồn, phát hiện khuyết tật/vết nứt, xác định biên miền nguồn) thành các phương trình và hệ phương trình tích phân tuyến tính hoặc phi tuyến tương ứng.
- Giải quyết tường minh bài toán duy nhất nghiệm:
- Chứng minh tính duy nhất nghiệm cho bài toán thác triển Cauchy bằng công cụ biến đổi Fourier và tích chập.
- Phát hiện và chứng minh cấu trúc không gian nghiệm afin của bài toán phân bố nguồn, khẳng định sự tồn tại và duy nhất của nghiệm có chuẩn $L^2$ cực tiểu.
- Chứng minh tính duy nhất nghiệm của bài toán xác định lỗ hổng/vết nứt bằng phương pháp phản chứng dựa trên thác triển giải tích và nguyên lý cực đại.
- Chứng minh tính duy nhất của biên miền nguồn 2D và 3D thông qua kỹ thuật tích phân thế vị và công thức vi tích phân Newton-Leibnitz.
- Phát triển phương pháp chỉnh hóa cho lớp bài toán tuyến tính: Xây dựng khung giải thuật cho bài toán moment tổng quát trên không gian Hilbert vô hạn chiều và áp dụng thành công để tìm nghiệm xấp xỉ ổn định kèm đánh giá sai số cho bài toán thác triển điều hòa và bài toán tìm nguồn.
- Phát triển phương pháp chỉnh hóa cho lớp bài toán phi tuyến: Đề xuất quy trình giải xấp xỉ ổn định cho bài toán xác định miền nguồn dựa trên định lý bậc tôpô J. Cronin kết hợp xấp xỉ không gian Banach hữu hạn chiều, xử lý được cả trường hợp dữ liệu chính xác và dữ liệu có sai số.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Luận án tập trung vào các miền hình học chuẩn tắc (nửa mặt phẳng, nửa không gian, miền bị chặn có biên đủ trơn) và giả định mật độ nguồn thuần nhất ($\rho = 1$) trong bài toán xác định hình học miền nguồn. Các kết quả giải tích chủ yếu dừng ở mức độ chứng minh toán học thuần túy và đánh giá sai số lý thuyết, chưa triển khai thử nghiệm số trên máy tính với dữ liệu đo đạc địa chấn/địa vật lý thực tế quy mô lớn.
- Hướng nghiên cứu tiếp: Mở rộng bài toán cho các môi trường vật chất không đồng nhất hoặc dị hướng; khảo sát bài toán xác định đồng thời cả hình dạng miền nguồn lẫn mật độ nguồn biến thiên $\rho(x,y,z)$; phát triển các thuật toán số song song giải hệ phương trình moment và phương trình toán tử phi tuyến nhiều chiều.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên khảo cho:
- Nghiên cứu sinh, học viên cao học và giảng viên ngành Toán Giải tích, Toán Ứng dụng, Phương trình Đạo hàm riêng và Vật lý Toán: cung cấp kỹ thuật giải tích xử lý bài toán không chỉnh, phương pháp biến đổi Fourier trong bài toán Cauchy, kỹ thuật bậc tôpô phi tuyến J. Cronin và lý thuyết moment không gian Hilbert.
- Các nhà nghiên cứu trong lĩnh vực Vật lý Địa cầu, Trắc địa và Thăm dò khoáng sản: tham khảo mô hình toán học giải đoán nguồn gốc trường trọng lực, trường từ, và thác triển số liệu đo đạc bị giới hạn về không gian địa hình.
- Kỹ sư và chuyên gia Cơ học, Kiểm tra không phá hủy (NDT): tham khảo cơ sở lý thuyết của bài toán xác định vị trí, hình dạng lỗ hổng và vết nứt ngầm trong lòng cấu trúc vật liệu bằng phương pháp trường điện từ.
Câu hỏi thường gặp
1. Tại sao các bài toán ngược trong lý thuyết thế vị thường không chỉnh theo nghĩa Hadamard?
Trong các bài toán này, toán tử truyền dẫn từ nguồn/biên tới số liệu đo đạc thường có tính chất làm trơn (toán tử compact tích phân loại 1). Khi giải ngược lại, toán tử nghịch đảo không bị chặn, dẫn đến việc nghiệm không phụ thuộc liên tục vào dữ kiện: một sai số đo đạc rất nhỏ (nhiễu) sẽ bị khuếch đại vô hạn, làm cho nghiệm xấp xỉ dao động mất kiểm soát hoặc làm bài toán mất nghiệm.
2. Tác giả giải quyết trường hợp vô số nghiệm của bài toán tìm phân bố nguồn (Bài toán 3 & 4) như thế nào?
Tác giả chứng minh tập hợp tất cả các nghiệm của bài toán lập thành một không gian afin $p + K$, trong đó $K$ là không gian con đóng của không gian Hilbert $L^2(\Omega)$. Dựa vào tính chất hình học của không gian Hilbert, tập lồi đóng $p + K$ chứa duy nhất một phần tử $\rho_0$ có chuẩn nhỏ nhất. Luận án chọn phần tử duy nhất $\rho_0$ này làm đối tượng mục tiêu để xây dựng nghiệm chỉnh hóa.
3. Công cụ giải tích nào được sử dụng để chỉnh hóa các bài toán ngược tuyến tính trong Chương 3?
Tác giả quy đổi các phương trình tích phân tuyến tính của bài toán thác triển điều hòa và bài toán tìm nguồn về bài toán moment tổng quát trên không gian Hilbert vô hạn chiều. Sau đó, tác giả thiết lập nghiệm chỉnh hóa dựa trên cấu trúc của toán tử moment và đưa ra các đánh giá sai số định lượng cho cả trường hợp dữ liệu chính xác và dữ liệu có sai số $\epsilon$.
4. Phương pháp nào được luận án áp dụng để chỉnh hóa bài toán xác định miền nguồn phi tuyến trong Chương 4?
Tác giả kết hợp Định lý của J. Cronin về mối liên hệ giữa bậc tôpô với số nghiệm của phương trình phi tuyến và phương pháp xấp xỉ không gian Banach bằng dãy các không gian con hữu hạn chiều để xây dựng nghiệm xấp xỉ ổn định, chứng minh sự hội tụ của dãy nghiệm xấp xỉ về nghiệm thực của bài toán khi dữ kiện bị nhiễu.
Kết luận
Luận án "Bài toán ngược trong Lý thuyết thế vị" của tác giả Chu Đức Khánh là công trình nghiên cứu hệ thống về lớp bài toán ngược không chỉnh phát sinh từ vật lý toán và giải tích điều hòa. Tác giả đã hoàn thành việc mô hình hóa toán học, khảo sát tính duy nhất nghiệm hoặc nghiệm có chuẩn cực tiểu cho cả 7 bài toán ngược tuyến tính và phi tuyến cốt lõi. Bằng việc vận dụng lý thuyết bài toán moment trên không gian Hilbert và lý thuyết bậc tôpô J. Cronin, luận án đã xây dựng thành công các phương pháp chỉnh hóa xấp xỉ ổn định kèm đánh giá sai số chặt chẽ, đóng góp cơ sở lý thuyết giá trị cho chuyên ngành Toán Giải tích và các ứng dụng địa vật lý liên quan.