Luận án Tiến sĩ: Phân tích ứng xử đất nền móng GBS dưới tải trọng tổ hợp

Tổng quan về luận án

Công trình biển trọng lực (Gravity Based Structure - GBS hay CTBTL) giữ vai trò then chốt trong phát triển hạ tầng ngoài khơi thế giới và chiến lược kinh tế biển - an ninh quốc phòng của Việt Nam. Khác biệt với công trình móng cọc thép truyền thống, công trình biển trọng lực dựa hoàn toàn vào trọng lượng bản thân và khối dằn (ballast) để duy trì độ ổn định trượt, lật và sức chịu tải của nền đất đáy biển dưới các tác động khắc nghiệt từ sóng, gió, dòng chảy và động đất. Tại vùng biển thềm lục địa Việt Nam, đặc biệt là Bể trầm tích Nam Côn Sơn ở độ sâu nước trên $110\text{ m}$, việc tính toán thiết kế gặp thách thức lớn do điều kiện địa chất trầm tích phức tạp đan xen và chế độ thủy động lực học ngẫu nhiên.

Khoảng trống nghiên cứu (Research Gap) chính nằm ở sự thiếu hụt các mô hình giải tích chính xác để tính toán tải trọng sóng cực đại tác động lên trụ công trình có tiết diện thay đổi đặt trong vùng nước sâu, cũng như việc xác định trường ứng suất - chuyển vị tổng quát của nền đất đáy biển dưới tổ hợp tải trọng lệch tâm (lực thẳng đứng $P$, lực ngang $H$ và mômen lật $M$). Các phương pháp truyền thống thường đơn giản hóa móng chịu tải đúng tâm hoặc sử dụng các công thức kinh nghiệm rời rạc, chưa thiết lập được bản đồ trường chuyển vị đồng mức (contour domain) cho bán không gian vô hạn khi chịu tải phức hợp. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Địa chất (Mã số: 62.01) của tác giả Đặng Xuân Trường (2019) tại Trường Đại học Bách khoa – ĐHQG TP. Hồ Chí Minh, dưới sự hướng dẫn của NGƯT. PGS. Đậu Văn Ngọ và PGS. Nguyễn Hữu Bảng, đã giải quyết trực tiếp khoảng trống này.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

• Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để mô hình hóa và lượng hóa chính xác các thông số sóng ngẫu nhiên cực hạn tác động lên trụ công trình biển trọng lực có tiết diện thay đổi tại thềm lục địa Nam Việt Nam?
• Câu hỏi nghiên cứu 2 (RQ2): Lời giải giải tích tường minh nào có thể xác định chính xác ứng suất và chuyển vị của bán không gian vô hạn dưới đáy móng tròn/phức hợp khi chịu tác động đồng thời của tải trọng phân bố đều $q$ và mômen $M$?
• Câu hỏi nghiên cứu 3 (RQ3): Mức độ tương thích và độ tin cậy của mô hình giải tích mới so với phương pháp phần tử hữu hạn phi tuyến và cơ học đất truyền thống tại khu vực Mỏ Đại Hùng là bao nhiêu?
• Giả thuyết khoa học 1 (H1): Ứng dụng lý thuyết quy hoạch thực nghiệm Box - Wilson và hàm hồi quy đa biến cho phép xác định chính xác các thông số sóng ngẫu nhiên phục vụ tính toán tải trọng cực đại.
• Giả thuyết khoa học 2 (H2): Bằng cách phân tích hàm dưới dấu tích phân $\int x \tan(ax) dx$ thành chuỗi chứa các hệ số Bernoulli $B_{2n}$, hoàn toàn có thể thiết lập lời giải giải tích giải bài toán chuyển vị đàn hồi Flamant mở rộng cho đáy móng CTBTL.
• Giả thuyết khoa học 3 (H3): Trường chuyển vị đồng mức thiết lập từ giải pháp giải tích phản ánh chính xác trạng thái biến dạng lệch tâm và nguy cơ mất ổn định trượt cục bộ của đất nền đáy biển so với mô hình số PLAXIS.

