The solvation structure of 1,3-butanediol (1,3-BD) in aqueous binary solvents of acetonitrile (AN), 1,4-dioxane (DIO), and dimethyl sulfoxide (DMSO) at various mole fractions of organic solvent xOS has been clarified by means of infrared (IR) and 1H and 13C NMR. The change in the wavenumber of O-H stretching vibration of 1,3-BD in the three systems suggested that water molecules which are initially hydrogen-bonded with the 1,3-BD hydroxyl groups in the water solvent (xOS = 0) are more significantly replaced by organic solvent molecules in the order of DMSO ≫ DIO > AN. This agrees with the order of the electron donicities of the organic solvents. The 1H and 13C chemical shifts of 1,3-BD also revealed the most remarkable replacement of water molecules on the hydroxyl groups by DMSO. In contrast to the DMSO system, the O-H vibration band of 1,3-BD in the AN and DIO systems suggested the formation of the intramolecular hydrogen bond between the two hydroxyl groups of 1,3-BD above xOS = ∼0.9. To further evaluate the intramolecular hydrogen bonding of 1,3-BD in AN-water binary solvents, molecular dynamics (MD) simulations and NMR experiments for spin-lattice relaxation times T1 and 1H-1H nuclear Overhauser effect (NOE) were conducted on 1,3-BD in the AN system. These results showed the intramolecular hydrogen bond within 1,3-BD in the AN-water binary solvents in the high AN mole fraction range of xAN > 0.9. Especially, the pair correlation functions g(r) of the OH-O interactions of 1,3-BD obtained from the MD simulations indicated that the intramolecular hydrogen bond remarkably increases in the AN solvent as the xAN rises to the unity.