A method of acquiring seismic survey data using phase encoded simultaneous sourcing technique for seismic acquisitions and using two to eight vibratory sources and one or more receivers, the method comprising: a) determining the number of vibratory sources; b) determining a number of simultaneous vibration sweeps, the vibratory sources being simultaneously activated during each of said sweeps; c) establishing a phase encoding scheme, wherein phase shifts are calculated using a Monte Carlo simulation; d) positioning said vibratory sources proximate to a predetermined location on the surface of the earth; e) energizing said vibratory sources according to said phase encoding scheme thereby vibrating the earth to produce a seismic signal; f) acquiring the reflected seismic signal data at said one or more receivers; and, g) processing said data to create a seismic survey; the method being characterized in that the phase encoding scheme being a non-orthogonal phase scheme; and step (c) involves choosing from the following pre-defined calculated phase encoding schemes indicating phase shifts, for said number of vibratory sources ("viv1, vib2, vib3, vib4, vib5, vib6, vib7, vib8") and vibration sweeps ("sweep"): sweep vib1 vib2 1 100 3 2 3 100 sweep vib1 vib2 vib3 1 148 194 29 2 194 148 194 3 29 194 148 sweep vib1 vib2 vib3 vib4 1 180 81 60 353 2 81 180 353 60 3 60 353 180 81 4 353 60 81 180 sweep vib1 vib2 vib3 vib4 1 137 294 238 245 2 347 9 65 252 3 274 200 295 27 4 38 144 239 151 5 356 290 103 17 6 232 346 159 253 sweep vib1 vib2 vib3 vib4 vib5 1 204 301 319 55 38 2 301 204 38 319 55 3 319 38 204 38 319 4 55 319 38 204 301 5 38 55 319 301 204 sweep vib1 vib2 vib3 vib4 vib5 1 277 343 27 304 194 2 194 277 343 27 304 3 304 194 277 343 27 4 27 304 194 277 343 5 343 27 304 194 277 sweep vib1 vib2 vib3 vib4 vib5 vib6 1 16 53 266 94 285 321 2 53 16 321 266 94 285 3 266 321 16 285 266 94 4 94 266 285 16 321 266 5 285 94 266 321 16 53 6 321 285 94 266 53 16 sweep vib1 vib2 vib3 vib4 vib5 vib6 1 352 277 322 155 61 302 2 302 352 277 322 155 61 3 61 302 352 277 322 155 4 155 61 302 352 277 322 5 322 155 61 302 352 277 6 277 322 155 61 302 352 sweep vib1 vib2 vib3 vib4 vib5 vib6 vib7 1 85 210 138 162 309 287 176 2 210 85 176 138 162 309 287 3 138 176 85 287 309 162 309 4 162 138 287 85 287 138 162 5 309 162 309 287 85 176 138 6 287 309 162 138 176 85 210 7 176 287 309 162 138 210 85 sweep vib1 vib2 vib3 vib4 vib5 vib6 vib7 1 159 45 196 225 309 249 80 2 80 159 45 196 225 309 249 3 249 80 159 45 196 225 309 4 309 249 80 159 45 196 225 5 225 309 249 80 159 45 196 6 196 225 309 249 80 159 45 7 45 196 225 309 249 80 159 sweep vib1 vib2 vib3 vib4 vib5 vib6 vib7 vib8 1 331 327 101 55 207 67 259 3 2 327 331 3 101 55 207 67 259 3 101 3 331 259 67 55 207 67 4 55 101 259 331 259 101 55 207 5 207 55 67 259 331 3 3 55 6 67 207 55 101 3 331 327 101 7 259 67 207 55 3 327 331 327 8 3 259 67 207 55 101 327 331 sweep vib1 vib2 vib3 vib4 vib5 vib6 vib7 vib8 1 298 122 193 138 247 180 93 78 2 78 298 122 193 138 247 180 93 3 93 78 298 122 193 138 247 180 4 180 93 78 298 122 193 138 247 5 247 180 93 78 298 122 193 138 6 138 247 180 93 78 298 122 193 7 193 138 247 180 93 78 298 122 8 122 193 138 247 180 93 78 298 a) determining the number of vibratory sources; b) determining a number of simultaneous vibration sweeps, the vibratory sources being simultaneously activated during each of said sweeps; c) establishing a phase encoding scheme, wherein phase shifts are calculated using a Monte Carlo simulation; d) positioning said vibratory sources proximate to a predetermined location on the surface of the earth; e) energizing said vibratory sources according to said phase encoding scheme thereby vibrating the earth to produce a seismic signal; f) acquiring the reflected seismic signal data at said one or more receivers; and, g) processing said data to create a seismic survey; the method being characterized in that the phase encoding scheme being a non-orthogonal phase scheme; and step (c) involves choosing from the following pre-defined calculated phase encoding schemes indicating phase shifts, for said number of vibratory sources ("viv1, vib2, vib3, vib4, vib5, vib6, vib7, vib8") and vibration sweeps ("sweep"):