Tổng quan nghiên cứu

Tiếng ồn đô thị đang trở thành một trong những vấn đề môi trường cấp bách nhất tại các thành phố hiện đại. Theo ước tính từ các cơ quan bảo vệ môi trường, mức áp suất âm tại các trục giao thông chính và khu công nghiệp thường xuyên chạm ngưỡng từ 80 dB đến 90 dB, vượt xa ngưỡng an toàn thính giác của con người là 70 dB và tiềm ẩn nguy cơ gây tổn hại sức khỏe nghiêm trọng khi tiến sát ngưỡng đau 130 dB. Trong môi trường không khí chuẩn ở nhiệt độ 20°C với áp suất khí quyển 101 kPa, sóng âm lan truyền với vận tốc xấp xỉ 344 m/s. Sự phân bố và lan truyền sóng âm trong không gian mở chịu sự chi phối phức tạp của hàng loạt hiện tượng vật lý như phản xạ, khúc xạ, tán xạ, giao thoa và nhiễu xạ qua các khối công trình cao tầng.

Nhằm giải quyết bài toán kiểm soát ô nhiễm tiếng ồn và hỗ trợ quy hoạch hạ tầng, luận văn thạc sĩ chuyên ngành Kỹ thuật Điện tử tại Trường Đại học Bách Khoa – Đại học Quốc gia TP. Hồ Chí Minh (hoàn thành vào tháng 06 năm 2014) đã tập trung nghiên cứu đề tài: "Áp dụng phương pháp ma trận đường dây truyền dẫn cho truyền sóng âm trong không gian mở". Mục tiêu trọng tâm của nghiên cứu là xây dựng khung mô phỏng số trong miền thời gian dựa trên phương pháp ma trận đường dây truyền dẫn (Transmission Line Matrix – TLM), kết hợp thiết lập các điều kiện biên hấp thụ không gian vô hạn (Perfectly Matched Layer – PML) và biên trở kháng phức đại diện cho bề mặt mặt đất.

Phạm vi nghiên cứu tập trung vào mô hình không gian mở đô thị 2D với cấu trúc các tòa nhà cao tầng và các dạng tường chắn âm thanh khác nhau. Nghiên cứu mang ý nghĩa khoa học và thực tiễn sâu sắc: không chỉ cung cấp một công cụ tính toán số có độ chính xác cao và tiết kiệm tài nguyên bộ nhớ máy tính hơn so với các phương pháp truyền thống, mà còn định lượng trực quan sự phân bố mức áp suất âm (Sound Pressure Level – SPL) với ngưỡng tham chiếu $p_0 = 2 \times 10^{-5}\text{ Pa}$, tạo tiền đề vững chắc cho việc thiết kế kiến trúc và giảm thiểu tiếng ồn đô thị tại Việt Nam.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự tích hợp giữa lý thuyết âm học cổ điển và phương pháp giải tích số trong kỹ thuật điện tử:

  • Lý thuyết sóng âm và các phương trình bảo toàn: Hệ thống phương trình sóng âm bậc hai trong không gian ba chiều được thiết lập từ định luật bảo toàn khối lượng (phương trình liên tục: $\frac{\partial p}{\partial t} + \rho c^2 \nabla \cdot \mathbf{u} = 0$) và định luật bảo toàn động lượng ($\frac{\partial \mathbf{u}}{\partial t} + \frac{1}{\rho}\nabla p = 0$), trong đó $p$ là áp suất âm, $\rho$ là mật độ môi trường và $\mathbf{u}$ là vận tốc dao động hạt. Các hiện tượng đặc trưng bao gồm nguyên lý Huygens về nguồn thứ cấp, định luật phản xạ góc tới bằng góc phản xạ, định luật khúc xạ dạng Snell ($\frac{\sin \theta_i}{\sin \theta_r} = \frac{c_1}{c_2}$), hiện tượng sóng đứng do giao thoa và hiệu ứng Doppler.
  • Phương pháp ma trận đường dây truyền dẫn (TLM): Khởi nguồn từ công trình nghiên cứu của Johns và Beurle trong điện từ trường, sau đó được phát triển sang lĩnh vực âm học bởi Saleh, Blanchfield và Kagawa. Không gian truyền sóng được rời rạc hóa thành lưới các nút nối tiếp (series node) hoặc nút song song (shunt node). Các đại lượng cơ học được ánh xạ tương đương sang mạch điện: áp suất âm $p$ tương ứng với điện áp $V$, vận tốc hạt $u$ tương ứng với dòng điện $I$, và trở kháng âm đặc trưng môi trường $Z_0 = \rho_0 c_0 \approx 415\text{ N}\cdot\text{s/m}^3$ tương ứng với trở kháng đường dây.
  • Mô hình hóa điều kiện biên: Thiết lập lớp hấp thụ PML theo mô hình của Gwenael và Picaut cũng như De Cogan nhằm triệt tiêu phản xạ giả tại biên vô hạn của bầu khí quyển. Đồng thời, mô hình hóa biên trở kháng mặt đất phức dựa trên lý thuyết vật liệu xốp của Miki và Zwikker - Kosten (ZW), được tối ưu hóa tham số thông qua bài toán quy hoạch toàn phương tuần tự (Sequential Quadratic Programming – SQP) kế thừa từ Cotte.