Khung lý thuyết của nghiên cứu tích hợp: Lý thuyết sóng Airy, công thức thủy động lực Morison, lý thuyết đàn hồi bán không gian vô hạn (Flamant/Boussinesq), lý thuyết quy hoạch thực nghiệm Box - Wilson và mô hình vật liệu đàn dẻo Mohr - Coulomb. Phạm vi nghiên cứu tập trung vào vùng Mỏ Đại Hùng (Lô 05-1a, Bể Nam Côn Sơn) với độ sâu cột nước $110\text{ m}$, khai thác bộ số liệu địa tầng khoan khảo sát 10 lớp đất, các đường cong nén lún $e-p$ và dữ liệu sóng quan trắc bão Phase 2. Toàn văn công trình dày 119 trang, 15 trang phụ lục, 14 bảng biểu, 51 hình vẽ và 03 phần mềm tự động hóa chuyên dụng lập trình trên ngôn ngữ Visual Basic.

Literature Review và Positioning

Lịch sử phát triển công trình biển trọng lực ghi nhận bước ngoặt năm 1973 khi hãng DORIS (Pháp) lắp đặt thành công bồn chứa dầu trọng lực bê tông đầu tiên tại mỏ Ekofisk (Biển Bắc) ở độ sâu $71\text{ m}$ (DORIS, 1973). Tiếp nối thành công này, hàng loạt dàn khoan trọng lực bê tông khổng lồ dạng CONDEEP được phát triển bởi Olav Olsen và Norwegian Contractors (NC) tại Biển Bắc như Beryl A ($118\text{ m}$, 1975), Brent B ($140\text{ m}$, 1975), Statfjord A/B/C ($145\text{ m}$, 1977-1984), Gullfaks C ($216\text{ m}$, 1989) và đỉnh cao là Troll A ($303\text{ m}$, 1995) với chiều cao kết cấu $369{,}4\text{ m}$, khối lượng bê tông $221.000\text{ m}^3$ và diện tích đáy móng $16.600\text{ m}^2$ (Norske Shell, 1995). Các thiết kế dạng SEA TANK, ANDOC (Dunlin A, $153\text{ m}$, 1977), CG Doris và Ove Arup (Malampaya, Philippines, $43\text{ m}$, 2000; Hebron, Canada, $109\text{ m}$, 2015) đã khẳng định tính khả thi vượt trội của giải pháp móng trọng lực ở các vùng biển khắc nghiệt.

Trong y văn thế giới tồn tại hai luồng quan điểm và tranh luận khoa học lớn:

1. Luồng quan điểm giải tích giải tích cổ điển và bán thực nghiệm (Morison et al., 1950; Hogben & Standing, 1974; Meyerhof, 1953; Vesic, 1975): Ưu tiên sử dụng các biểu thức giải tích đơn giản hóa dựa trên lý thuyết đàn hồi nửa không gian và cân bằng giới hạn. Nhóm tác giả này giả định phân bố áp lực đáy móng tuyến tính hoặc phẳng, bỏ qua tương tác phi tuyến cục bộ và phân bố chuyển vị không gian đa hướng.
2. Luồng quan điểm mô phỏng số phi tuyến liên tục (Zienkiewicz et al., 1978; Prevost, 1989; Potts et al., 2001): Cho rằng chỉ có phương pháp phần tử hữu hạn (FEM) hoặc sai phân hữu hạn phi tuyến với các mô hình đất phức tạp (như Hardening Soil, Cam-Clay biến tính) mới mô tả đúng ứng xử chu kỳ (cyclic loading) và sự suy giảm cường độ cắt của nền đất biển sâu. Tuy nhiên, phương pháp số đòi hỏi chi phí tính toán cực lớn, phụ thuộc vào quá nhiều thông số đầu vào khó xác định chính xác và thiếu tính tổng quát hóa trong giai đoạn thiết kế cơ sở.

Tại Việt Nam, các công trình biển trọng lực chủ yếu dừng lại ở quy mô nhỏ ven bờ hoặc móng đèn biển tại Trường Sa (Nam Yết, Sơn Ca, Sinh Tồn từ năm 2007) và móng tuabin phong điện tại Bạc Liêu, Ninh Thuận. Thực tiễn các Nhà giàn DK1 xây dựng trong thập niên 1980 trên nền san hô chịu rung lắc nghiêm trọng do kết cấu móng cọc thép phá vỡ cấu trúc tinh thể san hô, buộc phải gia cố bằng các khối bê tông trọng lực vào cuối thập niên 1990 để chuyển đổi nguyên lý làm việc sang móng trọng lực (Nguyễn Hữu Bảng, Đậu Văn Ngọ, 2010).

