Báo cáo nghiên cứu khoa học " Đánh giá tác động công trình đến bức tranh thủy động lực khu vực cửa sông ven bờ Bến Hải, Quảng Trị "

Chia sẻ: Nguyen Nhi | Ngày: | Loại File: PDF | Số trang:0

Thêm vào BST

Báo xấu

98
lượt xem 18
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Bài báo này giới thiệu một số kết quả nghiên cứu đánh giá tác động của các công trình đến các trường thủy động lực vùng cửa sông ven biển Cửa Tùng, sông Bến Hải tỉnh Quảng Trị bằng bộ mô hình MIKE 21 với các mô đun HD, SW và ST, từ đó phân tích nguyên nhân dẫn đến các

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: Báo cáo nghiên cứu khoa học " Đánh giá tác động công trình đến bức tranh thủy động lực khu vực cửa sông ven bờ Bến Hải, Quảng Trị "

Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 Đánh giá tác động công trình đến bức tranh thủy động lực khu vực cửa sông ven bờ Bến Hải, Quảng Trị Nguyễn Thọ Sáo1, Trần Ngọc Anh1,*, Nguyễn Thanh Sơn1, Đào Văn Giang2 1 Khoa Khí tượng Thủy văn và Hải dương học, Trường Đại học Khoa học Tự nhiên, ĐHQGHN, 334 Nguyễn Trãi, Hà Nội, Việt Nam 2 Trung tâm Hải văn, Tổng cục biển và hải đảo, Bộ Tài nguyên và Môi trường Nhận ngày 11 tháng 8 năm 2010 Tóm tắt. Bài báo này giới thiệu một số kết quả nghiên cứu đánh giá tác động của các công trình đến các trường thủy động lực vùng cửa sông ven biển Cửa Tùng, sông Bến Hải tỉnh Quảng Trị bằng bộ mô hình MIKE 21 với các mô đun HD, SW và ST, từ đó phân tích nguyên nhân dẫn đến các hiện tượng bồi xói bất thường trong khu vực, làm tiền đề cho việc đề xuất và quy hoạch các công trình chỉnh trị phục vụ phát triển bền vững kinh tế xã hội và môi trường. Từ khóa: MIKE, vùng cửa sông ven biển, Cửa Tùng, công trình chỉnh trị sông 1. Đặt vấn đề∗ công trình với các hiện tượng bồi xói bất thường ở các bãi biển trong khu vực. Khi có Trong những năm gần đây các hoạt động mặt các công trình như cầu, đê chắn cát, kinh tế, du lịch biển phát triển mạnh, nhiều các cảng,… trường thủy động lực sẽ có những biến công trình được xây dựng ở vùng cửa sông ven đổi nhất định và có tác động đến các quá trình biển như cầu, cảng, kè chắn sóng... đã phần nào vận chuyển bùn cát ven bờ do vậy có thể gây tác động đến các yếu tố thủy động lực trong ảnh hưởng đến hình thái vùng cửa sông ven vùng. Các tác động này có thể theo hướng tích biển. cực như chống bồi lấp luồng tàu, bảo vệ các Bài báo này giới thiệu một số kết quả vùng xung yếu hoặc cũng có thể tiêu cực như nghiên cứu đánh giá tác động của các công gây xói lở cục bộ, làm hư hại hoặc phá hủy các trình đến các trường thủy động lực vùng cửa công trình lân cận. Vùng cửa sông ven biển Cửa sông ven biển Cửa Tùng, sông Bến Hải tỉnh Tùng sông Bến Hải tỉnh Quảng Trị gần đây có Quảng Trị bằng công cụ mô hình toán, từ đó sự thay đổi mạnh mẽ về hình thái [1]. Một trong phân tích nguyên nhân dẫn đến các hiện tượng những nguyên nhân dẫn đến sự thay đổi đó có bồi xói bất thường trong khu vực, làm tiền đề thể là các công trình mới xây dựng ở khu vực cho việc đề xuất và quy hoạch các công trình này do có sự trùng hợp về thời điểm xuất hiện chỉnh trị phục vụ phát triển bền vững kinh tế xã _______ hội và môi trường kinh tế xã hội và môi trường. ∗ Tác giả liên hệ. ĐT: 84-4-38584943. E-mail: anhtnvnu.edu.vn 435
436 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 chảy sông và các tai biến lũ lụt; gió mùa đông 2. Khái quát khu vực nghiên cứu bắc, tây nam và bão, nền địa chất phức tạp trên Bãi biển Cửa Tùng, tỉnh Quảng Trị từng là khối bazan cùng với các chu kỳ biến đổi khí một trong những bãi tắm đẹp, thu hút nhiều du hậu toàn cầu và các hoạt động kinh tế xã hội khách. Hiện nay bãi biển phía bắc Cửa Tùng bị đang diễn ra mạnh mẽ trong khu vực. Những xói lở mạnh, chỉ còn một khoảng không gian năm gần đây, bãi tắm phía bắc Cửa Tùng ngày nhỏ. Về tự nhiên, Cửa Tùng nằm ở một vị trí càng bị thu hẹp về không gian do sự xâm thực địa lý phức tạp (hình 1), chịu ảnh hưởng của ngày càng gia tăng cả về quy mô lẫn cường độ, các điều kiện tự nhiên như: sóng biển, thủy dẫn tới các tổn thất về du lịch [2]. triều, hải lưu, nước dâng và dòng bùn cát; dòng Quang Tri Province Hình 1. Vị trí khu vực nghiên cứu. Nhằm đáp ứng các nhu cầu về giao thông, Cầu Tùng Luật được xây dựng vào năm giao thương giữa hai bờ Nam – Bắc sông Bến 2004 ngay nơi dòng sông gặp biển, nối hai Hải cũng như các hoạt động đánh bắt thủy sản, huyện huyện Gio Linh và Vĩnh Linh. Cầu có giao thương với huyện đảo Cồn Cỏ, thời gian thiết kế dài 461 m, rộng 9m, tải trọng H30- gần đây đã có một loạt công trình được xây XB80, khổ thông thuyền 50m, tĩnh không 8,5 m dựng (hình 2) mà tiêu biểu là cầu Tùng Luật, với 4 trụ cầu có kích thước mỗi trụ: 10,5mx 7m. cảng cá Cửa Tùng và đê chắn cát mố phía nam Kè Cửa Tùng bắt đầu được xây dựng từ cầu Tùng Luật (gọi tắt là kè). năm 2004, nằm phía bờ nam Cửa Tùng, bờ kè được xây dựng có chiều dài 430 m, cao 1,5 m,
437 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 rộng 6 m với kết cấu bằng đá hộc và cốt thép vực nghiên cứu đã được thu thập và cập nhật vươn dài ra biển nhằm mục đích chắn sóng, đến năm 2008 cùng các nguồn số liệu khác như chắn cát, giảm xói mòn trụ cầu Tùng Luật cũng tài liệu khảo sát địa hình năm 2000 và tài liệu như ngăn chặn bồi lấp cửa phục vụ giao thông phân tích cấp hạt của TEDI, ảnh hàng không thủy và an toàn hàng hải cho cảng cá Cửa Tùng. chụp năm 2003, ảnh vệ tinh Google... Cảng cá Cửa Tùng bắt đầu được xây dựng Bên cạnh đó nhằm bổ sung số liệu, hai đợt từ năm 2004, khi Sở Thuỷ sản Quảng Trị (nay khảo sát và đo đạc đã được thực hiện vào tháng thuộc Sở NN-PTNT Quảng Trị) thực hiện dự án 8/2009 và 4/2010 trong đó tập trung vào đo đạc khu neo đậu trú bão và hậu cần nghề cá với địa hình khu vực nghiên cứu phần dưới nước, mục tiêu là cung cấp bến neo đậu, tránh trú bão, phần trên cạn, dòng chảy, sóng, gió, trầm tích hậu cần nghề cá và hiện đại hoá khu sản xuất đáy và lơ lửng kết hợp với các điều tra thu thập nghề cá ven biển. Trước đây, tại vị trí cảng cá thông tin xã hội [1]. Số liệu khảo sát sau đó Cửa Tùng là một eo biển kín gió, được một cồn được xử lý và đưa về các dạng tương thích phục cát lớn nằm phía ngoài che chắn sóng biển. Khi vụ công tác thiết lập mô hình ở bước tiếp theo thực hiện dự án, hơn 200.000 m3 cát ở cồn này K1 bị múc đi đổ vào san lấp eo biển tạo thành một Cảng cá bãi cát bằng phẳng chính là mặt bằng cảng cá hiện nay. K2 Nhằm mục đích đánh giá tác động của các công trình cầu – kè – cảng đến bức tranh thủy động lực trong vùng, khu vực nghiên cứu được xác định gồm hai miền: miền trong sông từ cầu Kè Hiền Lương đến cửa và miền biển ven bờ (xem Cầu Tùng Luật hình 3 3. Phương pháp tiếp cận Hình 2. Sơ đồ bố trí công trình khu vực nghiên cứu. Để đánh giá tác động của công trình hoặc hệ thống công trình đến chế độ thủy động lực có nhiều phương pháp khác nhau, tuy nhiên sử dụng (các) mô hình thủy động lực là phương pháp rẻ và tin cậy đồng thời có thể cung cấp các giải pháp thay thế cũng như các dự báo cho tương tai. Hỗ trợ cho mô hình thủy động lực là các biện pháp điều tra, thu thập và xử lý số liệu để tăng cường độ tin cậy của các đầu vào yêu cầu cho mô hình cũng như phân tích nhằm làm sáng tỏ bức tranh thủy động lực trong khu vực trước và sau khi có công trình, từ đó đánh giá sơ bộ các luận điểm và sẽ được minh chứng lại Hình 3. Lưới phần tử hữu hạn dùng trong mô hình bằng các công cụ mô phỏng trong các bước tiếp MIKE 21 FM. theo. Số liệu khí tượng, thủy hải văn trong khu
438 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 100 m2, góc nhỏ nhất 23º; vùng tính toán được 4. Kết quả mô phỏng rời rạc hóa thành 10767 phần tử với 5717 nút Trên thế giới và trong nước hiện có nhiều lưới, độ phân giải thô nhất ở vùng ngoài khơi mô hình thủy động lực đang được áp dụng cho khoảng 400m, mịn nhất vùng gần bờ và các nhiều mục đích khai thác khác nhau như nghiên công trình khoảng 20m. Mỗi kịch bản tính toán cứu, quy hoạch và thiết kế hệ thống công về sự tồn tại hay không tồn tại công trình sử trình... tiêu biểu có thể kể đến SORBEK (Hà dụng một lưới tính khác nhau và hình 3 minh Lan), MIKE (Đan Mạch),…tuy nhiên, mỗi mô họa một ví dụ về lưới tính với trường hợp hiện hình đều có những ưu nhược điểm riêng và cho trạng (có cả cụm công trình cầu - kè - cảng). đến nay vẫn chưa có một đánh giá toàn diện và chi tiết về khả năng áp dụng trong thực tế của * Điều kiện biên các mô hình nói trên. Trong nghiên cứu này, Ba mô đun HD, SW và ST có thể tính toán với mục tiêu mô phỏng và tính toán ảnh hưởng riêng rẽ theo từng mục đích, tuy nhiên trong bài của các công trình lên trường thủy động lực này sử dụng kết quả mô hình FM liên hợp, tức vùng cửa sông ven biển Cửa Tùng - Quảng Trị, là tính toán trường thủy động lực có xét đến tác bộ mô hình MIKE 21 đã được lựa chọn do đáp động sóng, đồng thời với mô hình vận chuyển ứng được những tiêu chí: a) Là bộ phần mềm trầm tích và biến đổi đáy. Như vậy mô hình tích hợp đa tính năng (tính toán trường sóng, sóng được tính trước tiên với 3 hướng chính dòng chảy, vận chuyển trầm tích, diễn biến địa NE, E và SE. Để tiết kiệm thời gian tính toán, hình đáy; b) Đã được kiểm nghiệm thực tế ở sóng được coi phát triển hoàn toàn và ổn định nhiều quốc gia trên thế giới [3]; c) Giao diện trong suốt thời gian tính (ứng suất bức xạ không thân thiện, dễ sử dụng và tương thích với nhiều đổi theo thời gian). Đối với tính toán sóng, sử phần mềm GIS khác. dụng điều kiện biên là độ cao, chu kỳ và hướng 4.1. Thiết lập mô hình sóng tại biên phía biển, tại biên phía bắc và nam là điều kiện đối xứng. Trong tính toán mực * Miền tính và lưới tính nước và dòng chảy sử dụng điều kiện biên tại Miền tính được xác định từ khoảng vĩ độ Bến Hải là mực nước hoặc lưu lượng, tại phía 16º58’ N - 17º4’ N , kinh độ 107º3’ E - 107º10’ biển là các hằng số điều hòa thủy triều của 8 E, từ cầu Hiền Lương đến cầu Tùng Luật và từ sóng theo bản đồ đồng triều toàn cầu.. Điều cầu Tùng Luật ra biển cách bờ khoảng 3,5 km kiện biên đối với trầm tích lơ lửng là cân bằng (hình 3). Để tính toán các trường thủy động lực nồng độ, tham số trầm tích đáy lựa chọn là d50 và sau đó đánh giá tác động công trình đến bức = 0.27mm, độ chọn lọc cát là 1.4 [1]. tranh thủy động lực khu vực cửa sông ven bờ 4.2. Hiệu chỉnh và kiểm định mô hình Bến Hải, Quảng Trị trên cơ sở các tài liệu địa hình đã có, nghiên cứu này đã lựa chọn lưới Mô hình được hiệu chỉnh và kiểm định với phần tử hữu hạn với độ phân giải thay đổi tăng tài liệu đo đạc dòng chảy và mực nước tại hai dần từ ngoài biển vào sát bờ và cao nhất xung điểm trong vùng nghiên cứu (1 điểm trong vùng quanh các công trình. Trong toàn miền diện tích biển ven bờ K1 và 1 điểm trong vùng cửa sông của phần tử lớn nhất là 100000 m2, nhỏ nhất là K2, gần khu cảng cá) với chuỗi số liệu từ ngày
439 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 12-18/8/2009 và từ ngày 21-28/4/2010. Kết quả tính toán mô phỏng khá tốt đối với mực nước, nhưng chưa tốt lắm đối với vận tốc. Mặc dù vậy xu thế, dáng điệu và sự trùng pha tương đối phù hợp với thực đo cho phép đánh giá mô hình được thiết lập với bộ thông số đã hiệu chỉnh là đáng tin cậy và có thể sử dụng cho tính toán các kịch bản tiếp theo. Các kết quả tính toán dưới đây ứng với trường hợp 12-18/8/2009, các hình Hình 5. So sánh lưu tốc đo và tính toán tại K2. 4, 5, 6. 0.6 0.4 5. Đánh giá tác động các công trình 0.2 Trước hết cần nhận định rằng trong khu vực H (m ) 0 nghiên cứu, trường thủy động lực do sông và -0 .2 thủy triều không đáng kể so với trường sóng. Biên độ thủy triều chỉ khoảng 45cm, dòng chảy -0 .4 Obs. Cal. sông trung bình khoảng 20m3/s, trong khi độ -0 .6 12/8/2009 12:00 AM 14/8/2009 12:00 AM 16/8/2009 12:00 AM 18/8/2009 12:00 AM cao sóng có nghĩa xấp xỉ 3m tại ranh giới phía Ti me đông khu vực ngiên cứu Hình 6. So sánh mực nước đo và tính toán tại K2. Mặc dù điều kiện sóng tại Cồn Cỏ rất đa dạng, tuy nhiên Cồn Cỏ ở xa khu vực nghiên cứu khoảng 34km, do đó cần lấy sóng Cửa Tùng với ba hướng sóng cơ bản là NE, E và SE. Hình hoa sóng tại Cồn Cỏ và Cửa Tùng dưới đây cho thấy điều đó, lấy tháng 9-12 làm đại diện (hình 7), các tháng khác xu thế cũng tương tự Để làm cơ sở so sánh và đánh giá tác động của các công trình, ngoài việc tính toán với điều Hình 4. So sánh lưu tốc đo và tính toán tại K1. kiện hiện trạng cụm công trình (cầu-kè-cảng) để kiểm định mô hình ở trên, cần có kịch bản nền ứng với điều kiện tự nhiên khi chưa có các công trình.
440 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 Hình 7. Hoa sóng tại Cồn Cỏ và Cửa Tùng trong các tháng 9-12. Trước đây trong tình trạng tự nhiên, khu thượng lưu, dẫn đến thay đổi rõ nét về dòng vực cảng cá là một eo biển kín gió, được một chảy trong khu vực, đặc biệt, trục động lực cửa cồn cát lớn nằm phía ngoài che chắn sóng biển. sông lệch dần sang phải nếu nhìn ra biển. Cửa Khi thực hiện dự án, con người đã có những tác sông được mở rộng nên tạo cơ hội thuận lợi hơn động đáng kể đến địa hình đáy khu vực nghiên cho các trao đổi giữa hai thủy vực sông – biển cứu do vậy nghiên cứu này đã sử dụng tài liệu đồng thời tạo nên các xoáy cục bộ vùng lân cận địa hình 2000 [4,5] để làm kịch bản nền. Ngoài cửa làm tăng khả năng lắng đọng các hạt trầm ra, trên cơ sở các tài liệu về địa hình hiện trạng tích lơ lửng. và kích thước, hình dạng các công trình, thêm 3 kịch bản tính toán đươc thiết lập bao gồm: a) Kịch bản 1: địa hình 2010 + cảng cá; b) Kịch bản 2: địa hình 2010 + cảng + cầu; c) Kịch bản 3: địa hình 2010 + cầu + kè + cảng. Mỗi kịch bản tính toán được thực hiện đầy đủ cả 3 mô đun HD, SW và ST, ứng với các hướng sóng thịnh hành là Đông Bắc, Đông và Đông Nam. Tuy nhiên trong giới hạn của bài báo chỉ trình bày minh họa kết quả tính toán trường dòng chảy trong khu vực nghiên cứu ứng với trường hợp sóng Đông Nam (hình 8-11). Việc xây dựng cảng cá với các tác động trực tiếp đến bar cát chắn cửa sông chắc chắn sẽ mang lại nhiều thay đổi bức tranh thủy động lực. Các phân tích dựa trên kết quả tính toán với điều kiện địa hình 2000 và 2010 cho thấy, cảng cá đã làm tăng đáng kể bề rộng mặt cắt ngang Hình 8. Trường dòng chảy sóng Đông Nam trường cửa, nhưng cũng làm thu hẹp một phần phía hợp chưa có công trình.
441 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 Hình 9. Trường dòng chảy sóng Đông Nam kịch bản 1. Hình 11. Trường dòng chảy trong sóng Đông Nam kịch bản 3. Mặt khác, sau khi xây dựng cầu đã làm cho mặt cắt sông bị thu hẹp lai và kết hợp với các trụ cầu đã làm thay đổi hướng của dòng chảy và tạo lên các xoáy cục bộ phía sau trụ cầu. Mặc dầu vậy không quan sát thấy sự gia tăng đáng kể về vận tốc giữa các trụ cầu do ảnh hưởng của dòng chảy sông trong khu vực không lớn và do vậy các trụ cầu cũng không gây hiện tượng bồi lắng hay xói lở sau cầu. Ngược lại đối với trường sóng, sự xuất hiện các trụ cầu và mố cầu đã cản trở sự lan truyền của sóng vào phía trong sông, sóng gần như tắt hẳn ngay phía thượng lưu cầu Tùng Luật. Công trình kè (đê chắn cát) phía nam tác động mạnh mẽ đến cả hướng sóng và dòng chảy. Do hướng vuông góc với bờ, kè phía Nam Hình 10. Trường dòng chảy trong sóng Đông Nam đã có tác dụng chắn toàn bộ sóng Đông Nam, kịch bản 2. tạo nên vùng khuất sóng phía trong cửa. Đồng thời với việc chắn sóng, kè cũng có tác dụng ngăn dòng chảy do sóng dọc bờ hướng từ phía Nam lên phía Bắc và kết quả phân tích cho thấy kè đã giữ lại một lượng lớn cát ở phía Nam và
442 N.T. Sáo và nnk. / Tạp chí Khoa học ĐHQGHN, Khoa học Tự nhiên và Công nghệ 26, Số 3S (2010) 435‐442 do đó đã hạn chế hiện tượng bồi lấp luồng trong cảng cá là khá rõ ràng. Đây chính là những mùa sóng Đông Nam. Tuy vậy, kè cũng làm nguyên nhân dẫn tới những biến đổi bất thường xuất hiện một số xoáy cục bộ và làm tán xạ về mặt hình thái trong khu vực. Nghiên cứu này sóng khiến cho sóng đi thẳng vào mố cầu phía đặt cơ sở cho việc nghiên cứu đề xuất các biện Nam, có thể gây hiện tượng sạt lở mố cầu trong pháp giảm thiểu tác động bất lợi bảo vệ các mùa sóng Đông Bắc. Mặt khác, kè cũng có tác nguồn tài nguyên trong khu vực phục vụ phát dụng đẩy dòng chảy sóng dọc bờ ra xa hơn, triển kinh tế xã hội bền vững. khiến cho nguồn bùn cát phía Nam ít có cơ hội tiếp cận cửa và bãi biển phía bắc Cửa Tùng. Đó Tài liệu tham khảo cũng có thể là một phần nguyên nhân của hiện tượng thiếu hụt nguồn trầm tích bãi biển phía [1] Nguyễn Thọ Sáo, Nguyễn Minh Huấn, Ngô Chí bắc, làm mất cân bằng cán cân bùn cát và gây Tuấn, Đặng Đình Khá (2010), Biến động trầm tích và diễn biến hình thái khu vực cửa sông ven hiện tượng xói lở bãi biển bờ Cửa Tùng, Quảng Trị, Tạp chí Khoa học, ĐHQGHN, (đã nộp bản thảo tháng 8 năm 2010). [2] Cổng thông tin điện tử, UBND tỉnh Quảng Trị Kết luận www.quangtri.gov.vn. [3] Danish Hydraulics Institute (DHI). “MIKE21. Qua phân tích tài liệu và kết quả tính toán, User’s Mannual”. 2007 trường thủy động lực trong khu vực đã có [4] Báo cáo khảo sát địa hình-thủy văn. Công ty Tư những thay đổi mạnh mẽ khi có sự xuất hiện vấn Xây dựng Cảng đường thủy, Hà Nội, 2000. của các công trình. Tuy nhiên, các công trình [5] Báo cáo khảo sát địa chất công trình, Công ty Tư khác nhau có những tác động khác nhau trong vấn Xây dựng Cảng đường thủy, Hà Nội, 2000. đó tác động của công trình cầu ít rõ nét trong khi tác động của đê chắn cát và việc xây dựng Assessment of structural impacts on hydrodynamic patterns in Ben Hai river mouth and coastal zones, Quang Tri province Nguyen Tho Sao1, Tran Ngoc Anh1, Nguyen Thanh Son1, Đao Van Giang2 1 Faculty of Hydro-Meteorology & Oceanography, Hanoi University of Science, VNU, 334 Nguyen Trai, Hanoi, Vietnam 2 Marine Center, General Service for Sea and Islands, MONRE This paper describes the scientific researches of structural simpacts on hydrodynamic conditions in Ben Hai river mouth and Cua Tung beach, Quang Tri province. These hydrodynamic patterns were simulated by using MIKE 21 including hydrodynamics (HD), spectral wave (SW) and sand transport (ST) modules. The research results may take a part in analysis of the reasons leading to beach erosion and deposition in the study region. It may be served as the good information for planning and proposing the suitable training measures for social economical sustainable development.