Các khái niệm cốt lõi bao gồm:

  1. Mức áp suất âm (SPL): Thước đo logarit mức độ mạnh yếu của âm thanh tính bằng decibel: $\text{SPL} = 20 \log_{10}(p / p_0)$ với áp suất chuẩn $p_0 = 2 \times 10^{-5}\text{ Pa}$.
  2. Trở kháng âm học đặc trưng ($Z_0$): Tỷ số giữa áp suất âm và vận tốc hạt dao động trong môi trường truyền.
  3. Ma trận tán xạ ($S$) và ma trận liên kết ($C$): Hai ma trận toán học cốt lõi điều khiển quá trình phản xạ và truyền sóng giữa các cổng lân cận tại mỗi bước thời gian $k$.
  4. Lớp hấp thụ phối hợp hoàn hảo (PML): Cấu trúc biên nhân tạo làm suy hao sóng âm tiệm cận về 0 mà không tạo ra sóng phản xạ ngược vào miền tính toán.

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

Nghiên cứu sử dụng phương pháp mô phỏng số trên nền tảng MATLAB kết hợp đối sánh lý thuyết giải tích:

  • Nguồn dữ liệu và thiết lập tham số: Toàn bộ dữ liệu được tạo ra thông qua các bài toán mô phỏng lan truyền sóng xung Gaussian và sóng điều hòa hình sin ở các dải tần số 100 Hz, 500 Hz và 1000 Hz, tương ứng với dải tần nhạy cảm của thính giác con người (khoảng 20 Hz đến 20 kHz).
  • Cỡ mẫu và phân chia lưới không gian - thời gian: Miền tính toán 2D được rời rạc hóa với quy mô hàng chục nghìn nút mạng. Khoảng chia không gian $\Delta l$ được chọn nghiêm ngặt theo tiêu chí ổn định: $\Delta l \le \lambda / 10$ (với $\lambda$ là bước sóng ngắn nhất). Khoảng chia thời gian $\Delta t$ được khóa chặt với bước không gian theo vận tốc truyền sóng qua hệ thức $\Delta t = \Delta l / (c_0 \sqrt{2})$, đảm bảo thỏa mãn điều kiện ổn định Courant - Friedrichs - Lewy (CFL).
  • Lý do lựa chọn phương pháp phân tích: So với phương pháp phần tử hữu hạn (FEM) đòi hỏi giải hệ phương trình ma trận lớn hay phương pháp sai phân hữu hạn miền thời gian (FDTD) dễ mất liên kết trực tiếp giữa các bước không - thời gian, phương pháp TLM mang bản chất vật lý rõ ràng, theo dõi trực tiếp các xung sóng phản xạ và tán xạ tại từng nút theo thời gian thực, tiết kiệm đáng kể bộ nhớ RAM và thời gian tính toán của CPU.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng và khảo sát bằng thuật toán TLM đã mang lại 4 phát hiện quan trọng có giá trị thực tiễn cao:

  1. Hiệu năng vượt trội của lớp biên hấp thụ PML: Mô hình biên hấp thụ PML theo tiếp cận của Gwenael đã triệt tiêu sóng phản xạ tại biên khí quyển với sai số phản xạ dư chỉ dưới 1,2% ở tần số 1000 Hz. So với phương pháp của De Cogan (có mức sai số dao động từ 4,2% đến 5,1%), mô hình này giúp giảm độ dày lớp đệm tính toán đi khoảng 30% nhưng vẫn đảm bảo độ chính xác tương đương môi trường không gian tự do vô hạn.
  2. Độ chính xác cao trong xấp xỉ trở kháng phức mặt đất: Thuật toán tối ưu hóa SQP áp dụng cho mô hình Miki và Zwikker - Kosten đã xác định chính xác các hệ số xấp xỉ phân số hữu tỉ trong miền thời gian. Độ lệch biên độ trở kháng giữa mô hình xấp xỉ và lý thuyết thực nghiệm luôn duy trì ở mức dưới 2,5% trên toàn bộ dải tần số khảo sát từ 100 Hz đến 1000 Hz.
  3. Hiệu ứng suy giảm âm thanh của tường chắn chữ L vượt trội hơn tường chữ I: Khảo sát cho thấy bức tường chắn dạng chữ L tạo ra vùng bóng tối âm học (shadow zone) rộng hơn đáng kể phía sau vật cản. Mức áp suất âm SPL tại các điểm thu sau tường chữ L giảm từ 12 dB đến 16 dB (tương ứng mức suy giảm năng lượng sóng âm trên 50%) so với khi không có tường, trong khi tường thẳng dạng chữ I chỉ giúp giảm từ 6 dB đến 9 dB.
  4. Sự hình thành các nút và bụng sóng đứng trong hẻm đô thị (Street Canyon): Khi sóng âm truyền giữa các khối nhà cao tầng có bề mặt phản xạ cứng, hiện tượng giao thoa giữa sóng trực tiếp và sóng phản xạ tạo nên các vân sóng đứng rõ rệt. Tại thời điểm mô phỏng 60 ms, mức áp suất âm SPL tại các điểm bụng sóng ghi nhận mức tăng cục bộ từ 4 dB đến 7 dB so với trường sóng tự do.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự phân bố áp suất âm phức tạp trong đô thị xuất phát từ sự chồng chập pha của các tia sóng phản xạ từ mặt đất và mặt dựng công trình. Theo nguyên lý Huygens, mép trên của các tòa nhà và đỉnh tường chắn đóng vai trò như các nguồn phát thứ cấp, bẻ cong đường truyền sóng xuống vùng không gian phía sau thông qua hiện tượng nhiễu xạ.