So với các nghiên cứu quốc tế điển hình như Draugen Condeep (dàn 1 trụ lớn nhất thế giới ở độ sâu $251{,}3\text{ m}$) hay Hibernia CGS (Canada, $80\text{ m}$, chịu va chạm băng trôi), nghiên cứu của Đặng Xuân Trường định vị chính xác tại khoảng trống: Phát triển một công cụ giải tích toán học chuẩn xác kết hợp chuỗi Bernoulli để giải bài toán tích phân suy rộng của lý thuyết đàn hồi Flamant, đồng thời xây dựng thuật toán tích phân Simpson rời rạc hóa tải sóng Morison cho trụ tiết diện biến đổi trong vùng nước sâu $110\text{ m}$ tại Biển Đông. Đây là bước tiến quan trọng, cầu nối giữa lý thuyết giải tích giải tích và kiểm chứng phần tử hữu hạn số.

Đóng góp lý thuyết và khung phân tích

flowchart TD
    A["Điều kiện môi trường biển Mỏ Đại Hùng (Bể Nam Côn Sơn, H = 110m)"] --> B["Xác định tải trọng môi trường ngẫu nhiên"]
    A --> C["Khảo sát địa kỹ thuật 10 lớp đất & thí nghiệm e-p"]
    
    B --> B1["Quy hoạch thực nghiệm Box - Wilson & Hàm hồi quy sóng"]
    B --> B2["Lý thuyết sóng Airy & Công thức Morison mở rộng"]
    B2 --> B3["Tích phân số Simpson cho trụ tiết diện thay đổi"]
    
    B1 & B3 --> D["Tổ hợp tải trọng cực hạn (P, H, M)"]
    
    D & C --> E["Mô hình giải tích đàn hồi bán không gian vô hạn (Flamant)"]
    E --> E1["Khai triển chuỗi tích phân Elliptic qua hệ số Bernoulli B_2n"]
    E1 --> E2["Phần mềm chuyên dụng Visual Basic: Bản đồ miền đồng mức chuyển vị"]
    
    D & C --> F["Mô hình phần tử hữu hạn số PLAXIS (Mohr-Coulomb)"]
    
    E2 <-->|Đối chứng & Đánh giá sai số| F
    F --> G["Kết luận ứng xử nền móng & Tối ưu hóa kích thước đế móng CTBTL"]

Đóng góp cho lý thuyết

Luận án mang lại bước đột phá toán học trong cơ học môi trường liên tục và lý thuyết đàn hồi ứng dụng cho kỹ thuật địa kỹ thuật biển. Đóng góp lý thuyết độc bản bao gồm việc giải quyết tích phân suy rộng dạng Elliptic:

$$\int x \tan(ax) , dx$$

Trong cơ học vật rắn biến dạng, khi xác định chuyển vị của môi trường bán không gian vô hạn chịu tải trọng phân bố $q$ kết hợp mômen uốn $M$, tích phân trên xuất hiện trong biểu thức tích phân mặt giải tích Flamant mở rộng. Luận án đã chứng minh một cách giải tích tường minh rằng: “Tích phân $\int x \tan(ax) dx$ được lấy khi và chỉ khi phân tích hàm dưới dấu tích phân ra chuỗi, sau đó lấy tích phân thông qua các hàm chứa hệ số $B_{2n}$ của Bernoulli”.

Khai triển chuỗi được thiết lập:

$$\tan(ax) = \sum_{n=1}^{\infty} \frac{(-1)^{n-1} 2^{2n} (2^{2n}-1) B_{2n}}{(2n)!} (ax)^{2n-1}$$

Từ đó thiết lập biểu thức tích phân:

$$\int x \tan(ax) , dx = \sum_{n=1}^{\infty} \frac{(-1)^{n-1} 2^{2n} (2^{2n}-1) B_{2n}}{(2n)! (2n+1)} a^{2n-1} x^{2n+1} + C$$

Đây là lời giải giải tích đầu tiên được đề xuất và áp dụng hoàn chỉnh trong bài toán xác định chuyển vị nền đất móng công trình biển trọng lực, mở rộng lý thuyết đàn hồi cổ điển của Flamant và Boussinesq cho các trường hợp biên tải trọng phức tạp phi đối xứng.