Dữ liệu mô phỏng trong luận văn được trực quan hóa rất hiệu quả thông qua biểu đồ phân bố đường đẳng mức áp suất âm 2D (SPL Contour Map) trên toàn miền khảo sát và các bảng đối sánh dạng sóng áp suất tức thời theo thời gian. Các biểu đồ này thể hiện rõ nét biên độ dao động tại các thời điểm $t = 20\text{ ms}$, $40\text{ ms}$ và $60\text{ ms}$. Khi so sánh với các nghiên cứu trước đây của Cotte hay Hofmann, kết quả của luận văn đã chứng minh tính ưu việt của phương pháp TLM trong việc tái hiện đồng thời cả tác động hấp thụ của mặt đất xốp lẫn cấu trúc hình học đa giác của các vật cản đô thị mà không làm phát sinh sai số tích lũy lớn.

Đề xuất và khuyến nghị

Dựa trên các kết quả đạt được, luận văn đưa ra 4 nhóm giải pháp kỹ thuật và quy hoạch cụ thể:

  1. Tích hợp mô hình TLM vào phần mềm dự báo tiếng ồn đô thị: Đề nghị các Viện Quy hoạch Đô thị và Sở Xây dựng ứng dụng thuật toán TLM 2D/3D vào các công cụ phần mềm chuyên dụng nhằm đánh giá tác động tiếng ồn. Mục tiêu là đưa sai số dự báo âm học đô thị xuống dưới 5%, lộ trình triển khai thử nghiệm trong vòng 12 đến 24 tháng.
  2. Quy chuẩn hóa thiết kế tường chắn âm thanh dạng chữ L tại các tuyến cao tốc: Khuyến nghị Bộ Giao thông Vận tải và các ban quản lý dự án ưu tiên lắp đặt tường chắn âm dạng chữ L hoặc tường có gờ cong thay vì tường thẳng truyền thống dọc các tuyến đường vành đai. Mục tiêu là triệt tiêu từ 10 dB đến 15 dB tiếng ồn xâm nhập vào các khu dân cư lân cận trong giai đoạn quy hoạch 2024 - 2026.
  3. Ứng dụng vật liệu hấp thụ âm cho bề mặt mặt đất và mặt dựng tòa nhà: Khuyến cáo các nhà thầu và kiến trúc sư lựa chọn vật liệu lát vỉa hè xốp và vật liệu phủ tường có trở kháng âm tối ưu (phù hợp với mô hình Zwikker - Kosten). Mục tiêu là nâng cao hiệu suất tiêu âm bề mặt lên khoảng 35% đến 40%, giảm thiểu hiện tượng phản xạ cộng hưởng trong các hẻm phố hẹp.
  4. Xây dựng bản đồ phân vùng ô nhiễm tiếng ồn tự động: Đề xuất Sở Tài nguyên và Môi trường tại các thành phố lớn thiết lập hệ thống bản đồ số hóa giám sát mức SPL theo thời gian thực, chủ động khoanh vùng cảnh báo các khu vực có tiếng ồn vượt mức 80 dB trong thời gian 18 tháng tới.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo chuyên sâu và hữu ích cho 4 nhóm đối tượng chính:

  1. Kỹ sư âm học và chuyên gia viễn thông - điện tử: Nắm bắt phương pháp chuyển đổi bài toán truyền sóng âm sang mô hình mạch điện tương đương, ứng dụng ma trận tán xạ và lý thuyết đường dây truyền dẫn để phát triển các thuật toán xử lý tín hiệu và mô phỏng số.
  2. Kiến trúc sư và chuyên gia quy hoạch đô thị: Ứng dụng các kết quả khảo sát về hình học tường chắn và khoảng cách công trình để tối ưu hóa vị trí xây dựng khu dân cư, bệnh viện, trường học tránh xa các nguồn phát tiếng ồn lớn.
  3. Giảng viên, nghiên cứu sinh và sinh viên ngành Kỹ thuật: Sử dụng làm giáo trình chuyên đề và tài liệu nghiên cứu chuyên sâu về các phương pháp số trong miền thời gian (TLM, FDTD, FEM) cũng như kỹ thuật thiết lập biên PML và trở kháng phức trên MATLAB.
  4. Cơ quan quản lý nhà nước về tài nguyên và môi trường: Tham khảo các số liệu và mô hình phân bố áp suất âm để làm căn cứ khoa học xây dựng bộ quy chuẩn kỹ thuật quốc gia về giới hạn tiếng ồn trong đô thị.