Khung phân tích độc đáo

Khung phân tích của luận án là sự hợp nhất đa ngành giữa:

1. Lý thuyết thủy động lực học sóng phi tuyến: Áp dụng phương trình Morison kết hợp lý thuyết sóng Airy và công thức cầu phương Simpson để tích phân lực cắt $f_x(z)$ và mômen uốn lật $M(z)$ dọc theo chiều sâu trụ có đường kính biến đổi $D(z)$.
2. Lý thuyết tối ưu hóa và quy hoạch thực nghiệm Box - Wilson: Biến đổi các yếu tố khí tượng hải văn ngẫu nhiên (chiều cao sóng $H_s$, chu kỳ sóng $T$, vận tốc gió $V_w$) thành mô hình toán học dự báo cực trị đáng tin cậy.
3. Cơ học vật rắn biến dạng bán không gian vô hạn: Mô hình hóa sự tương tác giữa kết cấu móng tròn/chữ nhật độ cứng vô hạn với môi trường đất nền nhiều lớp quy đổi, thiết lập trường đẳng trị chuyển vị lún $S(r, \theta)$ và ứng suất tiếp xúc $\sigma_z(r, \theta)$.
4. Điều kiện biên giới hạn: Mô hình xác định rõ các điều kiện biên: tiếp xúc đáy móng không chịu kéo (no-tension interface), giới hạn góc ma sát trong $\varphi$, lực dính $c$, hệ số biến dạng đàn hồi $E_0$ và hệ số Poisson $\nu$ của từng phân lớp địa chất.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phản biện (Critical Realism). Thiết kế nghiên cứu là sự phối hợp chặt chẽ giữa suy diễn giải tích toán học (analytical deduction), mô phỏng số (numerical simulation) và kiểm chứng thực nghiệm dữ liệu địa chất thực địa mỏ Đại Hùng.

classDiagram
    class GravityBasedStructure {
        +Double Diameter_Base
        +Double Height_Total
        +Double Dead_Weight
        +Double Ballast_Weight
        +calculateHydrodynamicForces()
        +evaluateOverturningStability()
    }
    class WaveHydrodynamics {
        +Double Water_Depth_110m
        +Double Wave_Height_Hmax
        +Double Wave_Period_T
        +airyLinearTheory()
        +morisonEquationIntegration()
        +simpsonRuleSolver()
    }
    class SoilFoundationInteraction {
        +Double ElasticModulus_E
        +Double PoissonRatio_nu
        +Double Cohesion_c
        +Double FrictionAngle_phi
        +flamantSemiInfiniteSpace()
        +bernoulliSeriesExpansion()
        +generateSettlementContours()
    }
    class PlaxisVerification {
        +Mesh2D_Axisymmetric
        +MohrCoulombModel
        +StageConstructionLoad
        +extractSettlementProfile()
    }
    GravityBasedStructure --> WaveHydrodynamics : chịu tác động tải cực hạn
    GravityBasedStructure --> SoilFoundationInteraction : truyền ứng suất đáy móng
    SoilFoundationInteraction <--> PlaxisVerification : so sánh đối chuẩn kết quả

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu triển khai qua các giai đoạn nghiêm ngặt:

• Xử lý số liệu địa chất - địa tầng: Phân tích cột địa tầng tổng hợp Bể Nam Côn Sơn qua các giếng khoan thăm dò tại Mỏ Đại Hùng. Khu vực nghiên cứu trải qua 3 giai đoạn tiến hóa kiến tạo: Trước tạo rift (Paleocen - Eocen), Đồng tạo rift (Oligocen - Miocen sớm gồm Hệ tầng Cầu, Dừa) và Sau tạo rift (Miocen giữa - Đệ Tứ gồm Hệ tầng Thông - Mãng Cầu, Nam Côn Sơn và Biển Đông).
• Thu thập dữ liệu cơ lý đất nền: Xác định thông số của 10 lớp địa chất đặc trưng tại hố khoan Mỏ Đại Hùng (tọa độ địa lý trong Lô 05-1a). Khai thác các biểu đồ nén lún chuẩn $e-p$ từ thí nghiệm nén không nở hông oedometer để xác định hệ số nén lún $a$, chỉ số nén $C_c$, và môđun biến dạng $E_0$.
• Giải thuật tích phân số Simpson: Do trụ công trình biển có đường kính thay đổi theo chiều sâu (tiết diện lớn ở đáy đỡ thùng chứa dầu và thon dần lên mặt nước), phương trình Morison không thể lấy tích phân trực tiếp. Tác giả áp dụng công thức Simpson phân chia chiều sâu nước $H = 110\text{ m}$ thành $2n$ đoạn rời rạc để giải số lực sóng cực đại $F_{x\max}$ và mômen lật $M_{\max}$.
• Tự động hóa lập trình: Xây dựng trọn bộ mã nguồn Visual Basic để xuất ma trận chuyển vị và vẽ đường đồng mức chuyển vị đất nền.

Data và phân tích

Đặc tính địa tầng và cơ lý đất nền tại vị trí khảo sát Mỏ Đại Hùng được tổng hợp chính xác trong Bảng 1:

Dữ liệu tổ hợp tải trọng thiết kế tác động lên đáy móng công trình biển trọng lực tại mỏ Đại Hùng bao gồm:

• Tải trọng thẳng đứng cố định do trọng lượng bản thân và khối dằn: $P_v = 1.250.000\text{ kN}$ ($\approx 125.000\text{ tấn}$).
• Tải trọng ngang do sóng gió cực hạn 100 năm: $H_x = 185.000\text{ kN}$.
• Mômen lật tại đáy móng do sóng cực đại: $M_y = 6.450.000\text{ kN}\cdot\text{m}$.
• Đường kính đáy móng trọng lực giả định: $D = 90\text{ m}$ (diện tích đáy $A \approx 6.361{,}7\text{ m}^2$).

Mô hình phần tử hữu hạn PLAXIS 2D/3D Foundation được thiết lập với lưới phần tử tam giác 15 nút bậc cao, điều kiện biên ngàm cứng ở đáy và biên con lăn ở hai bên, mô phỏng quá trình thi công hạ thủy, dằn nước và đặt tải trọng bão cực hạn.

Phát hiện đột phá và implications

Những phát hiện then chốt

1. Quy luật phân bố bất đối xứng của trường chuyển vị dưới móng tròn: Dưới tác động đồng thời của tải thẳng đứng $P$ và mômen lật $M$, biểu đồ miền đồng mức chuyển vị lún không còn dạng đối xứng xuyên tâm tròn xoay mà chuyển dịch mạnh về phía mép chịu nén gia tăng. Độ lún cực đại tại mép móng phía đón sóng đạt $S_{\max} = 18{,}4\text{ cm}$, trong khi mép đối diện chỉ đạt $S_{\min} = 4{,}2\text{ cm}$, tạo ra độ lún lệch góc xoay $\theta = 0{,}00158\text{ rad}$.
2. Sự hình thành đới tập trung ứng suất cắt cục bộ ở tầng bùn sét Lớp 4: Mô hình giải tích và PLAXIS đồng quy kết luận rằng lớp bùn sét yếu (Lớp 4, dày $8{,}4\text{ m}$ ở độ sâu $15{,}5\text{ m} - 23{,}9\text{ m}$) là vị trí xung yếu nhất xuất hiện bề mặt trượt dẻo tiềm năng khi mômen sóng đạt đỉnh cực hạn, với hệ số an toàn trượt giảm từ $F_s = 2{,}45$ (khi tải đúng tâm) xuống $F_s = 1{,}38$ (khi có tổ hợp sóng bão 100 năm).
3. Độ chính xác vượt trội của chuỗi Bernoulli trong tích phân giải tích: So sánh giữa kết quả tính toán độ lún theo phương pháp tổng phân tố cổ điển ($S_{\text{tong}} = 19{,}8\text{ cm}$), giải tích chuỗi Bernoulli ($S_{\text{bernoulli}} = 18{,}1\text{ cm}$) và mô hình số PLAXIS ($S_{\text{plaxis}} = 17{,}6\text{ cm}$) cho thấy sai số giữa mô hình giải tích mới của luận án và mô hình số chỉ là $2{,}84%$, khẳng định độ tin cậy tuyệt đối của giải thuật toán học đề xuất.
4. Hiệu ứng giảm tải sóng cục bộ của trụ tiết diện biến đổi: Việc áp dụng tích phân Simpson chứng minh rằng việc thu hẹp đường kính trụ ở vùng nước dao động sóng mạnh (từ mặt nước tĩnh xuống $-30\text{ m}$) giúp giảm tổng lực cắt đáy tới $31{,}5%$ và giảm mômen uốn lật đáy tới $38{,}2%$ so với thiết kế trụ đường kính không đổi (Monotower trụ đứng thông thường).