Câu hỏi thường gặp

1. Phương pháp TLM có ưu điểm gì vượt trội so với FDTD trong mô phỏng âm học không gian mở?
Phương pháp TLM liên kết trực tiếp giữa bước không gian $\Delta l$ và bước thời gian $\Delta t$ thông qua vận tốc truyền sóng thực tế $c_0 = 344\text{ m/s}$, giúp phản ánh chính xác bản chất vật lý của sóng. Ngoài ra, TLM sử dụng ma trận tán xạ để tính toán trực tiếp các xung phản xạ tại nút, giúp giảm thiểu đáng kể dung lượng bộ nhớ RAM khi mô phỏng miền không gian lớn so với việc giải hệ sai phân trực tiếp của FDTD.

2. Tại sao cần phải thiết lập lớp hấp thụ PML trong bài toán mô phỏng sóng âm?
Trong mô phỏng số, máy tính chỉ có thể xử lý một miền không gian hữu hạn. Nếu không có lớp hấp thụ PML, sóng âm khi lan truyền chạm đến đường biên nhân tạo sẽ bị phản xạ ngược lại vào trong, gây sai lệch hoàn toàn kết quả giao thoa trường âm. Lớp PML của Gwenael hấp thụ triệt để năng lượng sóng tới với mức phản xạ dư dưới 1,2%, tái lập hoàn hảo môi trường bầu khí quyển vô hạn.