graph LR
    subgraph Tải trọng tác động
        P["Tải nén đúng tâm P"]
        H["Tải ngang H"]
        M["Mômen lật M (Sóng 100 năm)"]
    end
    
    subgraph Ứng xử nền đất mỏ Đại Hùng
        L1["Lớp 1-3: Cát pha chặt vừa (Biến dạng nén nhanh)"]
        L4["Lớp 4: Bùn sét yếu (Nguy cơ trượt dẻo cục bộ)"]
        L5_10["Lớp 5-10: Sét cứng & Cát kết (Tầng đỡ chịu lực chính)"]
    end
    
    subgraph Kết quả chuyển vị
        S_max["S_max = 18.4 cm (Mép đón sóng)"]
        S_min["S_min = 4.2 cm (Mép khuất sóng)"]
        Tilt["Độ lún lệch theta = 0.00158 rad"]
    end
    
    P & H & M --> L1 --> L4 --> L5_10
    L5_10 --> S_max & S_min --> Tilt

Implications đa chiều

• Về mặt lý thuyết: Cung cấp công cụ toán học giải tích mạnh mẽ, chuẩn hóa phương pháp giải chuỗi tích phân trong cơ học môi trường liên tục và lý thuyết đàn hồi cho nền móng công trình biển.
• Về mặt phương pháp luận: Thiết lập quy trình chuẩn tích hợp giữa quy hoạch thực nghiệm thống kê (xử lý tải trọng ngẫu nhiên) và phần mềm tự động hóa tính toán trường chuyển vị, rút ngắn thời gian tính toán thiết kế cơ sở từ hàng tuần xuống vài phút.
• Ứng dụng thực tiễn kỹ thuật dầu khí và năng lượng: Cho phép các kỹ sư Petrovietnam, Vietsovpetro tối ưu hóa kích thước khối đế móng, giảm khối lượng bê tông và cốt thép từ $10 - 15%$, tiết kiệm hàng chục triệu USD cho mỗi kết cấu dàn khoan ngoài khơi. Đồng thời ứng dụng trực tiếp cho thiết kế móng trọng lực trụ tuabin điện gió ngoài khơi (Offshore Wind Turbines) tại Ninh Thuận, Cà Mau, Bạc Liêu.
• Ý nghĩa đặc biệt đối với Quốc phòng - An ninh: Cung cấp cơ sở khoa học trực tiếp để tính toán, thiết kế và thi công các công trình Ụ chìm bê tông trọng lực dưới đáy biển phục vụ mục đích phòng thủ quốc phòng, kiểm soát vùng biển đặc quyền kinh tế và bảo vệ chủ quyền tại Biển Đông.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận các giới hạn học thuật:

1. Giả định môi trường đàn hồi tuyến tính: Lời giải giải tích Flamant mở rộng giả định đất nền là môi trường bán không gian vô hạn đồng nhất đàn hồi tuyến tính đẳng hướng, chưa phản ánh đầy đủ tính chất dị hướng và sự biến dạng dẻo phi tuyến tích lũy theo thời gian dưới hàng triệu chu kỳ tải sóng (cyclic degradation).
2. Mô hình tương tác đất - móng đơn giản hóa: Phân tích chưa xét đến hiện tượng hở khe tiếp xúc (gapping) cục bộ giữa đáy móng bê tông và đất nền khi chịu mômen uốn cực lớn, cũng như chưa tích hợp hiệu ứng chân khay (skirt foundation) cắm sâu vào đất đáy biển.
3. Giới hạn số liệu động đất cục bộ: Số liệu rung chấn động đất biển sâu tại rãnh trũng Nam Côn Sơn chưa được quan trắc thực đo dài hạn mà chủ yếu dựa trên các giả định tiêu chuẩn thiết kế chung.

Chương trình nghiên cứu 10 năm tiếp theo định hướng 4 trọng tâm:

• Hướng 1: Phát triển mô hình giải tích - số tương tác đất - kết cấu phi tuyến hoàn toàn (Fully Coupled Non-linear SSI) xét đến sự suy giảm độ bền cắt chu kỳ và áp lực nước lỗ rỗng thặng dư do bão kéo dài.
• Hướng 2: Mở rộng nghiên cứu móng trọng lực kết hợp chân khay hút (Suction Caisson Gravity Foundation) cho vùng nước sâu $>200\text{ m}$.
• Hướng 3: Nghiên cứu hiện tượng xói mòn cục bộ (scour) và hóa lỏng (liquefaction) xung quanh đế móng do tương tác dòng chảy đáy biển và sóng bão.
• Hướng 4: Tích hợp trí tuệ nhân tạo (AI/Machine Learning) và mạng nơ-ron vật lý (PINNs) để giải quyết các phương trình vi phân chuyển vị nền móng thời gian thực.

Tác động và ảnh hưởng

Công trình của Đặng Xuân Trường đã tạo ra tác động học thuật và công nghiệp sâu rộng:

• Tác động học thuật: Đóng góp 02 bài báo trên tạp chí quốc tế chuyên ngành Scopus/ISI, 01 bài báo tạp chí khoa học chuyên ngành trong nước và 5 báo cáo kỷ yếu hội nghị quốc tế uy tín. Khai mở hướng nghiên cứu ứng dụng chuỗi Bernoulli trong cơ học đàn hồi địa kỹ thuật biển.
• Chuyển dịch ngành công nghiệp dầu khí & năng lượng tái tạo: Cung cấp cẩm nang tính toán móng trọng lực cho các đơn vị tư vấn thiết kế hàng hải, mở đường cho việc nội địa hóa thiết kế và thi công dàn khoan bê tông trọng lực tại Việt Nam mà trước đây hoàn toàn phụ thuộc vào nhà thầu nước ngoài.
• Đóng góp chiến lược quốc phòng biển: Tạo tiền đề khoa học vững chắc cho dự án công trình Ụ chìm biển sâu, nâng cao năng lực công trình phòng thủ kiên cố của Hải quân Việt Nam trên thềm lục địa Biển Đông.

Đối tượng hưởng lợi

• Nghiên cứu sinh & Nhà khoa học Địa kỹ thuật biển: Tiếp cận phương pháp toán học mới để giải tích phân suy rộng trong lý thuyết đàn hồi và quy trình tích hợp mô hình số PLAXIS với giải pháp giải tích.
• Chuyên gia thiết kế công trình biển (Offshore Structural Engineers): Sử dụng các thuật toán Visual Basic và công thức tích phân Simpson rút gọn để tính nhanh lực sóng Morison và trường lún lệch của móng GBS.
• Tập đoàn Dầu khí & Năng lượng tái tạo: Tối ưu hóa chi phí đầu tư xây dựng khối đế trọng lực, kiểm soát rủi ro lún trượt công trình mỏ dầu khí Đại Hùng, Bạch Hổ, Cá Rồng Đỏ và các trang trại phong điện ngoài khơi.
• Bộ Quốc phòng & Quân chủng Hải quân: Sở hữu cơ sở kỹ thuật xây dựng công trình ngầm, nhà giàn thế hệ mới và ụ chìm bảo vệ chủ quyền biển đảo.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc bản nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc bản nhất là việc giải thành công tích phân dạng Elliptic $\int x \tan(ax) dx$ thông qua việc phân tích chuỗi hàm kết hợp hệ số Bernoulli $B_{2n}$. Nghiên cứu đã mở rộng trực tiếp Lý thuyết đàn hồi bán không gian vô hạn của Flamant và Boussinesq từ bài toán tải trọng tập trung/đúng tâm sang bài toán chịu đồng thời tải phân bố đều $q$ và mômen lệch tâm uốn lật $M$, thiết lập công thức tường minh tính chuyển vị tại mọi điểm dưới đáy móng tròn và chữ nhật.

2. Sự đổi mới trong phương pháp luận thể hiện như thế nào khi so sánh với các nghiên cứu trước đây?

So với các nghiên cứu của Box - Wilson (1951) trong quy hoạch thực nghiệm thuần túy và Morison (1950) trong lực sóng trên trụ đơn đều, luận án đã:

• Tích hợp quy hoạch thực nghiệm hàm hồi quy đa biến để xác định bộ thông số sóng ngẫu nhiên cực hạn tại vùng biển sâu Nam Côn Sơn.
• Ứng dụng tích phân số Simpson rời rạc hóa tải trọng Morison cho trụ có hình học biến đổi phức tạp theo phương đứng, vượt qua giới hạn của các công thức giải tích giải tích cổ điển vốn chỉ áp dụng cho trụ thẳng đứng tiết diện không đổi.

3. Phát hiện nào trong dữ liệu địa kỹ thuật Mỏ Đại Hùng gây bất ngờ nhất?

Phát hiện bất ngờ và quan trọng nhất là sự tồn tại của tầng bùn sét yếu Lớp 4 (dày $8{,}4\text{ m}$) ở độ sâu $15{,}5\text{ m}$ dưới đáy biển. Khi chịu mômen uốn do sóng bão cực hạn, ứng suất cắt truyền sâu tạo ra vùng biến dạng dẻo tập trung tại lớp này, khiến độ lún mép móng đón sóng tăng vọt lên $18{,}4\text{ cm}$ (gấp hơn 4 lần mép khuất sóng), biến cơ chế phá hủy tiềm năng từ trượt phẳng đáy móng thành trượt cung tròn sâu cắt qua Lớp 4.

4. Luận án có cung cấp đầy đủ quy trình lặp lại nghiên cứu (Replication Protocol) không?

Có. Luận án công bố chi tiết:

• Toàn bộ tọa độ, chiều dày và 7 chỉ tiêu cơ lý của 10 lớp đất tại hố khoan Mỏ Đại Hùng.
• Các bước toán học khai triển chuỗi Bernoulli và thuật toán tích phân Simpson.
• Lưu đồ thuật toán và mã nguồn chương trình Visual Basic tính toán thông số sóng, áp lực sóng cực đại và vẽ đường đồng mức chuyển vị.
• Các bước thiết lập thông số vật liệu và điều kiện biên trong phần mềm PLAXIS.

5. Khung chương trình nghiên cứu 10 năm tới được phác thảo ra sao?

Khung nghiên cứu 10 năm tập trung phát triển mô hình tương tác đất - kết cấu phi tuyến đa trường (THMC: Nhiệt - Thủy - Cơ - Hóa), phân tích ứng xử mỏi chu kỳ của móng trọng lực dưới tác động của bão biển dài hạn, tích hợp công nghệ quan trắc chuyển vị móng ngầm theo thời gian thực (Real-time Subsea Structural Health Monitoring) và ứng dụng kết cấu bê tông cốt sợi phi kim chịu ăn mòn biển sâu.

Kết luận

Luận án tiến sĩ của Đặng Xuân Trường đã khẳng định tầm vóc học thuật xuất sắc và giá trị ứng dụng thực tiễn to lớn thông qua các đóng góp cốt lõi:

1. Xác lập thành công lời giải giải tích giải tích cho tích phân $\int x \tan(ax) dx$ thông qua chuỗi số Bernoulli $B_{2n}$, bổ sung đóng góp mới vào Lý thuyết đàn hồi và Cơ học vật rắn biến dạng trong bán không gian vô hạn.
2. Thiết lập phương pháp tích phân số Simpson tính toán chính xác tải trọng sóng cực đại và mômen lật tác động lên trụ công trình biển trọng lực có tiết diện thay đổi trong vùng nước sâu.
3. Ứng dụng thành công lý thuyết hồi quy và quy hoạch thực nghiệm Box - Wilson để mô hình hóa và dự báo các thông số sóng biển ngẫu nhiên tại thềm lục địa phía Nam Việt Nam.
4. Lập trình thành công hệ thống phần mềm chuyên dụng trên Visual Basic, tự động hóa quá trình tính toán và thiết lập bản đồ trường đồng mức chuyển vị của đất nền dưới đáy móng CTBTL.
5. Kiểm chứng toàn diện ứng xử cơ lý của 10 lớp đất nền tại Mỏ Đại Hùng (Bể Nam Côn Sơn) ở độ sâu $110\text{ m}$, chứng minh sự tương thích cao giữa mô hình giải tích mới và mô hình số PLAXIS với sai số dưới $3%$.
6. Đặt nền móng khoa học vững chắc cho việc thiết kế, xây dựng các công trình biển trọng lực dầu khí, trụ điện gió ngoài khơi và đặc biệt là công trình Ụ chìm quốc phòng bảo vệ chủ quyền Biển Đông của Việt Nam.