3. Mô hình trở kháng đất Miki và Zwikker - Kosten đóng vai trò gì trong mô phỏng?
Mặt đất thực tế không phải là bề mặt phản xạ cứng hoàn toàn mà có tính xốp và đàn hồi. Hai mô hình Miki và Zwikker - Kosten mô tả chính xác trở kháng âm phức phụ thuộc tần số của đất. Khi kết hợp với thuật toán tối ưu SQP, các mô hình này cho phép phương pháp TLM tính toán chuẩn xác năng lượng sóng âm bị hút vào lòng đất và độ lệch pha của sóng phản xạ.

4. Điều kiện chia lưới không gian và thời gian trong TLM cần tuân thủ nguyên tắc nào?
Để đảm bảo thuật toán hội tụ và triệt tiêu hiện tượng tán sắc số học, kích thước ô lưới không gian bắt buộc phải thỏa mãn điều kiện $\Delta l \le \lambda / 10$, trong đó $\lambda$ là bước sóng ứng với tần số khảo sát cao nhất (ví dụ $\lambda = 0,344\text{ m}$ tại $1000\text{ Hz}$). Bước thời gian $\Delta t$ sau đó được cố định theo công thức $\Delta t = \Delta l / (c_0 \sqrt{2})$.

5. Kết quả khảo sát tường chắn âm thanh dạng chữ L có thể ứng dụng ngay vào thực tế không?
Hoàn toàn có thể. Mô phỏng đã chứng minh phần gờ ngang của tường chữ L tạo ra hiệu ứng cản nhiễu xạ mạnh mẽ, giúp giảm thêm từ 3 dB đến 7 dB so với tường phẳng chữ I cùng chiều cao. Đây là cơ sở thực chứng để các nhà thầu giao thông tối ưu kết cấu vách chống ồn mà không cần tăng quá nhiều chiều cao công trình.

Kết luận

  • Áp dụng thành công phương pháp TLM vào âm học: Luận văn đã chứng minh tính khả thi, độ chính xác và hiệu quả vượt trội của phương pháp ma trận đường dây truyền dẫn trong việc mô phỏng lan truyền sóng âm miền thời gian.
  • Xây dựng hoàn chỉnh hệ thống điều kiện biên: Thiết lập thành công lớp biên hấp thụ PML triệt tiêu phản xạ và biên trở kháng phức xấp xỉ chính xác bề mặt đất thực tế thông qua tối ưu hóa SQP.
  • Định lượng tác động của vật cản và tường chắn đô thị: Khảo sát chi tiết sự suy giảm áp suất âm qua các dạng tường chắn chữ I và chữ L, chứng minh hiệu quả giảm thiểu tiếng ồn từ 12 dB đến 16 dB của tường dạng chữ L.
  • Làm sáng tỏ hiện tượng sóng đứng trong đô thị: Phát hiện và mô hình hóa trực quan sự gia tăng cục bộ từ 4 dB đến 7 dB mức SPL do giao thoa sóng tại các hẻm phố và chân tòa nhà cao tầng.
  • Khai mở hướng nghiên cứu mới: Đặt nền móng lý thuyết và giải thuật số cho các công trình nghiên cứu sâu hơn về kiểm soát ô nhiễm âm thanh và thiết kế âm học công trình tại Việt Nam.

Trong các bước nghiên cứu tiếp theo (lộ trình 2024 - 2026), hướng phát triển trọng tâm sẽ là mở rộng không gian mô phỏng lên mô hình 3D hoàn chỉnh, tích hợp thêm các yếu tố khí tượng thực tế như gradient nhiệt độ, khúc xạ do gió và độ ẩm môi trường. Các chuyên gia âm học, kỹ sư xây dựng và nhà quy hoạch đô thị hãy ứng dụng ngay mô hình tính toán này để xây dựng những không gian đô thị văn minh, yên tĩnh và bền vững.